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Abstract 

This  dissertation  presents  a  new  model  for  computing  the  angle  dependent  performance 
measures  of  an  adaptive-optics  system.  By  incorporating  diffraction  caused  by  the  index-of- 
reff  action  variations  of  the  atmosphere,  die  phase  and  amplitude  fluctuations  of  die  propagating 
wave  are  computed.  New  theory  is  presented,  that  uses  the  diffraction-based  propagation 
model  to  yield  optical  transfer  function  (OTF)  expressions  that  are  more  accurate  as  compared 
to  current  theory  that  neglects  diffraction.  An  evaluation  method  far  calculating  the  OTF 
is  presented  that  utilizes  a  layered  atmospheric  model  and  normalized  OTF  expressions. 
The  diffraction  model  is  also  used  to  present  die  first  OTF  signal-to-noise  ratio  (SNR) 
expressions  that  are  a  function  of  separation  angle  between  the  beacon  and  the  object  in 
an  adaptive-optics  system.  An  evaluation  method  fa*  the  SNR  is  presented  that  utilizes 
normalized  correlation  functions  which  are  valid  over  a  wide  range  of  atmospheric  conditions 
and  correction  geometries.  An  analysis  of  die  angle  dependency  of  the  point  spread  function 
(PSF)  is  presented  using  the  derived  OTF  expression.  The  diffraction  model  is  then  used 
to  develop  a  new  adaptive-optics  wavefront  correction  algorithm  that  results  in  an  extended 
correctable  field -of- view  (FOV)  as  compared  to  current  correction  algorithms. 
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A  diffraction-based  model  of 


anisoplanatism  effects  in  adaptive  optic  systems 


/.  Introduction 

1.1  Overview 

The  heating  and  cooling  of  the  earth  along  with  the  mixing  effects  of  the  wind  causes 
die  atmosphere  to  be  non-homogeneous  and  turbulent  The  non-homogeneous  aspect  means 
that  the  atmosphere  at  one  point  in  space  is  different  from  the  atmosphere  at  another  point  in 
space.  In  particular,  temperature  fluctuations  in  the  atmosphere  cause  the  index-of-reffaction 
of  die  atmosphere  to  be  random  in  space  and  time.  The  study  of  how  turbulence  affects 
a  propagating  wavefront  is  well  advanced  with  one  of  die  most  complete  references  being 
written  in  1961  by  Tatarski  (46).  As  die  wavefront  from  an  object,  such  as  a  star,  travels 
through  the  atmosphere,  different  parts  of  die  wavefront  will  travel  through  different  indices 
of  refraction.  This  will  cause  different  parts  of  the  wavefront  to  experience  different  optical 
path  lengths  (20).  The  non-uniform  optical  path  lengths  result  in  bending  of  the  wavefront. 
This  bending  causes  different  portions  of  die  wavefront  to  travel  in  slightly  different  directions 
resulting  in  optical  field  perturbations.  These  optical  field  perturbations  cause  a  significant 
portion  of  the  blur  experienced  when  imaging  through  die  atmosphere  and  is  the  limiting  factor 
in  the  resolution  of  an  image  observed  from  a  large  aperture  ground  based  optical  imaging 
system.  However,  through  recent  advances  in  adaptive-optics,  it  is  now  theoretically  possible 
to  view  an  exo-atmospheric  object  at  resolutions  near  the  diffraction  limit  of  the  imaging 
system  (18). 

Wavefront  correction  through  the  use  of  adaptive  optics  has  been  widely  studied  as 
a  means  to  overcome  turbulence  effects  (1,  2,  17,  18,  19,  34).  Babcock  first  suggested 
‘die  possibility  of  compensating  astronomical  seeing’  in  19S3  (1),  and  die  relatively  recent 
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advances  of  Welsh  and  Gardner  in  the  areas  of  performance  analysis  and  the  use  of  laser  guide 
stars  (54,  55),  may  have  hastened  the  declassification  of  military  involvement  in  adaptive- 
optics  in  May  1991.  Since  that  time,  research  activity  has  increased  and  reports  concerning 
adaptive-optic  applications  have  made  their  way  to  the  open  literature  (8,  31,  45).  Recent 
work  at  AFTT  in  the  areas  of  adaptive-optics  and  image  reconstruction  include  die  work  of 
Johnston  in  the  area  of  multiconjugate  adaptive  optics  (22)  and  the  work  of  Pennington  in 
comparing  the  performance  of  the  shearing  interferometer  and  the  Hartman  wavefront  sensor 
for  extended  sources  and  large  subaperture  spacing  (38).  Also,  Stoudt  investigated  a  technique 
for  improving  the  image  quality  through  a  post  processing  algorithm  that  selected  the  higher 
quality  frames  for  the  image  reconstruction  (44).  VonNiederhausem  developed  an  analysis 
program  for  obtaining  system  transfer  functions  as  a  function  of  number  and  size  of  sensor 
subapertures,  sensor  noise,  and  atmospheric  turbulence  (51).  These  system  transfer  functions 
were  then  used  in  the  image  reconstruction  technique  identified  as  ‘Self-Referenced  Speckle 
Holography’  (53).  In  addition,  Koeffler  researched  a  technique  for  measuring  the  strength  of 
the  turbulence  and  wind  velocity  as  a  function  of  altitude  (26). 

The  theory  of  adaptive  optics  states  that  by  measuring  the  wavefront  phase  of  a  reference 
beacon  such  as  a  star,  a  laser  guide  star  ( 1 3, 54, 55, 56),  or  a  glint  off  of  the  object  to  be  imaged, 
and  then  subtracting  this  phase  from  the  object  wavefront,  an  improvement  in  image  quality 
will  result  Wavefront  measurements  are  made  with  a  wavefront  sensor  (WFS)  such  as  a 
Shack-Hartman  sensor  (18).  These  measurements  are  then  converted  into  deformable  mirror 
control  commands  (52)  which  are  used  to  drive  the  actuators  of  a  deformable  mirror  resulting 
in  the  desired  wavefront  correction.  Theory  indicates  that  perfect  correction  can  be  attained 
with  a  perfect  wavefront  correction  system  when  the  object  wavefront  propagates  along  the 
same  path  as  the  beacon.  This  perfect  correction  is  referred  to  as  diffraction  limited  correction. 

Over  the  past  several  years,  experiments  have  been  performed  using  actual  adaptive- 
optic  systems  (12,  36,  57).  The  results  have  been  encouraging  but  have  never  quite  matched 
the  theoretical  predictions  of  diffraction  limited  correction.  One  of  the  reasons  that  the  ideal 
diffraction  limited  correction  has  been  unattainable  is  that  a  wavefront  correction  system  can 
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never  be  perfect  There  will  always  be  errors  in  the  sampling  and  in  the  reconstruction  of  the 
wavefront  Much  work  is  being  done  to  minimize  these  errors,  but  as  this  analysis  will  show, 
for  a  large  aperture  optical  system  in  a  turbulent  atmosphere,  the  use  of  adaptive-optics  will 
never  allow  the  diffraction  limit  to  be  achieved.  An  understanding  of  the  problems  involved 
is  obtained  by  examining  a  typical  correction  geometry. 

In  the  most  general  case,  the  object  and  beacon  locations  are  separated  by  die  object 
angle,  t),  as  shown  in  Figure  1.  However,  the  performance  of  the  adaptive-optics  system 
is  known  to  degrade  as  d  increases  and  useful  performance  is  limited  to  a  small  range  of 
angles  about  the  beacon  (11, 12).  The  degradation  in  performance  results  from  the  turbulence 
effects  along  one  optical  path  being  only  partially  correlated  with  the  turbulence  in  the  other 
path.  This  partial  correlation  effect  is  referred  to  as  anisoplanatism  (10)  and  the  maximum 
angular  separation  where  good  correlation  is  maintained  is  called  the  isoplanatic  angle  (0O).  In 
addition,  the  effect  of  amplitude  variations  in  the  received  wavefront  degrades  die  performance 
of  an  adaptive-optics  system  which  inherently  only  corrects  for  phase  effects.  Experimental 
results  indicate  that  amplitude  as  well  as  phase  variations  are  a  real  concern  (5). 

Proper  accounting  of  phase  and  amplitude  variations  in  a  received  wavefront  requires 
an  analysis  based  on  diffraction.  Much  of  the  early  work  of  anisoplanatism  justifies  dropping 
diffraction  effects  part  way  through  the  analysis  (10).  This  analysis  approach,  commonly 
referred  to  as  geometric  optics,  is  widely  used  (23,  32,  33,  43,  49),  and  assumes  that  the 
only  effect  die  atmosphere  has  on  a  propagating  wave  is  in  die  form  of  wavefront  phase 
aberrations  due  to  integrated  optical  path  differences.  McKechnie,  for  example,  breaks  the 
atmosphere  into  phase  screens  which  simply  add  a  phase  to  a  propagating  wave  (32).  In 
good  seeing  conditions  where  the  atmospheric  profile  contains  little  high  altitude,  or  far-field 
turbulence,  the  use  of  a  geometric  optics  calculation  is  convenient  and  yields  accurate  results. 
However,  as  conditions  deteriorate  in  the  form  of  a  greater  amount  of  high  altitude  turbulence, 
a  performance  analysis  that  neglects  diffraction  becomes  less  accurate. 

The  analysis  presented  here  considers  an  adaptive-optics  system  which  is  able  to  per¬ 
fectly  measure  the  wavefront  phase  from  a  reference  beacon,  and  in  turn  perfectly  apply  this 
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Object  Direction 


Beacon  Direction 


correction  to  the  object  wavefront  The  constraint  of  an  ideal  adaptive-optics  system,  allows 
the  analysis  to  focus  on  atmospheric  effects,  which  for  any  large  diameter  optics  will  define 
the  upper  bound  of  performance.  A  diffraction  calculation  based  on  die  work  of  Lee  and 
Harp  (30),  and  Tatarski  (46)  is  used  to  account  for  both  phase  and  amplitude  perturbation 
effects  of  a  wavefront  in  the  pupil  plane.  The  diffraction  technique  consists  of  expanding 
the  atmosphere  into  thin  slabs  and  then  resolving  die  refractivity  field  within  each  slab  into 
Fourier  components  of  varying  wave  number.  Each  Fourier  component  is  represented  as  a 
sinusoidally  varying  grating  which  diffracts  the  incident  wave.  The  effect  on  a  wavefront 
from  each  Fourier  component  in  each  slab  can  then  be  calculated  and  summed  in  the  aperture 
to  yield  the  resulting  wavefront  The  adaptive-optics  system  applies  a  phase  correction,  and 
performance  measures  are  calculated  based  on  the  residual  wavefront  phase  and  amplitude 
perturbations  in  the  pupil.  Weak  turbulence  is  assumed,  which  implies  that  the  diffracted 
wave  is  relatively  weak.  Therefore,  die  effects  of  die  diffracted  wave  being  diffracted  again 
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by  subsequent  slabs  is  ignored.  Since  the  adaptive-optics  system  only  corrects  for  phase 
perturbations,  this  diffraction  analysis  allows  a  more  accurate  calculation  of  the  system  per¬ 
formance  than  the  current  theory  that  neglects  diffraction.  The  improvement  in  accuracy  is 
due  to  the  diffraction  analysis  accounting  for  both  phase  and  amplitude  effects.  This  analy¬ 
sis  becomes  essential  for  obtaining  accurate  performance  results  for  adaptive-optic  systems 
operating  in  an  atmosphere  containing  high  altitude  turbulence. 

Since  an  adaptive-optics  system  uses  information  measured  along  one  path  of  thr 
atmosphere  to  correct  a  wavefront  passing  through  a  different  path,  the  analyzed  performano 
measures  will  be  dependent  on  the  angular  separation,  or  ‘object  angle’  between  the  paths. 
One  performance  measure  of  an  adaptive-optics  system  is  die  average  optical  transfer  function 
(OTF).  The  OTF  represents  an  excellent  measure  since  other  measures  such  as  resolution,  the 
point  spread  function  (PSF)  and  die  Strehl  ratio  can  be  derived  from  the  OTF.  The  OTF  is  a 
measure  of  die  attenuation  of  each  object  spatial  frequency.  Knowledge  of  the  PSF  is  useful 
since  in  a  shift  invariant  imaging  system,  the  image  intensity  is  determined  by  convolving  the 
PSF  with  die  object  intensity.  However,  since  an  adaptive-optics  imaging  system  yields  an 
angle  dependent  PSF,  the  system  is  not  shift  invariant  The  image  intensity  in  a  system  that 
is  not  shift  invariant  involves  a  superposition  integral  of  the  shift  dependent  PSF  (in  our  case 
-  angular  dependent  PSF)  and  die  object  intensity  distribution.  The  Strehl  ratio  is  a  useful 
measure  in  adaptive-optics  since  it  defines  how  much  the  peak  of  the  PSF  has  been  reduced 
compared  to  a  diffraction  limited  peak.  An  equally  important  measure  is  the  OTF  signal-to- 
noise  ratio  (SNR).  Knowledge  of  the  SNR  is  particularly  important  when  image  reconstruction 
is  used  in  conjunction  with  adaptive-optics  imaging  (43).  In  image  reconstruction,  die  system 
OTF  can  be  used  in  conjunction  with  the  measured  image  intensity  to  obtain  an  estimate  of 
the  actual  object  intensity.  The  SNR,  which  is  a  measure  of  the  variability  or  randomness  of 
the  OTF  relative  to  the  mean  value  of  the  OTF,  provides  a  measure  of  the  quality  with  which 
the  spectral  components  of  the  OTF  can  be  used  to  reconstruct  the  object  spectrum. 

This  research  effort  is  the  first  to  explicidy  define  a  phase  and  amplitude  OTF  of 
an  adaptive-optics  system  (48)  and  the  first  to  yield  diffraction  based  OTF  SNR  results 
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that  are  a  function  of  the  angle  between  the  reference  beacon  and  the  object  (47).  Fried 
developed  expressions  for  the  angle  dependent  OTF  but  dismissed  diffraction  effects  (9),  and 
Roggemann  (42, 43)  produced  OTF  SNR  results  but  did  not  consider  the  angle  dependency. 

Additionally,  previous  analyses  utilize  the  Kolmogorov  power  spectrum  for  the  atmo¬ 
spheric  index  of  refraction  variations.  The  use  of  the  Kolmogorov  power  spectrum  explicitly 
contains  the  assumption  that  die  ‘outer  scale’,  or  die  largest  eddy  size  of  die  turbulence,  is 
infinite,  creating  a  pole  in  the  power  spectrum  at  the  index  variation  wavenumber,  k  =  0  (14). 
There  is  no  general  agreement  regarding  die  size  of  die  outer  scale,  L0.  Conventional  mea¬ 
surements  put  L0  at  a  value  greater  than  most  telescope  diameters  (5  -  10  meters),  while 
McKechnie  and  others  believe  that  the  value  of  La  is  much  smaller  (20  -  40  cm)  (6, 33).  The 
analysis  presented  here  is  completely  general  with  respect  to  die  refractive  index  power  spec¬ 
trum,  and  allows  for  calculations  of  die  effect  of  finite  inner  and  outer  scales.  In  Chapter  n, 
an  analysis  of  the  effect  of  the  inner  and  outer  scales  on  the  average  OTF  is  presented.  This 
research  effort  produced  the  first  analysis  of  this  type  on  the  inner  and  outer  scale  effects  (48). 

Finally,  a  diffraction  based  analysis  yields  additional  information  regarding  the  prop¬ 
agating  wavefronts.  This  information  is  used  to  develop  a  new  adaptive-optics  wavefront 
correction  algorithm  that  increases  the  correctable  field-of-view  (FOV)  of  an  adaptive-optics 
imaging  system. 

1.2  Problem  statement 

Establish  a  new  model  for  die  theory  and  numerical  evaluation  of  the  perfor¬ 
mance  measures  of  an  adaptive-optics  system  that  is  valid  over  a  wide  range  of 
atmospheric  turbulence  parameters  and  is  a  function  of  separation  angle  between 
the  reference  beacon  and  the  object  source.  The  new  model  is  to  account  for 
diffraction  effects  and  provide  for  an  effective  evaluation  of  the  average  optical 
transfer  function,  the  signal-to-noise  ratio,  and  the  point  spread  function. 

1.3  Contributions 

1.  Development  of  a  new  diffraction-based  model  of  atmospheric  propagation  for  use  in 
calculating  die  angle  dependent  performance  measures  of  an  adaptive-optics  system. 
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The  model  incorporates  diffraction  to  track  both  the  amplitude  and  phase  effects.  An 
evaluation  method  is  presented  that  results  in  rapid  evaluation  of  the  average  OTF 
under  changing  atmospheric  conditions  (48).  The  OTF  analysis  is  then  extended  by 
computing  the  angle  dependent  PSF  through  die  inverse  Fourier  Transform  of  the  OTF. 

2.  Development  of  a  method  for  calculating  the  angle  dependent  OTF  variance  and  SNR 
of  an  adaptive-optics  system.  This  method  uses  the  model  discussed  in  Item  1  above 
to  include  both  amplitude  and  phase  effects.  An  evaluation  method  is  presented  that 
results  in  the  rapid  evaluation  of  die  SNR  under  changing  atmospheric  conditions  as 
well  as  correction  geometries  (47). 

3.  Development  of  a  new  adaptive-optics  wavefront  correction  algorithm  that  extends 
the  correctable  FOV  of  the  adaptive-optic  system.  The  new  model  developed  in  Item 
1  above  yields  additional  information  regarding  die  propagating  wavefronts.  This 
additional  information  is  used  to  develop  the  new  correction  algorithm. 

1.4  Outline  of  Dissertation 

The  remainder  of  this  dissertation  is  organized  as  follows.  Chapter  n  presents  the 
derivation  of  die  object  angle  dependent  average  OTF.  This  derivation  consists  of  developing 
the  wave  perturbation  equations  in  Section  2.2  and  separately  deriving  the  expression  for  die 
amplitude  and  phase  OTFs  in  Sections  2.3  and  2.4  respectively.  Chapter  n  also  presents 
simplifications  necessary  to  rapidly  evaluate  the  OTF’s  under  changing  atmospheric  condi¬ 
tions  (48).  Chapter  HI  presents  an  analysis  of  the  object  angle  dependent  variance  and  SNR 
of  the  OTF,  along  with  the  required  simplifications  and  method  of  evaluation  (47).  Chap¬ 
ter  IV  presents  an  analysis  of  the  the  PSF  and  the  Strehl  ratio  utilizing  die  OTF  developed 
in  Chapter  IL  Chapter  V  considers  a  new  wavefront  correction  algorithm  that  has  the  effect 
of  extending  the  correctable  FOV  of  the  adaptive-optics  system.  This  correction  algorithm  is 
a  natural  utilization  of  the  wave  perturbation  equations  developed  in  Section  2.2.  Finally,  a 
summary  and  conclusions  of  this  research  effort  are  presented  in  Chapter  VI 
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II.  The  angle  dependent  optical  transfer  function  ( OTF ) 


2.1  Introduction 

The  purpose  of  this  chapter  is  to  develop  a  model  that  allows  for  the  derivation  and 
evaluation  of  an  angle  dependent  average  optical  transfer  function  (OTF)  for  an  ideal  adaptive- 
optics  system  (48).  The  model  is  to  be  valid  over  a  wide  range  of  atmospheric  conditions 
and  is  to  incorporate  diffraction  to  properly  track  phase  and  amplitude  effects  of  the  corrected 
wavefront  in  the  pupil  plane  of  the  adaptive-optics  system.  An  understanding  of  die  modeling 
of  atmospheric  turbulence  as  well  as  the  use  of  phase  screens  is  assumed.  A  model  of 
atmospheric  turbulence  is  presented  in  Appendix  A  and  the  concept  of  phase  screens  is 
discussed  in  Appendix  B.  The  analysis  begins  by  considering  a  linear  imaging  system. 

In  a  linear  imaging  system,  the  object  intensity  may  be  considered  to  be  composed 
of  a  sum  of  the  object  spatial  frequencies.  The  response  of  a  linear  system  to  each  spatial 
frequency  component  can  be  calculated  separately.  An  image  is  formed  by  summing  each 
spatial  frequency  at  die  output  of  the  system.  The  OTF  is  a  measure  of  the  attenuation  of  each 
object  spatial  frequency  as  it  passes  through  the  system.  The  linear  system  of  interest  in  this 
research  is  the  atmosphere  along  with  the  adaptive-optics  imaging  system. 

The  atmosphere  imposes  optical  distortions  on  a  propagating  wave  due  to  the  index- 
of-refr action  variations.  These  variations  are  random,  thereby  producing  an  optical  wave 
that  must  be  considered  a  random  process.  The  OTF  that  models  the  imaging  system  will 
also  be  a  random  process  where  the  quantity  of  interest  is  die  average  value  of  the  OTF.  We 
distinguish  between  die  phase  and  amplitude  variations  of  the  corrected  field  in  the  pupil 
plane  of  the  adaptive-optics  system  by  defining  Hp  as  the  average  OTF  due  to  die  residual 
phase  of  the  corrected  wavefront  in  die  pupil,  and  Ha  as  the  average  OTF  due  to  the  amplitude 
variations  of  the  corrected  wavefront  in  the  pupil  Based  on  the  workofTatarski  (46),  Fried  has 
shown  that  phase  and  amplitude  variations  of  a  wavefront,  resulting  from  propagating  through 
atmospheric  turbulence,  may  be  considered  as  independent  wide  sense  stationary  Gaussian 
random  variables  (9).  Therefore,  Ha  and  Hp  can  be  calculated  separately  and  multiplied 
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together  to  form  the  overall  OTF  due  to  die  atmosphere.  The  OTF  of  the  optics  will  not  be 
addressed  in  this  chapter,  but  the  atmospheric  OTF  can  always  be  multiplied  by  the  optics 
OTF  to  yield  an  overall  system  OTF. 

We  begin  by  calculating  the  wavefronts  from  a  beacon  and  an  object  that  have  propagated 
through  die  atmosphere  at  an  angle  t?  as  shown  in  Figure  1  on  page  4.  This  calculation  includes 
diffraction  which  tracks  the  phase  and  amplitude  variations  in  the  propagating  wave.  The 
adaptive-optics  system  subtracts  die  phase  of  die  beacon  field  from  the  phase  of  the  object 
field.  The  calculations  of  Ha  and  Hp  are  then  based  on  the  correlation  of  the  amplitude  and 
phase  variations  of  the  residual  field  in  the  pupil  plane.  The  OTF  curves  calculated  here  will 
be  expressed  in  terms  of  the  shift  variable  p  located  in  the  aperture  of  the  system.  Once  the 
optics  of  the  system  are  chosen,  die  shift  variable  can  be  related  to  spatial  frequency  v  by 
v  =  fj,  where  /  is  the  focal  length  of  the  optics  and  A  is  the  optical  wavelength. 

The  remainder  of  this  chapter  is  organized  as  follows.  In  Section  2.2.  expressions  are 
presented  for  die  amplitude  and  phase  perturbations  of  a  wave  that  has  propagated  through 
atmospheric  turbulence.  Much  of  the  actual  derivation  of  these  equations  is  deferred  to 
Appendix  C.  The  amplitude  and  phase  perturbation  equations  are  then  used  to  derive  the 
expressions  for  Ha  and  Hp  in  Sections  2.3  and  2.4,  respectively.  The  derived  OTF  expressions 
are  kept  general  with  respect  to  the  refractive  index  power  spectrum.  These  generalized 
expressions  allow  for  die  analysis  of  the  effect  on  the  OTF  due  to  die  scale  sizes  of  die 
turbulence  as  presented  in  Section  2.3.  In  Section  2.6,  a  new  method  of  modeling  the 
atmospheric  turbulence  through  the  use  of  discrete  layers  at  preselected  altitudes  is  presented. 
The  method  of  preselected  altitudes  is  essential  in  obtaining  an  OTF  that  can  be  scaled  to 
different  atmospheric  conditions.  Section  2.7  introduces  die  normalized  OTF  that  is  shown 
to  be  scalable  over  a  large  range  of  atmospheric  conditions.  Numerical  results  using  the 
normalized  OTF  are  then  presented  in  Section  2.8.  In  Section  2.9,  a  direct  comparison  is 
made  between  the  diffraction  method  of  OTF  calculation,  as  presented  in  this  chapter,  and  the 
geometric  optics  method  of  OTF  calculation.  Finally,  Section  2.10  summarizes  the  chapter. 
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2.2  Phase  and  amplitude  perturbation  equations 

This  section  develops  the  necessary  equations  to  describe  die  amplitude  and  phase 
variations  of  a  wave  that  has  propagated  through  the  atmosphere  at  some  angle  -9  with  respect 
to  the  reference  direction  (see  Figure  1  on  page  4).  We  begin  by  deriving  the  wavefront 
perturbations  resulting  from  a  single  Fourier  component  of  the  refr activity  field  in  a  single 
thin  slab  of  the  atmosphere  located  at  an  altitude  r?.  We  will  then  integrate  over  all  Fourier 
components  and  all  altitudes  to  obtain  the  total  wavefront  perturbations.  A  single  Fourier 
component  of  spatial  frequency,  k  ,  represents  a  sinusoidal  refractive  index  variation.  This 
single  frequency  component  slab  located  at  an  altitude  r?,  imposes  a  phase  variation  on  the 
incident  wavefront  given  by  A(tj,  k)  cos (k  •  x  +  (f>0(y, «)),  where  A(rj,  k)  and  <f>0(v, *)  are 
related  to  the  amplitude  and  phase  of  the  refractivity  field  Fourier  component  as  shown  in 
Appendix  B,  and  x  is  a  position  vector  in  the  plane  of  the  slab.  Propagating  the  perturbed 
wavefront  to  die  pupil  plane  creates  a  complex  wavefront  which  is  fully  described  by  its 
amplitude  and  phase  components. 

The  process  of  determining  die  amplitude  and  phase  variations  in  a  propagating  wave 
is  accomplished  using  the  approach  of  Lee  and  Harp  (30)  with  the  complete  derivation  given 
in  Appendix  C.  The  amplitude  variation  is  given  by 

Pa(?7,K,x,i9)  =  A(tj,  k)  sin(i;a)  cos(r/0  —  (k  •  x  +  <t>0(r], «))),  (1) 

and  the  phase  variation  is  given  by 

Pp(r},K,x,ti)  =  - A(r},  it)  cos(qa)  cos{tj -  (it  ■  x  +  <t>0(r),  it))),  (2) 

where,  a  =  0  =  k  ■  and  k  is  the  wavenumber  of  the  propagating  wave.  The  angle 

dependence  of  the  amplitude  and  phase  perturbations  is  given  as  t)  in  the  left  hand  side  of 
Eqs.  (1)  and  (2).  This  is  a  valid  representation  since,  as  shown  in  Appendix  C,  the  orientation 
of  t)  defines  the  x  axis  of  die  system  (Le.,  t?  =  fix).  The  amplitude  and  phase  perturbations 
of  Eqs.  (1)  and  (2),  contain  two  differences  from  those  given  by  Lee  and  Harp.  First  an 
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approximation  is  used  in  place  of  Lee  and  Harp’s  expression  for  the  wavenumber  of  the 
scattered  wave,  (\A2  -  k2  -  k).  Using  die  binomial  expansion  and  die  assumption  that 
the  wavenumber  of  die  scattered  wave  is  given  by 


where,  k  =  |k|.  Second,  Lee  and  Harp  only  considered  on-axis  propagation.  With  on- 
axis  propagation,  all  orientations  of  it  will  scatter  the  incident  wave  at  angles  with  the  same 
magnitude.  The  analysis  presented  here  considers  off-axis  propagation  where  the  orientation 
of  k  with  respect  to  die  direction  of  incident  propagation  now  determines  the  direction  of  the 
scattered  wave.  An  off-axis  propagation  analysis  produces  the  rj/3  term  of  Eqs.  (1)  and  (2), 
where  rj0  =  tjk  ■  x).  The  angle  dependence  of  Eqs.  (1)  and  (2)  yields  the  angle  dependence  of 
the  performance  measures  discussed  in  this  dissertation. 

The  idealized  adaptive-optics  system  subtracts  the  measured  phase  of  the  beacon  from 
the  phase  of  the  object  The  residual  phase  of  die  adaptive-optics  system  for  this  single  slab, 
designated  A Pp(r/,  it,  z,  t?),  is  the  difference  between  Pp{r),  k,  x,  )  and  Pp(r},  it,  x,  0): 

APp(»?,/c,z,i?)  =  —A(tj,  k)  cos(rja)  cos(i]f3  —  («  •  z  +  <t>0(ri,  k))) 

+A(t],  k)  cos (rja)  cos («  •  z  +  (f>0(r],  it)).  (4) 


The  total  residual  phase  is  found  by  integrating  over  all  altitudes  tj  and  all  frequencies  it. 

APp(z,i?)  =  -  J J  drjdit^A(T],K)cos(T]a)cos(T]/3  -  (it  •  x  +  <f>0(ri ,  it))) 

-A(t),  k)  cos(»7a)  cos(ic  •  z  +  <j>0(ri,  «))^ .  (5) 

This  analysis  is  simplified  by  use  of  complex  phasor  notation  where  the  amplitude  and  residual 
phase  variations  can  be  written  as 

Pa(x,  i?)  =  —  V,  jyy  6ij6kA(t),  sin(77a)e^| ,  (6) 
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and, 


A Pp(£,  t?)  =  ^  jyy  di/d/cA(^,K)e  cos(qa)(l  -  eJ^)| ,  (7) 

where,  ^(t/,  «)  =  A(»?,  /e)e_J*o(,,,*)  and  H  is  an  operator  that  takes  the  real  portion  of  a 
complex  quantity.  The  H  operator  is  now  dropped  and  work  continues  with  the  complex 
phasor  representations  Pa(x,  i9)  and  A Pp(x,  t?)  given  by 

Pa(x,  t?)  =  //d^(,,/c)e-^sin(,a)e^,  (8) 

and, 

APp(x,  i?)  =  J j 6t)<1kA(ti,  K)e~^x'^  cos(r/o)(l  —  (9) 

23  Amplitude  OTF  (Ha) 

The  derivation  of  i/a  treats  amplitude  effects  radio*  than  log-amplitude  effects  as  given 
in  other  research  (9,  14).  The  log-amplitude  can  be  a  useful  quantity  and  is  used  by  many 
researchers  in  the  study  of  atmospheric  propagation  (4,  9,  7,  28).  By  assuming  Gaussian 
statistics  and  small  amplitude  variations,  the  log-amplitude  leads  to  die  wave  structure  function 
as  discussed  in  (14).  However,  the  assumption  of  small  amplitude  variations  becomes  less 
valid  as  high  altitude  turbulence  strength  increases.  Since  the  analysis  of  the  amplitude  OTF 
will  not  rely  on  the  wave  structure  function,  the  final  result  of  this  section,  Eq.  (18),  will  be 
used  as  a  more  general  amplitude  OTF  expression  which  is  valid  under  stronger  turbulence 
conditions  than  a  log-amplitude  result 

The  derivation  of  Ha  begins  by  representing  the  amplitude  variations  of  the  wavefront 
in  the  pupil  plane  as  a  wide-sense  stationary  amplitude  transmittance  function  given  by 

tg(x,d)  =  t0  +  Pa(x,d),  (10) 
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where  in  general.  ta  is  some  real  and  nonnegative  bias  lying  between  zero  and  unity.  Here, 
we  assume  tQ  =  1.  The  average  OTF  of  this  transmittance  function  is  then  given  by  (14) 


Ha(p) 


1  +  Tg(p  ) 

i  +  r.io)  ’ 


(ii) 


where  p  is  a  shift  vector  in  die  pupil  plane  and  T a(p  )  =  (Pa(x,d)Pa(x  +  p,  t))).  The 
amplitude  correlation  function  is  related  to  the  complex  phasor  representation  of  the  amplitude 
perturbations  given  in  Eq.  (8)  by 


Up)  =  iR{r.(p)}  =  <)•£•(* -?,»)>•  (12) 


The  notation  ()  represents  the  ensemble  average.  Here  we  use  die  property  that  die  correlation 
of  die  complex  representation  is  twice  the  value  of  the  real  valued  correlation  (14).  Substituting 
Eq.  (8)  into  Eq.  (12)  gives 


Up)  =  j j j j 

x  sin(?7ia)sin(?/2a)e-7^T'I“/J,'2\  (13) 


The  relations  given  in  Appendix  B  and  Appendix  D  are  now  used  for  propagation  through 
long  distances  of  atmospheric  turbulence  to  yield 

(A(7/i,«i)i4*(j72,«2))  =  8t rk26(Ki  -  k2)6(th  -  (14) 

where,  4>(k,  tj)  is  die  three  dimensional  refractive  index  power  spectrum  with  k2  =  k\  +  «2, 
and  6  is  die  Dirac  delta  function.  The  6(7/2  -  hi)  part  of  Eq.  (14)  is  a  derived  relation  and 
should  not  be  considered  as  a  statement  of  independence  of  turbulent  layers.  Substituting 
Eq.  (14)  into  Eq.  (13)  gives 

Ta(p)  =  Snk2  J j  d7/dK$(K,7/)sin2(r/a)e-j('7'p).  (15) 
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Note  that  Eq.  (IS)  is  not  a  function  of  either  or  t?.  Taking  half  the  real  portion  of  Eq.  (IS), 
the  final  form  of  the  amplitude  correlation  function  is  given  by 

ra(p  )  =  4nk2  J  j  dr}dii$(K,  tj)  sin 2{rja)  cos(/c  •  p  ).  (16) 


Substituting  Eq.  (16)  into  Eq.  (11)  yields  an  expression  for  the  average  amplitude  OTP, 


Ha(p) 


1  -f  4irk2ff  dT)dK$(ii,  tj )  sin^a)  cos (k  •  p  ) 
1  +4irk2ff  drjdK$(K,T})  sin2 (r/a) 


and,  substituting  in  the  expression  for  a,  where  a  =  we  are  left  with 


Ha(p) 


1  +  4irk2ff  drjdK${it,  tj)  sin2  cos(/c  -  p  ) 

1  +  4 irk2ff  drjdK$(K,  tj)  sin2  (»7^) 


(17) 


(18) 


Equation  (18)  represents  the  final  result  of  this  section.  An  example  of  an  amplitude  OTF 
is  shown  in  Figure  2.  This  figure  is  simply  an  evaluation  of  Eq.  (18)  where  an  optical 
wavelength  of  0.5  microns  and  a  Hufhagel- Valley  turbulence  profile  with  a  54  mph  upper 
altitude  wind  (21)  has  been  assumed.  This  turbulence  profile  is  labeled  HV-54  in  Figure  5.  In 
this  calculation,  the  Von  Karman  power  spectrum  is  used,  as  discussed  in  Appendix  A,  with 
the  inner  scale,  Lm,  set  at  1  millimeter  and  the  outer  scale,  Lot  set  at  5  meters.  The  use  of 
turbulence  profiles  is  discussed  in  more  detail  in  Section  2.6  and  an  analysis  of  the  inner  and 
outer  scale  effects  is  presented  in  Section  2.5.  Here,  it  is  important  only  to  note  the  initial 
sharp  roll  off  and  then  the  constant  value  attained  by  the  transfer  function.  This  is  consistent 
with  the  amplitude  transmittance  functions  given  by  Goodman  (14)  and  will  be  characteristic 
of  all  amplitude  OTFs. 


2.4  Residual  phase  OTF  ( Hp) 

The  analysis  of  the  average  phase  OTF,  (Hp),  begins  much  the  same  as  the  last  section. 
Here  file  result  for  the  total  residual  phase  given  by  Eq  (5),  and  file  corresponding  complex 
phasor  representation  given  by  Eq.  (8),  is  used  to  define  a  phase  only  transmittance  func- 
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Figure  2.  Amplitude  OTF  calculated  for  a  wavelength  of  0.5  microns,  a  Hufnagel- Valley 
turbulence  profile  with  a  54  mph  upper  altitude  wind  (21),  and  a  Von  Karman 
power  spectrum  with  Lm  =  1mm  and  L0  =  5m. 


tion  (14).  The  residual  phase  is  considered  as  a  single  phase  screen  with  a  transmittance  of 
tp(x,  d)  =  exp(j  APp(x,  j?)).  The  OTF  of  file  phase  screen  is  then  calculated  using  (14) 


HP(pJ) 


(tp(x,#)t;(x-p,#)) 

(exp(jAPp(x,  0))  exp(—jAPp(x  -  p  ,  i?))) 

(exp(;  APp(x,  i?))  exp(->  APp(x,  i?))) 
(exp(j(APp(x,  )  -  APp(x  -  p  ,  d)))) 

1 


Assuming  that  the  residual  phase  obeys  Gaussian  statistics  gives, 

Hp(p  ,  *)  =  exp  { -|<( APp(x,  tf )  -  APp(x  -  p  ,  tf))2) }  . 


(19) 


(20) 
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Expanding  the  square  within  the  expected  value,  and  using  the  assumption  that  the  residual 
phase  is  wide-sense  stationary,  the  following  expression  for  the  residual  phase  OTF  is  obtained: 


=  exp{r,(p  ,«)  -  r,(M)}.  (21) 

As  with  the  amplitude  correlation  function,  the  complex  phasor  representation  of  die  phase 
correlation  function  is  defined  as 

£p(/M)  =  (AP^x^AP^ix-pJ)).  (22) 

Substituting  Eq.  (9)  into  Eq.  (22)  gives 

Tp(p,fl)  =  ////  dVidKidr)2dK2{A(T)UK1)A*{Ti2,K2))e~3{*1'*~*2'{*~p)) 

x  cos(r]ia)(l  —  e^m)cos(rj2a)(l  —  e~J0Tt2).  (23) 

Using  the  same  independence  relation  given  by  Eq.  (14),  yields 

Tp(pJ)  =  8rk2  J J  dT]dK$(K,  T])e~^x  ^^  cos2(r)a)(l  —  e^r,)(l  —  e~i0r)) 

=  8irk2  JJ  d?/d«$(«,  77)  cos2(i/a)|l  —  e^r,\2e~^'^ .  (24) 

For  the  final  form  of  the  correlation  function,  half  the  real  part  of  Tp(p  ,  t?)  is  considered, 

Tp(p  ,1))  =  4nk2  J J  d?;dK$(«,  i;)cos(k  •  ^)cos2(j;a)(2  —  2cos(t7/?)),  (25) 

and,  substituting  in  the  expressions  for  a  and  /?,  where  a  =  0  =  k  ■  d,  results  in 

Tp(p  ,d)  =  8nk2  J  J  cos2(-^-)  cos(k  •  p)[l  _  cos(t]k  ■  d)]dt)dii.  (26) 


Substituting  Eq.  (26)  into  Eq.  (21),  yields  the  final  expression  for  Hp: 


Hp(p  ,i?)  =  exp  8tt k2  J  J  cos2(^^-)[1  —  cos(/c  •  p)][l  —  cos(r?/c  •  d)]drfdii 

(27) 

This  is  die  OTF  due  to  die  residual  phase  of  an  ideal  adaptive-optics  system  where  the  beacon 
and  the  object  are  separated  by  an  angle  t?.  Note  that  this  OTF  is  not  circularly  symmetric 
about  the  zero  frequency.  An  example  of  a  phase  OTF  is  shown  in  Figure  3.  Although  the 
OTF  equations  allow  for  complete  two  dimensional  calculations,  the  results  of  Figure  3  are 
calculated  for  p  along  die  axis  parallel  to  d  which  defines  the  x  axis  of  the  system.  This  figure 
is  an  evaluation  of  Eq.  (27)  with  an  optical  wavelength  of  O.S  microns  and  a  Hufhagel- Valley 
turbulence  profile  with  a  54  mph  wind  (21)  as  discussed  in  Section  2.6.  The  object  angle,  d, 
equals  2.4  /jrad,  which  is  equal  to  the  isoplan ahc  angle  as  defined  by  Fried  (10).  Once  again, 
the  Von  Karman  power  spectrum  with  die  inner  scale,  Lm,  set  at  1  millimeter  and  die  outer 
scale,  L0,  set  at  5  meters  is  used.  The  effect  of  inner  and  outer  scale  on  Ilv  is  discussed  in 
section  2.5.  Here,  note  die  characteristic  shape  of  a  phase  OTF.  The  initial  sharp  roll  off  is 
followed  by  a  much  more  gradually  decreasing  value  of  the  OTF. 

25  Effect  of  inner  and  outer  turbulence  scale  sizes  on  the  OTF 

The  purpose  of  this  section  is  to  quantify  die  significance  of  the  turbulence  scale  sizes 
on  die  OTF.  As  previously  stated,  this  analysis  refrains  from  taking  a  stand  on  the  value  of 
the  turbulence  scale  sizes  of  the  atmosphere,  but  rather  presents  a  method  of  calculating  the 
effect  Up  to  this  point  a  specific  form  for  the  turbulence  power  spectrum  $(£)  has  not 
been  assumed.  The  power  spectrum  can  be  written  as  a  constant  times  a  normalized  power 
spectrum,  $„(«).  Details  of  the  constant  are  discussed  in  Section  2.6.  Most  current  research 
utilizes  the  Kolmogorov  power  spectrum  where  $0(«)  =  |«|-n/3.  The  Kolmogorov  power 
spectrum  yields  analytic  results  to  many  problems.  Also,  this  form  is  known  to  match  well  with 
actual  measured  data  over  die  region  of  the  inertial  subrange  which  is  the  range  between  the 
inner  and  outer  scale  values  of  |  k  | .  What  sometimes  becomes  lost  with  use  of  the  Kolmogorov 
spectrum  is  die  implication  that  the  inner  scale  equals  zero  and  the  outer  scale  is  infinity.  A 
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Figure  3.  A  phase  OTF  calculated  for  a  wavelength  of  O.S  microns,  a  separation  angle 
of  2.4  //rad,  a  Hufhagel- Valley  turbulence  profile  with  a  34  mph  high  altitude 
wind  (21),  and  a  Von  Karman  refractive  index  power  spectrum  with  Lm  =  1mm 
and  L0  =  5  m. 


non-zero  inner  scale  and  a  finite  outer  scale  is  incorporated  into  the  calculations  by  using  a 
modification  to  die  Kolmogorov  spectrum  known  as  the  Von  Karman  power  spectrum  given 
by 


*.(*>  =  7T7 


(1*1 2  +  $>n/6’ 


(28) 


where  Lm  and  La  represent  the  inner  and  outer  turbulence  scale  sizes  respectively.  To 
determine  die  effect  of  the  parameter  L0,  the  OTF  value  is  evaluated  for  various  values  of 
L0  at  a  single  spatial  frequency.  As  noted  before,  this  spatial  frequency  is  related  to  the  shift 
parameter  p  by  p  =  v\f.  Figure  4  shows  die  effect  of  die  outer  scale  size  on  the  phase  and 
amplitude  OTFs  for  a  typical  Hufiiagel- Valley  atmosphere  (21)  evaluated  at  \p\  =  1  m  and 
oriented  parallel  to  t?.  From  Figure  4,  it  is  seen  that  the  effect  on  the  phase  OTF  caused  by 
differences  in  outer  scale  dimension  are  significant  enough  to  warrant  concern  about  the  proper 
selection  of  L0.  hi  a  similar  manner,  the  effect  of  the  inner  scale  on  the  OTF  was  studied.  The 
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inner  scale  was  found  to  cause  insignificant  changes  in  the  OTF  over  all  conceivable  scale 
sizes.  The  remainder  of  this  work  utilizes  scale  sizes  of  La  =  5  m  and  Lm  =  1  mm. 


Figure  4.  //o(1.0)  and  Hp(  1.0)  versus  La.  This  figure  gives  a  measure  of  outer  scale 
significance  on  the  calculated  OTF  values.  The  OTFs  are  evaluated  for  p  =  1 
meter,  i?  =  2/rrad,  and  a  Hufnagel- Valley  atmosphere  (21)  with  a  54  mph  upper 
atmosphere  wind. 


2.6  Discrete  atmospheric  layers  at  preselected  altitudes 

Although  Eqs.  (18)  and  (27)  are  useful,  they  are  time  consuming  to  compute  and  do 
not  lend  a  great  deal  of  insight  to  the  parameter  dependences.  The  purpose  of  this  section 
is  to  begin  simplifying  the  OTF  expressions  by  creating  a  layered  atmosphere  approximation 
to  the  continuous  atmosphere.  The  result  of  this  process  will  serve  two  purposes.  First,  by 
simplifying  the  OTF  with  respect  to  the  atmospheric  turbulence  profile,  an  OTF  is  obtained 
that  can  be  scaled  with  changing  atmospheric  parameters  rather  than  completely  recalculated. 
Second,  notice  that  the  integration  within  die  OTFs  of  Eqs.  (18)  and  (27)  are  over  the  altitude 
parameter  r),  and  tire  two  dimensional  frequency  vector  k.  By  creating  a  layered  atmospheric 
model,  the  integration  over  rj  can  be  replaced  with  a  summation  over  the  layers.  Layering 
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die  atmosphere  can  reduce  the  computation  time  by  over  90%  and  still  yield  results  that  are 
within  1%  of  that  obtained  using  a  continuous  atmosphere. 

The  subject  of  layering  the  atmosphere  begins  by  considering  the  turbulence  power 
spectrum,  $(/c,  77).  The  power  spectrum  contains  information  regarding  the  frequency  dis¬ 
tribution  of  die  index  of  refraction  variations  and  die  distribution  of  turbulence  strength  as  a 
function  of  altitude.  These  two  functional  dependencies  are  commonly  broken  up  as  a  turbu¬ 
lence  strength  profile  represented  by  C*(rj)  and  a  normalized  power  spectral  density  $0(£) 
as  (46) 

$(«,»?)  =  0.033Cn2(J?)$o(«).  (29) 

Several  C2(r/)  turbulence  profiles  are  shown  in  Figure  5  where  C*(rj)  is  given  in  units 
of  meters-2^3.  The  goal  is  to  determine  a  discrete  model  representation  of  C2( tj),  where 
properly  weighted  turbulent  layers  are  placed  at  specified  altitudes. 

Creating  discrete  turbulence  layers  is  a  standard  function  approximation  problem  where 
die  function  to  be  approximated  is  die  C2(r?)  turbulence  profile.  The  only  issue  is  deciding 
when  a  good  approximation  has  been  obtained.  For  the  diffraction  calculations  presented  in 
this  research,  higher  order  moments  of  the  layered  model  of  the  atmospheric  turbulence  must 
match  the  higher  order  moments  of  the  continuous  turbulence  profile.  For  example,  the  first 
moment,  called  the  center  of  mass,  could  be  matched  by  either  placing  a  single  layer  at  exactly 
the  center  of  mass  of  the  continuous  C^(r} )  profile,  or  by  properly  weighted  layers  somewhere 
on  either  side  of  die  known  center  of  mass.  Higher  order  moments  can  be  similarly  matched. 
Two  general  techniques  exist  for  accomplishing  this  type  of  approximation.  The  first,  which 
utilizes  Gaussian  Quadrature  (39),  specifies  both  die  location  and  the  strength  of  the  layers 
and  for  n  layers  will  match  2n  moments  (die  integrated  value,  or  zeroth  moment,  plus  the  first 
2n  -  1  moments).  For  the  case  where  a  C2(7?)  turbulence  profile  is  known,  and  a  layered 
equivalent  atmosphere  is  desired,  this  is  the  best  method  available.  This  method  was  used  to 
model  many  turbulence  profiles,  and  it  was  found  that  by  using  only  two  discrete  layers,  OTFs 
were  obtained  that  were  within  1%  of  the  OTFs  calculated  using  a  continuous  C2(r/)  profile. 
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Altitude 


Figure  5.  C^(r/)  Profiles:  Model  SLC-N  represents  a  fit  to  the  ARPA  Maui  Optical  Station 
(AMOS)  night  data  (IS,  35),  and  SLC-D  represents  a  fit  to  the  AMOS  night  data 
with  the  addition  of  an  altered  boundary  layer  to  simulate  daytime  conditions  (15, 
35).  Model  Greenwood  is  Greenwood’s  ‘good  seeing’  model  (15),  and  TRW  is 
a  TRW  high  turbulence  model  for  the  Capistrano  Test  Site  (CTS)  environment. 
HV-21  and  HV-54  represent  the  Hufhagel- Valley  model  calculated  with  a  21  and 
54  mph  upper  atmospheric  wind,  respectively  (21). 
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One  problem  with  using  the  Gaussian  Quadrature  method  is  that  for  a  new  C*  turbulence 
profile,  the  algorithm  must  be  reaccomplished  to  determine  new  layer  strengths  and  altitudes. 
In  order  to  obtain  the  full  2 n  matched  moments,  Gaussian  Quadrature  requires  the  extra  degree 
of  freedom  of  determining  the  altitude  of  die  layers.  This  process  is  accurate  but  it  does  not 
allow  the  transfer  function  to  be  scaled  for  different  atmospheric  conditions.  To  accomplish 
the  goal  of  a  scalable  transfer  function,  the  altitudes  of  die  turbulent  layers  must  be  fixed 
and  then  weights  separately  determined  that  will  properly  approximate  the  desired  continuous 
CI(tj)  turbulence  profile.  This  is  accomplished  by  picking  die  altitudes  and  then  solving  die 
matrix  equation 
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where  //,  represents  the  altitude  of  the  ith  layer,  is  the  weight  (or  percent  of  die  integrated 
Cl  turbulence  profile)  for  the  ith  layer,  and  m,  is  the  normalized  ith  moment  of  die  turbulence 
profile.  This  method  allows  us  to  specify  the  altitudes  of  the  layers,  but  n  layers  only  match  die 
first  n  moments  (the  integrated  value,  or  zeroth  moment,  plus  die  first  n  —  1  moments).  There 
are  also  constraints  placed  on  the  altitudes  so  as  to  ensure  positive  weights,  but  a  discussion 
of  these  constraints  are  beyond  die  scope  of  this  research.  Here,  it  is  sufficient  to  state  that  the 
choice  of  altitudes  is  flexible  but  not  completely  arbitrary.  Using  this  method,  it  was  found 
that  the  proper  placing  and  weighting  of  four  turbulent  layers  (altitudes  of  200  m,  2  km,  10 
km,  and  18  km)  could  model  most  standard  turbulence  models  and  also  yield  OTFs  that  were 
within  1%  of  the  OTFs  calculated  using  a  continuous  atmosphere.  The  calculated  weights  for 
several  different  models  are  shown  in  Table  1  with  die  corresponding  profiles  being  given  in 
Figure  5. 
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Table  1 .  Weights  for  a  4-layered  atmospheric  model.  Layer  1  is  placed  at  200  meters,  layer  2 
at  2  km,  layer  3  at  10  km,  and  layer  4  at  1 8  km.  W,  represents  the  weight  of  layer  t  as 
a  percentage  of  the  integrated  turbulence.  Cl.  Model  SLC-N  represents  a  fit  to  the 
AMOS  night  data  (IS,  35),  and  SLC-D  represents  a  fit  to  the  AMOS  night  data  with 
the  addition  of  an  altered  boundary  layer  to  simulate  daytime  conditions  (IS,  35). 
Model  Greenwood  is  Greenwood’s  ‘good  seeing’  model  (IS),  and  TRW  is  a  TRW 
high  turbulence  model  for  the  CTS  environment  HV-21  and  HV-54  represent  the 
Hufnagel- Valley  model  calculated  with  a  21  and  54  mph  upper  atmospheric  wind, 
respectively  (21). 
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43 

The  discrete  layers  allow  file  integral  over  rj  in  Eq.  (27)  to  be  replaced  with  a  summation 
over  the  layers.  This  leaves  a  phase  OTF  expression  of 

{m  r  n/f2 

-8tt0.033A:2^C^(»7,)  J$0{k)  cos2(-^-) 

X  [1  —  COs(k  •  ^)][1  —  COS (f]iK  •  t9)]d/c  j 

=  IIexp{  ~  8tt0.033A;2C2(77,)  J $„(k)cos2(^-) 

x  [1  —  cos(k  •  p)][l  -  cos(rjiK  ■  t?)]d«j,  (31) 

where  Cl(r}t)  represents  the  turbulence  strength  of  the  ith  layer  located  at  altitude  r?t,  and  m 
is  the  number  of  turbulent  layer.  Cl(r),)  can  be  replaced  by  a  percentage  of  the  integrated 
turbulence  W,C2: 

CfM  =  W,  J  ClWdri  =  W,Cl  (32) 

where  the  integration  is  over  the  path  from  the  source  to  the  aperture  plane.  In  the  case  of  a 
source  outside  the  atmosphere  and  an  aperture  on  file  ground,  file  integration  is  over  the  entire 
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vertical  C2  profile.  This  allows  the  phase  OTF  to  be  expressed  as 


Hp(p,0)  =  nexpj-  SvO.OZWWiClj  $0(«)cos2(^) 
x  [1  —  cos(k  •  p)][l  —  cos(i;,7C  •  d)]<f/cj. 


(33) 


Following  a  similar  analysis,  the  amplitude  OTF,  Ha,  can  be  written  as 

ff  r  1  +  47r0.033fc2C2  ££ i  WJ  *0{Z)  sin 2(Via)  cos (k  •  p)dH 
a{P)  1  +  47r0.033ifc2C2  YZi  sin2(r/,a)d« 


(34) 


2.7  Simplified  OTF  calculations 

Equations  (33)  and  (34)  represent  die  complete  phase  and  amplitude  OTFs  for  a  layered 
atmosphere.  However,  they  still  require  recalculation  if  the  integrated  turbulence  or  the 
distribution  of  turbulence  changes.  In  this  section,  Eqs.  (33)  and  (34)  are  simplified  by 
incorporating  Fried’s  coherence  parameter  (9),  rot  and  introducing  die  concept  of  a  scalable 
transfer  function.  This  will  allow  OTFs  to  be  quickly  computed  undo-  changing  atmospheric 
conditions. 

Fried’s  coherence  parameter  is  commonly  used  to  describe  the  ‘seeing’  conditions  of 
the  atmosphere.  It  is  known  that  die  OTF  of  the  atmosphere  is  roughly  equivalent  to  the  OTF 
resulting  from  an  aperture  with  diameter  r0,  where  r0  is  calculated  from  the  relation 

ro  =  0.185  .  (35) 

By  incorporating  Eq.  (35)  into  Eq.  (33),  the  ‘seeing’  condition  is  directly  incorporated  into  the 
OTF: 
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#p(/M)  = 


m  (  W  r  n  k2 

JJexpj  -  32x30.033(0.185)5/,3-r-4  /$0(k)cos2(^-) 
i=l  (  To  J  I* 

X  [1  —  COs(k  •  p)}[l  —  COs(tfiK  •  $) )]dic  j 

n[e*p{  - L9€?  /  $<>(«)  cos2(^-) 

w. 

X  [1  —  COs(k  •  p  )][1  —  COs(j)iK  ■  t?)]d/c  jj  r° 


n H^p^yT, 

1=1 


(36) 


where. 


Hpi(p  ,t?)  =  exp 


(- 


967  J  $0(/c)  cos2(^^“^^  —  cos(*  ‘  P  )][1  ~  cos(T;i/c  • 


*}■ 

(37) 


and,  die  parameter  roi  is  defined  as  Fried’s  coherence  parameter  for  layer  i,  where  r~5/3  = 
W,r~5/3.  In  Eq.  (37),  the  term  Hpx  can  be  thought  of  as  the  phase  OTF  for  layer  i  calculated 
for  an  roi  of  1  meter.  Therefore,  by  calculating  Hpi  a  single  time,  we  have  a  general  phase 
OTF  that  can  now  be  scaled  to  the  actual  turbulence  profile  through  Eq.  (36).  The  strength  of 
the  ith  turbulence  layer  can  be  adjusted  by  varying  r^.  The  overall  value  of  r0  for  all  of  the 
layers  is  related  to  the  individual  ’s  by  (16) 


.=i 


(38) 


The  complete  phase  OTF  is  then  found  by  scaling  and  multiplying  die  OTFs,  Hpx,  for  each 
layer.  Results  obtained  with  this  method  are  presented  in  Section  2.8. 

No  analogy  to  the  simplifications  obtained  for  the  phase  OTF  exist  for  the  amplitude 
OTF.  This  is  due  to  the  form  of  the  amplitude  OTF  in  Eq.  (34).  However,  a  similar  analysis 
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yields 


H*(p)  = 


1+E£i 

i  +  Eli  -feiUo)  ’ 

ro* 


(39) 


with,  f  ai(p  )  representing  a  normalized  amplitude  correlation  function  for  layer  i  and  given 
by 


tai(p)  =  16jr30.033(0.185)5/l3y  $0(/c)sin2(»?,a)cos(K  •  p)dn 

—  0.983 j  $<,(«)  sin2(i?,a)  cos(ac  •  p)dZ.  (40) 

Amplitude  OTF  results  are  also  presented  in  Section  2.8. 

2.8  Numerical  results 

The  purpose  of  this  section  is  to  provide  representative  results  and  to  demonstrate  the 
utility  of  the  OTF  equations  derived  in  Section  2.7.  In  Sections  2.3  and  2.4,  the  amplitude  and 
phase  OTFs  were  derived  for  a  continuous  atmospheric  C2  profile.  These  OTF  equations  were 
derived  using  diffraction  to  account  for  both  phase  and  amplitude  variations  in  the  corrected 
wavefront  The  OTF  equations  accurately  model  the  performance  of  the  adaptive-optics 
system.  However,  the  three  dimensional  integration  is  very  time  consuming.  In  Section  2.6,  it 
was  shown  that  by  modeling  the  continuous  atmosphere  with  discrete  layers,  one  dimension 
of  die  OTF  integration  could  be  removed.  Finally,  Section  2.7,  demonstrated  how  to  create 
a  set  of  normalized  phase  OTFs,  Hpi,  and  normalized  amplitude  correlation  functions,  foi, 
that  can  be  used  to  compute  the  overall  OTF  of  file  system  through  the  use  of  the  scaling 
parameters  r^.  The  benefit  of  using  HPi  and  fa,  is  that  once  these  functions  are  calculated, 
the  OTFs  resulting  from  a  wide  range  of  atmospheric  conditions  can  be  calculated  through  the 
very  simple  scaling  operations  of  Eq.  (36)  and  Eq.  (39). 

This  section  demonstrates  the  ease  of  computing  the  OTFs  by  first  generating  a  family 
of  file  curves  for  and  fa<.  A  family  of  Hpi  curves  is  shown  in  Figure  6.  These  curves 
are  created  for  an  optical  wavelength  of  0.5  microns  and  an  object  angle,  t?,  of  2.4  /rrad 
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for  a  4  layered  atmosphere  with  the  inner  scale  set  at  Lm  =  1  millimeter  and  die  outer 
scale  set  at  La  =  5  meters.  The  object  angle  of  2.4  ^rad  is  chosen  to  be  equal  to  the 
isoplanatic  angle  as  defined  by  Fried  (10)  for  a  Hufnagel- Valley  atmosphere  with  a  54  mph 
upper  atmosphere  wind  (21).  This  angle  is  chosen  to  allow  the  direct  comparison  with  a 
geometric  optics  calculation  discussed  in  Section  2.9.  The  curves  in  Figure  6  are  computed 
such  that  when  used  in  conjunction  with  Eq.  (36)  to  obtain  the  complete  phase  OTF,  rot  in 
Eq.  (36)  is  evaluated  in  units  of  centimeters.  Once  the  tfp,  curves  are  created,  the  phase  OTF 
for  a  wide  range  of  atmospheric  conditions  can  be  computed  by  raising  Hpt  to  the  r  ~5^3  power, 
and  then  multiplying  the  four  OTFs  together.  This  becomes  a  fast  and  convenient  method  for 
determining  how  a  change  in  atmospheric  conditions  effects  die  phase  OTF. 


(A) 


fXm) 


(B)  p(m) 


Figure  6.  Hv%{p)  for  turbulent  layer  located  at  A)  200  m,  B)  2  km,  C)  10  km,  and  D)  18  km. 
These  OTFs  are  calculated  for  A  =  0.5/im,  t?  =  2.4/xrad,  and  roi  -  1  cm. 
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For  example,  we  want  to  calculate  the  phase  OTF  for  a  Hufnagel- Valley  atmosphere 
with  a  54  mph  upper  atmosphere  wind  and  a  i)  of  2.4  ^rad.  We  simply  use  the  Hpi  curves 
plotted  in  Figure  6  and  the  values  of  r0,  given  for  HV-54  in  Table  1  and  evaluate  the  equation 


HM 


n 

i=i 

(41) 


This  calculation  produces  a  curve  that  matches  the  OTF  curve  found  in  Figure  3  which  required 
many  computational  hours  to  calculate  using  a  continuous  atmospheric  C\  turbulence  profile. 

In  a  similar  manner,  a  family  of  curves  can  be  created  for  calculating  the  amplitude 
OTF  utilizing  Eq.  (39).  Unfortunately,  file  amplitude  OTF  can’t  be  scaled  with  the  same  ease 
as  the  phase  OTF.  However,  by  creating  a  family  of  curves  for  f  „,(/>),  a  simple  method  of 
calculating  the  overall  amplitude  OTF  is  obtained.  Also,  file  process  is  somewhat  simplified 
since  the  amplitude  OTF  is  independent  of  i).  By  evaluating  Eq.  (39)  at  a  wavelength  of  0.5 
microns,  the  family  of  curves  given  in  Figure  7  is  obtained.  Note  that  f  „,(/>)  has  a  positive 
value  at  p  =  0,  and  then  quickly  decays  to  zero.  Knowledge  of  the  general  shape  of  fot(p) 
and  the  amplitude  OTF  can  be  used  to  considerably  simplify  the  numerical  integration.  It  is 
known  that  Ha( 0)  =  1.0  and  then  falls  off  to  a  relatively  constant  value  as  shown  in  Figure  2 
on  page  15.  Greater  amounts  of  turbulence  will  not  significantly  change  file  cutoff  point,  but 
rather  will  change  the  constant  value  to  which  Ha(p)  decays.  A  comparison  of  Figures  7 
and  2  shows  that  when  the  constant  value  of  Ha  is  attained,  fa,(p)  «  0.  Therefore,  a  good 
approximation  of  the  constant  value  is  found  using 


The  constant  value  attained  by  the  amplitude  OTF,  designated  Ha(c),  for  a  Hufhagel- 
Valiey  atmosphere  with  a  54  mph  upper  atmosphere  wind  is  determined  using  the  plotted 
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Figure  7.  The  normalized  correlation  function,  rat(p),  for  turbulent  layer  located  at  A) 
200  m,  B)  2  km,  Q  10  km,  and  D)  18  km.  These  functions  are  calculated  for 
A  =  0.5am. 


values  for  f  ai(0)  and  the  values  of  rot  given  for  HV-54  in  Table  1  and  evaluating  Eq.  (42): 


Ha(c) 


1 


i  +  £,m= 


i  r  . 

1  T73  1  Ol 


(0) 


=  0.8082. 


(43) 


Again,  this  value  compares  very  well  with  the  constant  value  attained  by  Ha(p)  in  Figure  2. 

By  using  a  small  amount  of  memory  to  store  die  data  points  contained  in  Hpi(p)  and 
the  values  of  fa,(0),  complete  OTFs  can  be  calculated  using  a  programmable  calculator  or 
plotted  using  a  lap-top  computer.  The  utility  of  this  method  is  in  rapid  computations  under 
changing  atmospheric  conditions. 


2.9  Comparison  between  the  diffraction  and  geometric  optics  OTF 

In  this  section,  two  useful  comparisons  are  made  between  the  diffraction  method  and  the 
geometric  optics  method  of  calculating  the  OTF  of  an  adaptive-optics  system.  First,  a  direct 
comparison  is  made  of  the  OTFs  for  a  given  atmosphere.  In  this  comparison,  the  additional 
phase  and  amplitude  information  obtained  through  the  diffraction  method  is  highlighted. 
Second,  the  high  frequency  limit  of  the  OTFs  as  a  function  of  object  angle,  t).  is  examined.  It 
is  common  to  define  the  isoplanatic  angle  as  file  angle,  d,  where  file  OTF  high  frequency  limit 
has  fallen  to  1/e  of  the  d  =0  value  (10).  In  this  comparison  it  is  shown  that  the  diffraction 
method  of  calculation  predicts  a  more  optimistic  isoplanatic  angle. 

In  his  discussion  of  atmospheric  turbulence,  Goodman  (14)  defines  the  ‘near  field’ 
region  where  geometric  optics  is  valid  as  the  region  where  7?  <  jfp .  This  is  the  equivalent  of 
stating  that  coS!(3gi)  »  1  in  the  phase  OTF  expression  of  Eq.  (27)  and  that  sin2  (#)«0 
in  the  amplitude  OTF  expression  of  Eq.  (18).  With  file  cos2  term  of  Eq.  (27)  set  equal  to  one, 
and  assuming  a  strict  Kolmogorov  power  spectrum  where  $(/e,  rj)  oc  «T  file  integrations 
over  k  can  be  solved  analytically.  As  shown  in  Appendix  E,  file  substitution  of  cos2  ( )  =  1 
in  the  phase  OTF  expression  of  Eq.  (27)  yields  the  same  OTF  solution  given  by  Fried  using 
the  geometric  optics  analysis  (10).  Also  note  that  setting  sin2  <¥>-  0  in  Eq.  (18)  results  in 
Ha(p)  =  1  which  is  consistent  with  a  geometric  optics  calculation. 
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One  of  the  advantages  to  using  a  diffraction  method  to  calculate  die  system  OTF  is 
that  it  allows  for  the  separation  of  phase  and  amplitude  effects.  This  is  important  because  die 
adaptive-optics  system  is  only  correcting  for  phase  effects.  Figure  8  shows  a  direct  comparison 
between  the  diffraction  method  of  computing  die  OTF  and  the  geometric  optics  method.  In  this 
figure,  the  4  layered  turbulence  model,  as  discussed  in  die  previous  section,  has  been  used  with 
an  optical  wavelength  of  A  =  0.5  microns,  a  Hufnagel- Valley  turbulence  profile  with  a  54  mph 
upper  altitude  wind  (21),  an  object  angle  of  t?  =  2.4/irad,  and  a  Von  Karman  refractive  index 
power  spectrum  with  Lm  =  1mm  and  L0  =  5m.  The  geometric  optics  OTF  was  calculated 
using  the  same  turbulence  model  and  the  same  Von  Karmon  power  spectrum  but  with  the 
approximation  of  cos2(!Z^-)  =  1.  The  diffraction  method  gives  a  better  understanding  of  the 
atmospheric  correction  problem  and  helps  quantify  die  level  of  improvement  possible  in  die 
OTF.  For  example,  from  Figure  8  it  is  seen  that,  for  this  atmosphere,  no  matter  how  good  die 
phase  correction  is,  there  will  always  be  about  a  20%  attenuation  due  to  the  amplitude  effects. 


Figure  8.  Comparison  between  diffraction  theory  and  geometric  optics  theory.  The  diffrac¬ 
tion  OTFs  are  calculated  with  a  4  layered  atmosphere  at  200m,  2km,  10km,  and 
18km,  A  =  0.5/im,  a  Hufnagel- Valley  turbulence  profile  with  a  54  mph  upper 
altitude  wind  (21),  d  =  2.4/zrad,  and  a  Von  Karman  refractive  index  power  spec¬ 
trum  with  Lm  =  1mm  and  L0  =  5m.  The  geometric  optics  OTF  is  calculated  with 
the  same  atmospheric  conditions  but  using  the  geometric  optics  approximation. 
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In  adaptive-optics,  the  size  of  the  isopianadc  angle  ex'  isopianatic  patch  of  the  atmo¬ 
sphere  is  of  interest  This  is  commonly  thought  of  as  an  area  of  die  atmosphere  that  has 
high  correlation  in  the  index  of  refraction.  Therefore,  a  ray  of  light  propagating  at  any  angle 
through  this  portion  of  the  atmosphere  will  receive  approximately  the  same  wavefront  distor¬ 
tions.  However,  no  two  paths  through  the  atmosphere  will  induce  identical  distortions.  The 
isopianatic  angle  is  therefor  the  angle  over  which  die  difference  in  the  distortions  are  ‘not 
too  large’.  Fried  established  a  method  of  quantifying  this  measure  by  defining  the  isopianatic 
angle  as  the  angle,  t),  where  the  OTF  high  frequency  limit  has  fallen  to  e-1  of  the  t?  =  0 
value  (10).  This  definition  can  be  interpreted  two  different  ways.  First,  since  for  the  idealized 
system  modeled  here,  geometric  optics  always  calculates  an  OTF  value  of  1.0  at  d  =  0,  the 
isopianatic  angle  can  be  considered  to  be  the  angle  at  which  the  high  frequency  limit  of  the 
overall  OTF  has  fallen  to  a  value  of  e-1.  Second,  recall  that  the  adaptive-optics  system  is 
concerned  only  with  the  phase  correctability.  In  this  case,  only  the  phase  portion  of  the  OTF 
is  considered,  and  the  isopianatic  angle  is  defined  as  the  angle  at  which  the  high  frequency 
limit  of  Hp  has  fallen  to  a  value  of  e~l.  Figure  9  shows  a  comparison  between  the  OTF 
high  frequency  limit  using  diffraction  theory  and  the  OTF  high  frequency  limit  using  geo¬ 
metric  optic  theory  for  the  same  atmospheric  conditions  discussed  above.  We  see  that  where 
geometric  optics  predicts  an  isopianatic  angle  of  2.4/xrad  for  this  atmosphere,  the  diffraction 
theory  indicates  a  greater  value  regardless  of  how  the  isopianatic  angle  is  defined.  When  the 
isopianatic  angle  definition  is  taken  to  be  the  angle  at  which  the  high  frequency  limit  of  the 
overall  OTF  has  fallen  to  a  value  of  e-1,  the  isopianatic  angle  equals  3.1/rrad.  When  the 
definition  is  taken  to  be  the  angle  at  which  the  high  frequency  limit  of  Hp  has  fallen  to  a  value 
of  e_1 ,  the  isopianatic  angle  equals  3.65/xrad.  It  is  easy  to  show  that  any  atmosphere  will  yield 
a  greater  isopianatic  angle  when  a  diffraction  calculation  is  performed  than  when  a  geometric 
optics  analysis  is  used. 
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Figure  9.  Comparison  of  the  OTF  high  frequency  limit  (p  —  oo)  between  diffraction  theory 
and  geometric  optic  theory  as  a  function  of  separation  angle,  t).  The  diffraction 
OTFs  are  calculated  with  a  4  layered  atmosphere  at  200m,  21cm,  10km,  and  18km, 
A  =  0.5pm,  a  Hufhagel- Valley  turbulence  profile  with  a  54  mph  upper  altitude 
wind  (21),  and  a  Von  Karman  refractive  index  power  spectrum  with  Lm  =  1mm 
and  L0  =  5m.  The  geometric  optics  OTF  is  calculated  with  the  same  atmospheric 
conditions  but  using  the  geometric  optics  approximation,  (a)  represents  the  iso- 
planatic  angle  calculated  for  geometric  optics,  (b)  and  (c)  represent  die  isoplanatic 
angle  computed  with  two  different  definitions  for  die  diffraction  method  . 
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2.10  Summary  and  Conclusions 

In  this  chapter  a  model  is  developed  that  allows  for  die  derivation  and  evaluation  of 
an  angle  dependent  average  OTF  for  an  ideal  adaptive-optics  system.  An  expression  is 
derived  for  a  wavefront  propagating  through  a  turbulent  atmosphere  at  an  arbitrary  angle, 
d.  The  derived  expression  includes  diffraction  effects  to  properly  account  for  phase  and 
amplitude  perturbations  in  die  propagating  wave.  Expressions  are  then  derived  for  the  phase 
and  amplitude  correlation  functions,  ra(p*)  and  Tp(p,0),  and  the  average  amplitude  and 
phase  OTFs,  ( Ha(p ))  and  { Hp(p ,  0)).  These  expressions  are  used  in  the  evaluation  of  the 
SNR  in  Chapter  HI,  die  analysis  of  the  PSF  in  Chapter  IV,  and  the  development  of  a  new 
wavefront  correction  algorithm  in  Chapter  V.  All  of  the  above  expressions  are  generalized 
with  respect  to  the  refractive  index  power  spectrum,  $(77,  /c).  This  generalization  allows  an 
analysis  of  performance  effects  as  a  function  of  the  turbulent  scale  sizes.  A  useful  method  of 
computing  ( Ha(p ))  and  {Hp(p,  0))  is  also  presented  in  this  chapter.  A  layered  atmosphere  is 
created  that  results  in  an  expression  that  can  be  scaled  to  arbitrary  turbulent  conditions.  An 
actual  OTF  is  then  determined  through  a  scaling  term  of  r~^3,  where  roi  represents  Fried’s 
coherence  parameter  related  to  the  ith  layer  (9, 16).  This  scaling  method  of  evaluation  allows 
rapid  evaluation  of  the  effect  of  changing  atmospheric  conditions.  The  advantage  of  using  a 
diffraction  theory  is  demonstrated  in  a  comparison  with  geometric  optics.  It  is  shown  that  the 
diffraction  method  separates  phase  and  amplitude  effects  and  allows  for  determining  the  effect 
of  various  CI(tj)  profiles.  When  the  profile  contains  far-field  turbulence,  amplitude 
effects  become  significant  and  this  method  becomes  essential  for  accurate  calculations. 
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111.  The  angle  dependent  OTF  signal-to-noise  ratio  (SNR) 

3.1  Introduction 

In  this  chapter,  the  work  of  Chapter  n  is  extended  by  calculating  the  object  angle 
dependent  variance  of  the  OTF.  The  variance  is  then  used  in  conjunction  with  die  average 
OTF  to  obtain  a  signal-to-noise  ratio  (SNR)  performance  measure  (47).  The  SNR  performance 
measure  is  particularly  important  when  image  reconstruction  is  used  in  conjunction  with 
adaptive-optics  imaging  (43).  In  image  reconstruction,  knowledge  of  the  system  OTF  can 
be  used  along  with  the  measured  image  intensity  to  obtain  an  estimate  of  the  actual  object 
intensity. 

To  fully  understand  die  motivation  for  computing  the  SNR,  consider  the  linear,  locally 
space-invariant,  system  model  for  incoherent  adaptive-optics  imaging.  Let  t(x,t?)  represent 
the  measured  image  intensity,  where  x  designates  a  point  in  the  image  and  d  is  the  angle 
between  the  reference  source  and  the  object  The  image  i(x,  0)  can  be  related  to  die  object 
being  imaged,  o(x),  by 

i(x,  i?)  =  o(x)  *  h(x,  •d ),  (44) 

where  h(x,  d)  is  the  instantaneous,  object  angle  dependent  point  spread  function  (PSF)  of 
the  system  and  die  notation  *  designates  2 -dimensional  convolution.  Note  that  the  object  is 
considered  to  be  the  same  object  at  any  angle  and  is  therefore  not  a  function  of  tf.  The  problem 
of  imaging  an  extended  field-of-view  with  regard  to  the  angle  dependent  PSF  is  examined  in 
Chapter  IV. 

The  relationship  between  the  image  and  the  object  can  be  equivalently  stated  in  the 
frequency  domain  as 

I(v,4)  =  0(v)H(v,ti),  (45) 

where  V  is  the  spatial  frequency  variable  and  the  capital  letters  7,  O,  and  //  are  the  Fourier 
transforms  of  i,  o,  and  h,  respectively.  H(v,  t?)  is  the  instantaneous  OTF  of  the  imaging 
system.  In  a  deterministic  system  (Le.,  a  system  in  which  the  OTF  is  perfectly  known). 
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the  object  spectrum  0(i 7),  can  be  recovered  out  to  the  cutoff  frequency  of  the  system  by 
multiplying  the  measured  spectrum  I(v,  t?)  by  the  inverse  of  H(v ,  d): 

0{y)  =  l{vJ)H~\W).  (46) 


However,  imaging  through  the  atmosphere  produces  an  instantaneous  OTF  that  is  random. 
This  randomness  will  be  manifested  as  noise  in  the  image  measurement  and  will  always  tend  to 
degrade  the  quality  of  die  reconstructed  image.  The  SNR.  which  is  a  measure  of  the  variability 
or  randomness  of  die  OTF  relative  to  it’s  mean  value,  provides  a  measure  of  the  quality  with 
which  die  spectral  components  of  H(i 7,  r?)  can  be  used  to  reconstruct  the  object  spectrum  in 
Eq.  (46). 

In  Chapter  n.  it  was  shown  how  to  compute  die  average  OTF.  However,  the  average 
OTF  gives  no  indication  as  to  how  a  single  realization  of  die  OTF  can  vary  from  the  mean 
value.  The  SNR  is  used  to  quantify  the  fluctuations  in  the  OTF  and  is  defined  as 


SNR  {H(P,d)} 


\<W*))\ 

V'vfcWMj}’ 


(47) 


where  Var{H(P,  t?)}  is  the  variance  of  H(v,  d). 

This  chapter  presents  the  first  analysis  to  yield  OTF  SNR  results  that  are  a  function 
of  the  angle  between  die  reference  beacon  and  the  object  The  analysis  considers  an  ideal 
adaptive-optics  system  which  is  able  to  perfectly  measure  a  wavefront  phase  from  a  reference 
beacon,  and  in  turn  perfectly  apply  this  correction  to  die  object  wavefront  The  diffraction 
based  analysis  of  the  average  OTF  and  the  corresponding  correlation  functions  presented  in 
Chapter  n  are  fully  utilized  in  this  analysis. 

The  remainder  of  this  chapter  is  organized  as  follows.  Section  3.2  presents  an  analysis 
of  die  object  angle  dependent  variance  and  SNR  of  die  OTF.  In  Section  3.3,  simplifications 
of  die  numerical  computation  process  are  presented.  These  simplifications  allow  solutions  to 
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otherwise  computationally  prohibitive  problems.  Section  3.4  contains  numerical  results  for  a 
typical  atmosphere.  Finally,  conclusions  are  presented  in  Section  3.S. 

3.2  Derivation  of  variance  and  SNR  expressions 

In  this  section,  an  expression  is  derived  for  the  variance  of  the  OTF  of  an  adaptive-optics 
system.  The  variance  is  designated  Var{tf(p,  t?)},  where  H(p,  t?)  represents  the  instantaneous 
system  OTF.  Once  an  expression  for  the  variance  is  obtained,  the  SNR  is  calculated  using 
Eq.  (47). 

Rather  than  computing  Va r{H(p,  ■d) }  directly,  the  second  moment  of  the  OTF  is  estab¬ 
lished,  and  the  variance  is  computed  using 

Var{//(£t?)}  =  (\H(p,0)\>)  -  (H(p,0))> 

=  (\H(p,d)f)  -  Hl{p)(Ha{p  ))*(Hp(p,0))\  (48) 

where,  Ha(p)  is  die  OTF of  the  optics,  and  ( Ha(p ))  and  (Hp(p,  i?)}  are  the  average  amplitude 
and  phase  OTFs,  respectively.  The  average  phase  and  amplitude  OTPs  have  been  computed 
in  Chapter  EL  However,  the  second  moment  of  die  OTF,  (\H(p,d)\2),  still  needs  to  be 
calculated.  A  combined  phase  and  amplitude  transmission  screen,  t(x,  t?).  is  used  to  represent 
the  residual  phase  and  amplitude  perturbations  of  the  corrected  wavefront  in  die  pupil  of  the 
adaptive-optics  system.  The  transmission  screen  is  given  by 

=  ta{x)tp(x,d)  =  (1  +  a(x))exp{j0(£,t?)},  (49) 

where  a(x)  represents  the  amplitude  perturbations  and  <f>(x,  t?)  represents  the  angle  dependent 
residual  phase  in  the  pupil  of  the  adaptive-optic  system.  The  second  moment  of  the  OTF  can 


37 


now  be  written  as  (14) 


(\mm 


d iP(g)P(2~? 


) 


J  d2P(Z)*\t(2j)\3  J 


(50) 


)■ 


where  P(x)  is  the  real  valued  pupil  function  with  a  value  of  1  when  x  is  inside  the  aperture 
and  0  for  x  outside  die  aperture.  The  variable  a(x)  in  Eq.  (49),  is  the  zero  mean,  Gaussian, 
random  variable  representation  of  the  amplitude  perturbations  of  die  field  in  the  pupil.  The 
variable  <f>(x,  t?)  in  Eq.  (49),  is  the  zero  mean,  Gaussian,  random  variable  representation  of 
the  phase  difference  between  the  wavefronts  of  the  reference  beacon  and  die  object  Note  that 
the  phase  difference  is  dependent  on  the  object  angle,  d. 

Next  the  expected  value  of  the  fraction  in  Eq.  (SO),  is  assumed  to  be  approximately  equal 
to  the  expected  value  of  the  numerator  divided  by  the  expected  value  of  die  denominator  (14). 
The  integrations  over  the  pupil  function  in  Eq.  (50)  are  represented  in  simplified  form  by 


JJ  -) dfd^  =  //  *x&x'P{x)P{x  -  p  )P{x)P(x,  -p).  (51) 


Substituting  Eq.  (49)  into  Eq.  (50)  and  using  the  assumption  of  independence  between  a(x) 
and  <t>(x,  t?)  gives 


(\H(pJ)\>)  = 

f  j  dfdf'((l  +  a(x))(l  +  a(x  -  p))(  1  +  a(x'))(l  +  a(x'  -  p})) 

J  JP(i-p) 


dx&r'((l  +  a(x))2(l  +  o(x'))2) 


(52) 


We  now  compute  the  three  expected  values  in  Eq.  (52).  Note  that  the  third  order  moment 
of  jointly  Gaussian,  zero  mean  random  variables,  (/„  is  equal  to  zero  and  that  the  fourth  order 
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moment  is  given  by  (40) 


(UXU3U3UA)  =  (UtU2)(U3U4)  4  (UXU3)(U2U4)  +  {UXUA)(U2U3).  (53) 

It  can  then  be  shown  using  Eq.  (53)  that  the  fourth  order  amplitude  expectation  in  the 
numerator  of  Eq.  (52)  is  given  by 

((1  +  «(f»(l  +  a(x  -  p))(l  +  a(x'))(l  +  «(*  -  p)))  = 

(1  +  ra(p))2  +  (1  +  ra(Af))2  +  (1  +  ra(Af +p))(i  +  ra(A  s-  p))  -  2,  (54) 

where  Ax  =  x  —  x1  and  the  correlation  function,  ra(p),  is  defined  as 

ra(x-x/)  =  (a(x)a(x')).  (55) 

In  a  similar  manner,  the  fourth  order  amplitude  expectation  in  die  denominator  of  Eq.  (52)  is 
computed  to  be 

((1  +  a(x))2(l  +  a(x'))2)  =  (1  +  To(0))2  4-  4ro(Ax)  +  2r2(Ax).  (56) 

Next,  the  fourth  order  moment  is  computed  for  of  the  exponential  term  in  the  numerator  of 
Eq.  (52).  Since  <f>(x,  t?)  is  a  zero  mean,  Gaussian  random  variable,  the  fourth  order  moment 
can  be  written  as 

exp  {-i([«*,0)  -  Ml- pj))  -  (MU)  -  Ml  -  *  *))]’)  } 

=  exp{  -  -  Ml  -  PJ)f  +  (MU)  ~  Ml  ~  U)f 

-  2(MU)  -  Ml  -  p,mMU)  -  Ml  -  P.  *)))]}•  (57) 
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Expanding  die  first  two  squared  term  in  the  exponent  of  Eq.  (57)  yields 


mi  *)-*(*-£  I?))2)  =  2rp(o,  t?)  -  2rp(p,  t?),  (58) 


where  the  phase  correlation,  r P(p,  t?),  is  defined  as 

rp(f-z',t?)  =  (59) 


The  last  term  in  die  exponent  of  Eq.  (57)  can  be  expanded  and  written  in  terms  of  phase 
correlations  in  a  straight  forward  manner.  Equation  (57)  is  now  only  a  function  of  second 
order  phase  correlation  functions  and  is  given  by 

=  exp |2rp(p,  i?)  -  2rp(0,  t?)  +  2rp( Ax,  t?) 

-  rp(Ax  +  p,  t?)  -  rp(Af  -  p,  7?) J,  (60) 

where.  Ax  is  the  separation  between  correlation  points  (i.e.,  Ax  —  x  —  a?).  Equa¬ 
tions  (54),  (56),  and  (60)  are  now  combined  to  yield  an  expression  for  (|//(p,  t?)|2).  For 
convenience,  the  result  is  broken  into  numerator  (Num)  and  denominator  (Den)  expressions: 


Num(|//(p,  i?)|2)  =  (1  +Ta(p  ))2exp{2rp(p,t?)  -2rp(0,t?)} 
x  //  dxdx-fl  +  (1  ±  ra(Ax))2  +  (1  +  r,( Ax +/))(!  ±  ro(Ax  -  p ))  -  2N 

JJp(s-p)  \  (i  +  ra(P  ))2  / 

x  exp{2rp( Ax,  t?)  -  rp(Ax  +  p,  t?)  -  Tp(Ax  -  p,ti)}, 

(61) 

and. 


Den(|ff(ft  t»P>  =  (1  +  r„(0))2  J fpir)  (l 


+ 


4ro(Ax)  +  2r2(A 
(i  +  ra(o))2 


2), 


(62) 


where  the  terms  not  dependent  on  x  or  x'  have  been  factored  outside  the  integration.  Notice 
that  when  Eqs.  (61)  and  (62)  are  combined  as  numerator  and  denominator,  there  are  terms 
outside  die  integrals  that  correspond  to  { Ha(p ))2  and  (Hp(p,  t?))2  as  shown  in  Eqs.  (11) 
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and  (21).  Also  note  that  the  integration  in  Eq.  (62)  can  be  evaluated  independent  of  p  and  t). 
This  independence  is  useful  for  the  numerical  analysis  in  Section  3.4. 

Equations  (61)  and  (62),  are  now  combined  and  substituted  into  Eq.  (48),  to  yield  a  final 
expression  for  the  OTF  variance  as 

Var {H(p,t)}  =  (63) 

where  T(/T,t?)  is  defined  by 


T  (pj)  = 

Jj  did?!  (l  + 

x  exp{2r,(Af,d)  -r,(Af  +pj)  -  r„(Ax  -p,d)} 


4ra(Aj?)+2r;(Ag) 

(i+ra(0))i 


It  is  now  a  simple  matter  to  substitute  Eq.  (63)  into  Eq.  (47)  to  obtain  the  desired  object  angle 
dependent  SNR: 

smwM-  ™ 


Note  that  die  SNR  is  a  function  of  the  phase  and  amplitude  correlation  functions,  Tp 
and  ra.  The  amplitude  correlation,  ro,  was  derived  in  Chapter  n  with  the  result: 

ra(p  )  =  4vk2jJ d?7dK$(«,i7)sin2(^^-)cos(«  •  p).  (66) 


A  typical  amplitude  correlation  function  is  shown  in  Figure  10.  The  correlation  function  in 
Figure  10  is  calculated  for  A  =  0.5  microns  using  a  Hufnagel- Valley  turbulence  profile  with  a 
54  mph  upper  atmospheric  wind  (21),  and  a  Von  Karman  index  of  refraction  power  spectrum 
with  inner  scale  set  at  1  millimeter  and  die  outer  scale  set  at  5  meters. 
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Figure  10.  Amplitude  correlation  function,  ra(p),  plotted  versus  p.  This  plot  was  calculated 
for  A  =  0.5  microns  using  a  Hufnagel- Valley  turbulence  profile  with  a  54  mph 
upper  atmospheric  wind.  This  turbulence  profile  has  a  Fried  coherence  cell  size 
of  ra  =  5  cm  and  an  isoplanatic  angle  of  60  —  2.4  //rad. 


The  phase  correlation,  rp,  was  also  derived  in  Chapter  II  with  the  result: 

r p(Pi'O)  =  8nk2  J J  4>(ic,  t])  cos2  j  cos(*  •  )[l  —  cosmic  •  i?)]d7?d/c.  (67) 


A  typical  phase  correlation  function  is  shown  in  Figure  11.  The  correlation  function  in 
Figure  1 1  is  calculated  for  A  =  0.5  microns  using  a  Hufnagel- Valley  turbulence  profile  with  a 
54  mph  upper  atmospheric  wind  (21).  The  object  angle  ■d  is  set  equal  to  the  isoplanatic  angle, 
0O,  as  defined  by  Fried  (10). 


33  Simplifications  for  evaluation  of  the  SNR 

In  this  section,  simplifications  are  presented  that  are  required  to  numerically  evaluate 
the  SNR  expression  given  in  Eq.  (65),  and  specifically  the  evaluation  of  T(£  t?)  in  Eq.  (64). 
In  the  numerator  of  T (p,  d),  the  integration  over  x  and  £  is  a  four  dimensional  integration 
over  an  overlapping  region  of  shifted  pupil  functions.  Since  the  evaluation  of  the  correlation 
functions,  T a(p)  and  Tp(p,  3),  involve  three  dimensional  integrations,  as  shown  in  Eqs.  (66) 
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Shift  Distance  in  the  Aperture  Plane  ~  p(m) 

Figure  11.  Phase  correlation  function,  rp(/>,  i9),  plotted  versus  p.  This  plot  was  calculated 
for  A  =  0.5  microns  and  d  =  0o  using  a  Hufnagel- Valley  turbulence  profile  with 
a  54  mph  upper  atmospheric  wind.  This  turbulence  profile  has  a  Fried  coherence 
cell  size  of  r0  =  5  cm  and  an  isoplanatic  angle  of  0„  =  2.4  /irad. 

and  (67),  the  complete  evaluation  of  T(p,  $)  requires  a  seven  dimensional  integration.  Even 
with  present  day  computational  power,  seven  dimension  numerical  integration  is  extremely 
intensive  for  all  but  very  simply  functions.  In  addition,  file  integrands  of  Eqs.  (66)  and  (67) 
are  oscillatory  and  require  relatively  fine  subdivisions  to  obtain  convergence.  It  therefore 
becomes  computationally  prohibitive  to  numerically  integrate  the  function  in  its  derived  form. 

The  simplification  process  will  be  presented  in  two  steps.  First,  a  normalized  version  of 
the  correlation  function  will  be  derived.  This  will  allow  for  scaling  under  changing  atmospheric 
conditions.  Next  an  approximating  function  will  be  determined  for  the  normalized  correlation 
functions.  This  process  will  allow  for  scaling  under  changing  atmospheric  conditions  as  well 
as  changing  geometrical  conditions  (i.e.,  changing  values  of  the  object  angle,  i?). 

3.3.1  The  normalized  correlation  function  for  layered  turbulence.  The  first  simplifi¬ 
cation  involves  using  a  layered  atmosphere  as  discussed  in  Section  2.6.  By  properly  modeling 
the  atmospheric  turbulence  with  relatively  few  discrete  layers,  the  integration  over  the  altitude 
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variable,  r?,  in  Eqs.  (66)  and  (67)  can  be  replaced  with  a  summation  over  the  number  of  layers. 
This  is  a  significant  simplification,  but  die  evaluation  of  Bq.  (64)  still  requires  six  dimensions 
of  integration. 

The  next  simplification  is  achieved  by  realizing  that  the  correlation  functions,  Ta  and 
Tp  in  Eq.  (64),  are  simply  functions  of  the  shift  vectors  p,  Ax,  Ax  +  p,  and  Ax  -  p .  Since 
the  correlation  functions  are  themselves  integrations,  they  tend  to  be  relatively  smooth  as  a 
function  of  their  argument  The  smoothness  of  the  correlation  functions  is  shown  in  the  plots 
of  ro(p)  and  rp(p,i?)  in  Figures  10  and  11.  Due  to  the  smooth  nature  of  the  correlation 
functions,  they  can  be  represented  over  their  complete  argument  range  by  a  relatively  sparse 
grid  of  precomputed  values.  The  precomputed  grid  greatly  reduces  computation  times  and 
the  smooth  nature  of  the  correlation  functions  allows  simple  bilinear  interpolation  between 
grid  points  to  yield  accurate  approximations.  The  phase  correlation  function  is  found  to  be 
symmetric  about  file  x  and  y  components  of  the  shift  variable,  p,  and  the  amplitude  correlation 
function  is  circularly  symmetric.  These  two  symmetries  reduce  file  number  of  grid  points 
needed  to  completely  characterize  ra(/j)  and  Tp(p,  i9).  The  evaluation  of  T(p)  then  involves 
four  dimensional  integration  over  an  integrand  that  is  precomputed  on  a  grid  of  points.  The 
integration  in  the  numerator  of  Eq.  (64)  must  be  recomputed  for  each  p,  but  the  denominator 
is  independent  of  p  and  is  only  computed  once. 

Although  this  method  of  precomputing  file  correlation  functions  was  necessary  to  obtain 
SNR  results,  creating  the  grid  of  points  is  still  computationally  intensive.  Also,  a  disadvantage 
of  precomputing  the  correlation  functions  is  that  the  results  are  only  valid  for  a  single  scenario, 
defined  by  the  object  angle,  t?,  and  the  particular  atmospheric  C^(tj)  profile.  The  quantity 
C%(tj)  is  the  altitude  dependent  structure  constant  of  the  refractive  index  fluctuations  and  is  a 
measure  of  the  strength  and  distribution  of  the  turbulence.  It  would  be  useful  to  have  a  method 
of  quickly  computing  the  SNR  that  is  not  limited  to  a  single  scenario  but  is  robust  in  both  the 
choice  of  object  angle,  d,  and  atmospheric  C*(ri)  profile. 

The  first  step  in  developing  a  robust  analysis  method  is  to  model  file  atmosphere  with  a 
set  of  preselected  turbulent  layers.  Section  2.6  demonstrated  that  four  discrete  layers  placed  at 
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200m,  2km,  10km,  and  18km,  is  adequate  for  modeling  a  wide  range  of  atmospheric  profiles. 
Using  a  4  layered  atmospheric  turbulence  model,  a  normalized  index  of  refraction  power 
spectrum  given  by,  $D(/c)  =  g^f^r,  and  an  analysis  method  similar  to  the  development  of 
Eqs.  (36)  and  (39)  in  Chapter  n,  the  correlation  functions  given  in  Eqs.  (66)  and  (67)  can  be 
rewritten  as 

r  =  1-967  -jfi  J  ^0(«)cos2(^^-)cos(k  •  p)[l  -  cos(r?,/c- i?)]d/c 

«=1  roi 

=  E-gJsWJ),  (68) 

<=1  roi 

and, 

ra(p)  =  0.9835  -gyj  J  $0(«)sin2(^^-)cos(«  •  p  )d/c 

'=1  roi 

=  E4j3^(^  (69) 

i=i  rot 

where  the  parameter  rot  is  defined  as  Fried’s  coherence  parameter  (9)  for  layer  i  located  at 
altitude  and  is  a  measure  of  die  turbulence  strength  at  that  altitude.  Values  of  rot  for  die 
4  layer  atmospheric  model  of  several  common  atmospheric  profiles  are  found  in  Table  1  on 
page  23. 

3.3.2  Approximating  the  normalized  correlation  function.  In  Eqs.  (68)  and  (69),  die 
terms  rpt  and  fat  represent  the  normalized  correlations  for  layer  i  which  are  then  scaled  for 
the  particular  turbulence  distribution  through  r0i .  This  scaling  allows  a  new  SNR  to  be  quickly 
recomputed  when  atmospheric  conditions  change  (i.e.,  a  changing  turbulence  distribution). 
However,  the  desired  robustness  with  respect  to  die  object  angle  has  not  yet  been  attained. 
Due  to  the  d  dependence  of  the  phase  correlation  function,  fp,(p,  tf)  must  be  recomputed  for 
every  new  object  angle,  d.  The  robustness  in  object  angle  is  accomplished  by  first  computing 
rpi(p,  d )  over  a  grid  of  points  for  several  values  of  i?.  The  grid  points  are  then  used  to  obtain 
a  single  analytic  function  that  approximates  fp,(p,  d)  over  the  range  of  t?’s.  Note  that  this 


45 


approximation  for  fp,(/T,  t?)  will  be  a  three  dimensional  function  in  t?  ,  p  ,  and  0P,  where  p 
and  0P  represent  the  magnitude  and  angle  of  the  shift,  p.  The  amplitude  correlation  function  is 
independent  of  0  and  circularly  symmetric,  consequently  it  is  a  much  easier  process  to  obtain 
a  one  dimensional  fit  in  p  that  approximates  the  amplitude  correlation  data.  For  rp,(p,  0),  the 
form  of  the  approximating  function  found  to  yield  accurate  results  is  a  rational  polynomial 
given  by 


/i,W  +  /2,(gP,d)/)  +  /3t(gp,t?)p2 

1  +  f4i(0p,ti)p  +  f5i(0p,0)p2  +  f6i(0p,  O'jp3 


(70) 


where. 


/«(*) 
MW) 
MW) 
MW) 
MW) 
/«(*,.  *) 


aU  +  i0  +  a3t  0*  +  fl4t^3  +  a5t^4? 
06t  +  <*7t^  +  (08«  +  09tl f)0p, 
oiot  +  +  (ai2«  +  Ol3  i0)0p, 

ct14*‘  + 

ai8«  +  ai9i#  -f  (a2o i  4-  a2u0)0p, 
a22l  +  a23li?  +  (a24,  +  a25i^)^p, 


(71) 


and,  a„i  are  the  coefficients  tiiat  must  be  found  to  best  approximate  r pi(p,  0P,  0).  The  values 
of  the  coefficients  are  found  using  a  non-linear  least  squares  fitting  routine  that  is  slightly 
modified  from  the  standard  algorithm  given  in  (39)  for  the  inclusion  of  multiple  dimensions. 
The  coefficient  values  for  each  of  the  four  layers  are  presented  in  Table  2  for  a  refractive  index 
power  spectrum  that  has  been  normalized  to  an  r0  of  1  cm.  It  is  determined  that  for  accurate 
SNR  results,  it  is  best  to  limit  the  range  of  0 ’s  in  the  approximating  function.  The  coefficients 
in  Table  2  are  accurate  for  a  range  of  object  angles  between  2  and  14  /xrads  and  p  out  to  2 
meters.  Accurate  SNR  results  can  be  obtained  over  a  much  larger  range  of  i?  by  creating  a 
piecewise  set  of  approximating  functions  out  to  the  desired  0. 
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Table  2.  Coefficients  for  the  phase  correlation  approximating  function,  f  6P,  d)  for  t 

layered  atmosphere.  Layer  1  is  placed  at  200  meters,  layer  2  at  2  km,  layer  3 
10  km,  and  layer  4  at  18  km.  The  units  of  p,  8P,  and  tf,  are  meters,  radians,  and 
micro-radians,  respectively.  Using  these  units,  r*,  in  Eq.  (68)  is  given  in  units  of 
centimeters. 


Coefficient 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

<*i 

0.67219E-02 

0.40004E-01 

0.45689E+01 

-0.24508E+01 

a  2 

-0.49126E-02 

-0.71964E-01 

-0.32579E+01 

0.23486B+01 

<*3 

0.42528E-02 

0.20869E+00 

0.34852E+01 

0.72594B+01 

a4 

-0.69582E-04 

-0.45049E-02 

-0.14006B+00 

-0.36121E400 

a5 

-0.14925E-05 

-0.307 50E-04 

0.74663E-03 

0.45284E-02 

a6 

0.85394E-03 

0.10410E+02 

0.20650E+03 

a7 

0.24471E-01 

-0.57929E+01 

-0.33821E*01 

-0.15031E+03 

a s 

-0.43777E+00 

-0.12087E402 

-0.19866E+03 

-0.23592E+03 

Og 

0.30305E+00 

0.79010B+01 

0.15233E+03 

0.16623E403 

aio 

0.57525E+00 

-0.21106E+01 

0.32366B+03 

-0.16138E+02 

an 

-0.43740E+00 

-0.11245E+00 

-0.20472E+03 

0.16954E+02 

a12 

-0.27234E+00 

0.27593E+01 

-0.18859B403 

0.54170E+02 

a13 

0.20937E+00 

-0.13095B+01 

-0.39813B+02 

Oi4 

0.18175B+03 

0.14392E+02 

0.63275E+01 

«15 

-0.41943E+01 

-0.21243E+01 

-0.10887E+01 

-0.68352E+00 

a16 

-0.30247E+02 

-0.12206E+02 

0.77047E+01 

0.15286E+01 

Ol7 

-0.36634E+00 

0.43943E+00 

0.141 12E+00 

«18 

0.63805E+02 

-0.81229E+01 

0.13558B+03 

0.352 14E+02 

«19 

-0.87324E+01 

0.31468E+01 

-0.4915jE401 

-0.93767E+00 

a20 

0.35184B+03 

0.483 18E+02 

-0.11003E402 

-0.76625E+01 

a21 

-0.13716E+02 

-0.32610E+01 

-0.26197E400 

0.12417E+00 

a22 

0.14140E+04 

0.14184E+03 

0.15687E+03 

0.25237E+01 

a23 

-0.66265E+02 

-0.10011E+02 

-0.38129E+00 

a24 

-0.48690E+03 

-0.35202E402 

-0.14124E+02 

0.40961 E+01 

a25 

0.23843E+02 

0.33274B+01 

0.10753E+01 

-0.11481E+00 
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Table  3.  Coefficients  for  die  amplitude  correlation  approximating  function,  fa,(p)  for  a  4 
layered  atmosphere.  Layer  1  is  placed  at  200  meters,  layer  2  at  2  km,  layer  3  at 
10  km,  and  layer  4  at  18  km.  With  the  shift  parameter,  p,  in  units  of  meters,  rot, 
is  given  in  units  of  centimeters.  Coefficients  are  valid  for  all  shift  values  and  all 
object  angles. 


Coefficient 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

h 

0.15000E-01 

0.5000E-01 

0.1000E+00 

0.1800E400 

h 

0.28791E+00 

0.19884E401 

0.75052E+01 

0.12246E+02 

h 

-0.9734SE+02 

-0.301 12E+02 

-0.15856E+02 

-0.72605E401 

h 

-0.53476E+02 

-0.11641E+03 

-0.19928E+03 

-0.23363E403 

h 

0.10012E+05 

0.11210E404 

0. 32861 E403 

0.13730E+00 

h 

0.22359E+04 

0.15350E+04 

0.12089E+04 

0.97764E+03 

67 

0.12556E+03 

-0.37670E+04 

-0.16727E+04 

0.10287E+04 

For  die  amplitude  correlation  function,  rat  ,  die  same  rational  polynomial  form  worked 
well  for  die  approximating  function.  The  form  used  is  given  by 


f «(/?)  «  (u(p)  “  U(p  -  bu ))  - 


h i  +  hiP  +  hip2 


+  hip  +  hip2  +  hiP3  ’ 


(72) 


where  u(p)  is  the  unit  step  function  widi  a  value  of  one  when  die  argument  is  greater  than 
zero.  The  purpose  of  die  combination  of  step  functions  is  to  limit  the  required  range  of  the 
fitting  function.  It  is  determined  that  the  value  of  the  amplitude  correlation  function  beyond 
p  =  hi  is  approximately  zero.  The  values  of  the  seven  coefficients  for  each  of  the  four  layers 
are  given  in  Table  3  for  a  refractive  index  power  spectrum  normalized  to  an  roi  of  1  cm.  The 
coefficients  in  Table  3  are  accurate  for  all  values  of  the  shift  parameter,  p. 

The  end  result  is  approximating  functions  given  by  fpt  (p,  0P,  d)  and  f  „,(/>)  in  Eqs.  (70) 
and  (72)  that  completely  represent  the  phase  and  amplitude  correlation  functions  of  each  of 
the  four  layers.  These  functions  are  valid  over  a  wide  range  of  object  angle,  t?,  independent  of 
the  atmospheric  C2(rj)  pi  -die  and  yield  extremely  accurate  SNR  results.  The  approximating 
functions  developed  here  for  the  purpose  of  evaluating  the  SNR  of  an  adaptive-optics  system 
can  also  be  used  to  yield  accurate  calculations  of  (i/p(p,  t?))  and  ( Ha(p )).  Chapter  n 
demonstrated  how  to  combine  and  scale  OTFs  for  individual  layers  for  rapid  evaluation  of 
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the  overall  average  OTF  under  changing  atmospheric  conditions.  The  use  of  rpt(p,  0P,  t?) 
and  fol(p),  along  with  r allow  { Hv(p ,  d))  and  (Ha(p))  to  be  calculated  through  Eqs.  (11) 
and  (21).  These  OTPs  can  be  rapidly  evaluated  under  changing  atmospheric  conditions  as 
well  as  die  choice  of  object  angle,  t). 

3.4  Numerical  results 

In  this  section  numerical  results  are  presented  for  the  OTF  variance  and  SNR  as  defined 
in  Section  3.2  by  Eqs.  (63)  and  (65)  respectively.  The  theory  developed  in  Sections  3.2  and  3.3 
allow  for  calculations  of  the  variance  and  SNR  over  a  wide  range  of  atmospheric  conditions. 
However,  in  this  section,  results  are  calculated  using  a  Hufnagel- Valley  turbulence  profile 
with  a  54  mph  upper  atmospheric  wind,  rQ  =  5  cm,  and  0O  =  2.4  prad  (21).  The  expressions 
in  Eqs.  (63)  and  (65)  allow  calculations  for  any  orientation  of  p,  however,  to  the  results 
presented  here,  the  orientation  of  p  is  fixed  to  be  along  the  x-axis  of  the  pupil  (Le.  9P  —  0). 
Recall  that  d  is  also  assumed  to  be  directed  along  the  x-axis.  Subsequent  results  will  also 
consider  p  directed  along  the  y-axis  of  die  pupil  (Le.  9P  =  f ). 

Figures  12  and  13  plot  the  variance  of  the  system  OTF  as  a  function  of  where  D  is 
the  diameter  of  the  imaging  system  aperture.  The  variance  plotted  in  Figure  12  is  for  D  =  1 
meter  and  in  Figure  13  to  D  =  2  meters.  Each  figure  plots  the  variance  to  object  angles, 
i),  of  1,  2.  and  3  times  the  isoplanatic  angle,  0O.  Note  the  general  decrease  in  variance  with 
increasing  i).  By  examining  Eqs.  (63)  and  (64),  it  is  clear  that  die  decrease  is  mainly  due 
to  the  reduction  in  die  phase  OTF,  ( Hp(p ,  t))).  Also  notice  that  the  variance  decreases  with 
increasing  aperture  diameter.  As  the  aperture  diameter  increases,  the  OTF  is  formed  by  the 
correlation  of  points  that  are  separated  by  proportionately  larger  distances.  The  value  of  the 
OTF  can  be  thought  of  as  calculated  through  a  type  of  averaging  over  all  correlation  values. 
As  the  separation  between  correlation  points  increases,  the  correlation  value  will  decrease  and 
tend  toward  zero  with  very  little  variance  between  individual  realizations.  Therefore,  as  the 
aperture  diameter  increases,  the  result  of  die  ‘averaging’  over  all  die  correlation  contributions 
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will  be  a  decrease  in  the  overall  variance.  Again,  by  examining  Eqs.  (63)  and  (64).  it  can  be 
shown  that  the  variance  approaches  zero  as  the  aperture  diameter  approaches  infinity. 

The  corresponding  SNR  calculations  are  plotted  in  Figures  14  and  IS.  Once  again,  the 
SNRs  are  plotted  as  a  function  of  £  for  aperture  diameters  D  =  1  and  2  meters  and  for  object 
angles,  d  =  1 , 2,  and  3 90.  In  addition,  Figure  14  shows  a  comparison  between  calculating  the 
SNR  using  the  method  of  precomputing  the  correlation  function  for  each  particular  t),  and  the 
method  utilizing  die  approximating  functions  as  discussed  in  Section  3.3.  The  two  methods 
produce  nearly  identical  SNR  results. 

One  use  of  the  SNR  is  found  by  defining  the  effective  cutoff  frequency  of  an  adaptive- 
optics  imaging  system.  The  effective  cutoff  frequency,  pc,  can  be  defined  as  the  highest  spatial 
frequency  for  which  SNR(p)  is  greater  than  or  equal  to  (me: 

pc  —  p  such  that  SNR(p)  =  1.0  .  (73) 

The  effective  cutoff  frequency  can  be  thought  of  as  the  highest  spatial  frequency  where  the 
OTF  is  known  with  a  minimum  acceptable  degree  of  certainty.  To  illustrate  pc  for  die  results 
presented  in  Figures  14  and  15,  a  horizontal  line  is  drawn  at  die  level  where  SNR {H}  =  1. 
Vertical  lines  indicate  pc  for  the  particular  values  of  t?.  Notice  that  at  d  =  60,  pc  is  near  the 
diffraction  limit  for  die  particular  aperture.  This  is  consistent  with  Fried’s  definition  of  die 
isoplanatic  angle  where  the  wavefront  phase  of  the  object  is  expected  to  be  well  correlated 
with  the  wavefront  phase  of  the  beacon  and  therefore  relatively  correctable  (10). 

Note  that  the  SNR  significantly  decreases  as  t?  increases.  This  decrease  is  due  to  the 
dominant  effect  of  die  average  OTF  in  the  SNR  calculation  for  a  given  aperture  diameter. 
Recall  that  the  SNR  is  the  ratio  of  the  average  OTF  to  the  square  root  of  the  variance.  The 
average  OTF  decreases  with  increasing  d  much  faster  than  the  decrease  in  the  variance.  These 
two  figures  also  show  that  for  the  cases  of  6  equal  to  1  and  2  90,  an  aperture  diameter  of  2 
meters  yields  a  significantly  increased  SNR  over  an  aperture  diameter  of  1  meter.  For  the 
t)  =  26 o  case,  die  SNR  =  1  frequency  cutoff  increases  from  34  %  to  55  %  of  the  diffraction 
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limit  as  D  increases  from  1  to  2  meters.  At  t?  =  2 9ot  the  variance  decreases  faster  than 
(Hp(p,d))  in  going  from  D  —  1  to  2  meters.  However,  note  that  at  t?  =  3 90,  the  SNR  =  1 
cutoff  frequency  decreases  in  going  from  D  =  1  to  2  meters.  At  t?  =  3 0o,  { Hp(p ",  t?))  is 
dominating  the  SNR  calculation.  This  type  of  evaluation  is  important  in  system  design  issues. 

Figure  16  is  a  plot  of  ^  as  a  function  of  t ?/0o.  This  plot  includes  the  SNR  =  1  cutoff 
frequencies  of  Figures  14  and  15,  and  extends  the  cutoff  frequency  results  out  to  0  =  60 a. 
Note  that  beyond  0  =  30o,  die  SNR  =  1  cutoff  frequency  levels  off  at  a  value  below  10  % 
of  die  diffraction  limit 

Figure  17  shows  die  relation  between  the  SNR  as  calculated  for  values  of  p  along  the 
i -axis  and  the  SNR  with  p  along  the  y-axis.  The  SNRs  are  plotted  as  a  function  of  ^  for 
an  aperture  diameter  of  2  meters  and  for  object  angles,  t?,  of  1,  2,  and  30o.  In  each  case  the 
SNR  has  increased  in  moving  from  die  shift  axis  parallel  to  t?,  to  the  shift  axis  perpendicular 
to  0.  The  increased  SNR  along  the  perpendicular  axis  is  expected  since  the  phase  correlation 
as  well  as  the  average  phase  OTF  is  larger  along  the  perpendicular  axis. 


Figure  12.  Variance  of  H(p,d )  versus  £  for  A  =  0.5  microns  using  a  Hufnagel- Valley 
turbulence  profile  with  a  54  mph  upper  atmospheric  wind,  r0  =  5  cm,  and  D  =  1 
meter.  Individual  plots  are  for  d  —  1, 2,  and  30o. 
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Figure  13.  Variance  of  H(p,  i?)  versus  £  for  A  =  0.5  microns  using  a  Hufnagel- Valley 
turbulence  profile  with  a  54  mph  upper  atmospheric  wind,  r0  =  5  cm,  and  D  =  2 
meters.  Individual  plots  are  for  t?  =  1, 2,  and  36 0. 


Figure  14.  OTF  SNR  versus  £  for  A  =  0.5  microns  using  a  Hufnagel- Valley  turbulence 
profile  with  a  54  mph  upper  atmospheric  wind,  ra  =  5  cm,  and  D  =  1  meter. 
Individual  plots  are  for  0  =  1,  2,  and  39 0.  Note  that  double  lines  at  each  value 
of  0  indicate  relative  accuracy  of  calculation  method  as  stated  in  Section  3.4. 
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Figure  15.  OTF  SNR  versus  £  for  A  =  0.5  microns  using  a  Hufnagel- Valley  turbulence 
profile  with  a  54  mph  upper  atmospheric  wind.  rQ  =  5  cm,  and  D  =  2  meters. 
Individual  plots  are  for  t?  =  1. 2.  and  30o. 


Figure  16.  Normalized  cutoff  frequency,  vs.  Data  is  calculated  for  A  =  0.5  microns 
using  a  Hufnagel- Valley  turbulence  profile  with  a  54  mph  upper  atmospheric 
wind  and  r0  =  5  cm. 
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Figure  17.  OTF  SNR  versus  £  for  A  =  0.5  microns  using  a  Hufiiagel- Valley  turbulence 
profile  with  a  54  mph  upper  atmospheric  wind,  r0  =  5  cm,  and  D  =  2  meters. 
Individual  plots  are  for  d  =  1,  2,  and  30o.  This  figure  compares  the  SNR 
calculated  with  p  along  the  z-axis  to  the  SNR  calculated  with  p  along  the  y-axis. 

35  Summary  and  Conclusions 

An  analysis  of  the  object  angle  dependent  OTF  signal-to-noise  ratio,  SNR(p,  d),  of  an 
adaptive-optics  system  is  presented  in  this  chapter.  This  analysis  allows  for  diffraction  in 
the  propagating  wavefronts  which  accounts  for  both  residual  phase  and  amplitude  effects.  It 
is  found  that  the  SNR  can  be  written  in  terms  of  second  order  correlations  of  the  residual 
phase  and  amplitude  in  the  system  pupil  given  by  Tp(p,  tf)  and  r„(/>),  respectively.  The 
SNR  analysis  results  in  seven  dimensional  integrations.  Simplifications  are  presented  that 
utilize  a  layered  atmospheric  model  and  normalized  approximation  functions  of  rp(p,  t?)  and 
r„(p).  The  approximating  functions  are  robust  in  the  choice  of  t?,  and  can  be  scaled  through 
Fried’s  coherence  parameter  related  to  the  ith  layer,  roi.  The  layered  atmospheric  model  and 
normalized  approximating  functions  allow  for  relatively  quick  calculations  under  a  wide  range 
of  atmospheric  conditions  and  correction  geometries.  Numerical  results  are  presented  that 
demonstrate  die  utility  of  an  SNR  analysis  as  applied  to  an  adaptive-optics  imaging  system. 
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IV.  The  angle  dependent  point  spread  function  ( PSF )  and  Strehl  ratio 

4.1  Introduction 

This  chapter  presents  an  analysis  of  the  angle  dependent  point  spread  function  (PSF) 
and  Strehl  ratio.  Knowledge  of  the  PSF  is  useful  in  a  conventional  linear,  shift  invariant 
imaging  system,  because  the  image  intensity  is  determined  by  convolving  the  PSF  with  the 
object  intensity.  However,  since  an  adaptive-optics  imaging  system  yields  an  angle  dependent 
PSF,  the  system  is  not  shift  invariant  The  image  intensity  of  a  non-shift  invariant  system  is 
a  superposition  integral  of  the  shift  dependent  PSF  with  the  object  intensity.  In  the  somewhat 
simpler  case  of  astronomical  imaging  where  the  image  may  be  many  stars  in  a  wide  field-of- 
view  star  field,  each  star  may  be  considered  a  separate  point  source.  Assuming  that  each  star  is 
separated  from  surrounding  stars  such  that  the  corresponding  PSFs  do  not  overlap  (me  another, 
the  image  of  each  star  is  simply  die  point  source  convolved  with  die  PSF  corresponding  to  the 
star’s  angular  location.  The  angle  dependent  PSF  results  in  each  star  within  a  star  field  having 
a  slightly  different  spread  or  blur.  The  Strehl  ratio  is  also  a  useful  measure  in  adaptive-optics 
since  it  defines  how  much  the  peak  of  the  PSF  has  been  reduced  compared  to  the  peak  of  a 
diffraction  limited  PSF. 

The  organization  of  this  chapter  is  as  follows.  A  discussion  of  the  analysis  procedure 
for  the  PSF  and  Strehl  ratio  is  presented  in  Section  4.2  followed  by  the  results  of  the  analysis 
for  a  particular  atmospheric  turbulence  profile  in  Section  4.3.  Finally,  Section  4.4  summarizes 
this  chapter  and  discusses  the  significance  of  the  results. 

4.2  Analysis 

In  general,  die  PSF  is  defined  as  the  inverse  Fourier  Transform  of  the  OTF.  In  Chapter  n, 
the  angle  dependent  OTF  was  derived.  Therefore,  the  angle  dependent  PSF  can  be  defined 
here  as  die  inverse  Fourier  Transform  of  the  angle  dependent  OTF.  No  attempt  is  made  in  this 
analysis  to  define  an  analytical  expression  for  the  angle  dependent  PSF.  Rather,  for  a  particular 
atmospheric  turbulence  profile,  the  OTF  expression  of  Chapter  n  is  used  to  create  discrete 
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values  of  the  angle  dependent  OTF  which  are  then  inverse  Fourier  transformed  using  standard 
Fast  Fourier  Transform  (FFT)  techniques.  The  resulting  PSF  is  then  examined  with  regard 
to  the  angular  dependence.  Since  the  OTF  is  not  circularly  symmetric,  the  PSF  analyzed 
here  will  not  be  circularly  symmetric.  To  explore  the  non-symmetric  property  of  the  PSF. 

1-D  slices  of  the  2-D  PSF  are  plotted  along  the  axis  parallel  to  die  object  angle  separation 

— ♦ 

direction,  i),  as  well  as  along  the  perpendicular  axis.  The  angular  dependence  of  die  spread  of 
the  PSF  is  shown  by  plotting  the  Full-Width-Half-Max  (FWHM)  as  a  function  of  the  object 
angle,  i?.  As  the  name  implies  FWHM  is  a  measure  of  how  wide  a  function  is  when  the  value 
of  the  function  has  reached  half  of  the  maximum  value. 

The  Strehl  ratio  is  defined  as  the  ratio  of  die  peak  of  the  PSF  compared  to  the  peak  of 
a  diffraction  limited  PSF.  The  Strehl  ratio  yields  a  single  number  that  indicates  the  general 
performance  of  die  adaptive-optics  system.  Also,  when  die  Strehl  ratio  is  plotted  against 
object  angle,  it  is  easily  seen  how  performance  degrades  with  increasing  d. 

43  Results 

The  atmospheric  profile  used  in  this  analysis  is  the  same  Hufnagel- Valley  atmospheric 

profile  with  a  54  mph  upper  atmospheric  wind  (HV-54)  used  in  previous  chapters.  Individual 

points  of  the  atmospheric  OTF  were  computed  using  the  method  discussed  in  Chapter  II 

and  multiplied  by  the  diffraction  limited  OTF  for  a  particular  aperture  to  obtain  an  array  of 

individual  points  for  the  overall  system  OTF.  The  resulting  array  of  OTF  values  is  inverse 

Fourier  transformed  using  standard  FFT  techniques  to  obtain  a  2-D  array  of  PSF  values. 

— * 

Figure  18  is  a  plot  of  die  PSF  for  the  axes  parallel  and  perpendicular  to  i?  for  a  1  meter 
aperture.  Figure  18  includes  the  diffraction  limited  PSF  and  the  angle  dependent  PSF  for 
object  angles  ranging  from  d  =  1  to  60o.  Recall  that  0O  is  the  isoplanatic  angle  for  the  given 
atmospheric  turbulence  profile.  The  isoplanatic  angle  for  the  HV-54,  C\  profile  is  2.4//rad. 
The  difference  between  the  PSFs  in  Figure  18  for  die  parallel  and  perpendicular  directions 
is  difficult  to  see  for  the  smaller  ranges  of  t?.  However,  for  the  larger  values  of  d,  a  definite 
widening  of  the  PSF  can  be  seen  along  the  parallel  axis.  The  widening  along  the  parallel  axis 
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will  cause  the  image  of  off-axis  stars  to  contain  an  oval  shape  with  the  major  axis  aligned 
parallel  to  t?.  Figure  19  plots  the  PSF  for  an  aperture  of  D  =  2  meters.  A  comparison  of 
Figures  18  and  19  indicates  that  the  PSF  peak  for  the  2  meter  aperture  falls  off  faster,  as  a 
function  of  the  object  angle,  than  for  the  1  meter  aperture  case.  For  example,  at  d  =  4 0o,  the 
PSF  peak  for  the  1  meter  aperture  case  is  0.027  and  the  PSF  peak  for  the  2  meter  case  is  0.01 1 . 
Recall  that  at  d  =  0,  die  PSF  is  only  limited  by  the  amplitude  OTF  attenuation.  As  a  result, 
the  PSF  peak  for  the  2  meter  aperture,  at  d  =  0,  will  be  much  larger  than  the  PSF  peak  for  the 
1  meter  aperture.  However,  as  d  increases,  the  coherence  length  of  the  atmospheric  turbulence 
begins  to  dominate  die  PSF  result.  The  dominance  of  the  atmosphere  will  tend  to  drive  die 
PSF  peak  for  the  1  and  2  meter  cases  to  the  same  values.  The  end  result,  after  normalizing 
to  the  diffraction  limited  peak,  is  the  faster  fall  off  of  the  PSF  peak  for  the  2  meter  aperture 
case.  Figure  20  is  a  plot  of  the  Strehl  ratio  versus  This  plot  again  demonstrates  the  faster 
fall  off  of  the  PSF  peak  for  the  2  meter  aperture  case  as  a  function  of  increasing  i).  Note  that 
even  at  an  object  angle  of  i)  =  0,  the  Strehl  ratio  is  down  to  0.8  which  reflects  the  attenuation 
in  the  amplitude  OTF.  Figure  21  plots  the  FWHM  of  the  PSF,  which  has  been  normalized  to 
the  diffraction  limited  PSF,  versus  4-.  The  FWHM  is  the  full  width  of  the  PSF  as  evaluated 
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at  the  point  equal  to  half  the  PSFs  maximum  value.  The  FWHM  value  has  been  normalized 
such  that  the  diffraction  limited  FWHM  equals  one.  Figure  21  clearly  shows  die  difference 
in  the  spread  of  the  PSF  for  the  axis  parallel  to  d  compared  to  the  axis  perpendicular  to  v 
for  both  the  1  and  2  meter  case.  This  figure  indicates  that  the  PSF  has  a  wider  spread  in  the 
direction  parallel  to  x).  It  is  also  seen  that  the  spread  of  the  PSF  increases  at  a  faster  rate  for 
the  2  meter  case  than  for  the  1  meter  case.  Again,  the  increasing  spread  is  an  indication  that 
at  larger  separation  angles,  atmospheric  turbulence  dominates  the  system  performance. 

4.4  Summary  and  Conclusions 

An  analysis  of  the  angle  dependent  PSF  and  Strehl  ratio  for  a  particular  atmospheric 
turbulence  profile  is  presented  in  this  chapter.  The  work  presented  in  this  chapter  represents 
the  first  time  an  analysis  of  an  angle  dependent  PSF  that  accounts  for  both  amplitude  and 
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phase  effects  in  an  adaptive-optics  system  has  been  presented.  The  results  demonstrate  die 
performance  of  the  adaptive-optics  system  with  respect  to  object  angle  out  to  a  separation  of 
6  =  66  0 .  It  is  seen  that  die  spread  of  die  PSF  is  wider  along  the  axis  parallel  to  d  than  along  die 
axis  perpendicular  to  6.  The  non-symmetric  spreading  of  the  PSF  results  in  a  non-symmetric 
image  of  an  off  axis  point  source. 


i  p  0  1 


Figure  18.  The  system  PSF  plotted  against  a  normalized  spatial  shift  in  die  image  plane  for 
a  1  meter  aperture  and  A  =  0.5  microns,  using  a  Hufnagel- Valley  turbulence 
profile  with  a  54  mph  upper  atmospheric  wind  and  ra  =  5  cm.  When  p  —  1,  the 
spatial  shift  is  equal  to  where  fi  is  the  focal  length  of  the  optics.  Individual 

plots  are  for  the  diffraction  limited  PSF  as  well  as  for  6  =  0  to  36a  in  the  upper 
plots  and  for  d  =  3  to  66 0  in  the  lower  plots.  This  figure  compares  the  PSF 
calculated  with  p  parallel  to  6  to  the  PSF  calculated  with  p  perpendicular  to  6. 
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Figure  19.  The  system  PSF  plotted  against  a  normalized  spatial  shift  in  the  image  plane  for 
a  2  meter  aperture  and  A  =  0.5  microns,  using  a  Hufnagel- Valley  turbulence 
profile  with  a  54  mph  upper  atmospheric  wind  and  r0  =  5  cm.  When  p  —  1,  the 
spatial  shift  is  equal  to  where  //  is  the  focal  length  of  the  optics.  Individual 

plots  are  for  the  diffraction  limited  PSF  as  well  as  for  t?  =  0  to  36a  in  the  upper 
plots  and  for  i?  =  3  to  66  0  in  file  lower  plots.  This  figure  compares  the  PSF 
calculated  with  p  parallel  to  t?  to  the  PSF  calculated  with  p  perpendicular  to  9. 
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Figure  20.  The  Strehl  Ratio  as  plotted  against  £■  which  represents  multiples  of  die  isoplanatic 
angle  for  A  =  0.5  microns  using  a  Hufhagel- Valley  turbulence  profile  with  a  54 
mph  upper  atmospheric  wind  and  rQ  =  5  cm.  Individual  plots  are  results  for  1 
and  2  meter  apertures. 
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Figure  21.  FWHM  of  the  PSF,  normalized  to  die  diffraction  limited  PSF,  as  plotted  against 
f-  which  represents  multiples  of  die  isoplanatic  angle.  The  PSF  is  calculated 
with  A  =  0.5  microns  using  a  Hufhagel- Valley  turbulence  profile  with  a  54  mph 
upper  atmospheric  wind  and  r0  =  5  cm.  This  figure  compares  the  FWHM  of  the 
PSF  calculated  with  p  parallel  and  perpendicular  to  t?  for  both  the  1  and  2  meter 
aperture  case. 
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V.  Improved  performance  through  a  new  wavefront  correction  algorithm 

5.1  Introduction 

The  purpose  of  this  chapter  is  to  develop  a  new  wavefront  correction  algorithm  for 
adaptive-optics  systems  that  extends  the  correctable  field -of-view  (FOV).  The  FOV  is  defined 
as  the  maximum  angle  where  good  correction  is  attainable.  Therefore,  the  FOV  is  defined  by 
the  isoplanatic  angle  as  discussed  in  Chapter  IL  The  extension  of  the  FOV  is  accomplished 
using  information  obtained  through  the  diffraction-based  propagation  model  developed  in 
Chapter  II. 

Current  correction  algorithms  involve  measuring  the  reference  wavefront  phase  in  the 
pupil  and  then  reconstructing  die  negative  of  the  measured  phase.  The  reconstructed  phase  is 
then  applied  as  a  correction  to  an  object  wavefront  In  such  a  correction  algorithm,  the  same 
correction  will  be  applied  regardless  of  the  object  angle,  i?.  The  development  in  Chapter  II 
indicates  that  information  is  known  regarding  the  difference  between  the  phase  of  the  reference 
wavefront  and  the  phase  of  the  object  wavefront  This  information  can  be  used  to  formulate  a 
better  method  of  applying  a  wavefront  correction.  The  quality  of  the  new  correction  algorithm 
will  be  judged  by  the  resulting  phase  OTF  of  fee  adaptive-optics  system.  Recall,  from 
Chapter  n,  that  the  amplitude  OTF  will  not  be  affected  by  the  phase  correction.  The  increase 
in  the  FOV  will  be  determined  by  the  size  of  (he  resulting  isoplanatic  angle  as  defined  in 
Chapter  H 

The  remainder  of  this  chapter  is  organized  as  follows.  Section  3.2  develops  the  new 
wavefront  correction  algorithm.  This  development  identifies  the  difference  between  the 
reference  and  object  phase,  and  uses  knowledge  of  the  difference  to  apply  a  better  correction 
to  the  object  wavefront  Section  3.3  tests  the  new  correction  algorithm  by  calculating  the 
OTFs  and  corresponding  isoplanatic  angles  for  two  different  atmospheric  profiles.  Finally,  a 
summary  of  this  chapter  is  presented  in  Section  5.4. 
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5  2  New  wavefront  correction  derivation 

Currently,  a  wavefront  is  corrected  by  subtracting  die  measured  phase  of  a  reference 
wavefront  from  die  phase  of  an  object  wavefront  The  subtraction  is  accomplished  in  the 
pupil  of  die  adaptive-optics  system.  The  phase  of  the  reference  and  object  wavefronts  due 
to  a  single  spatial  frequency  component  of  a  single  phase  screen  located  at  altitude  77  was 
developed  in  Chapter  n  and  is  again  presented  here: 

Pre/(»7, «,  *)  =  ->1(77,  /c)  cos(Tja)  cos(/c  •  x  +  <j>0(rj,  «)),  (74) 

and, 

P0b]{r),K,x,ti)  =  -A(i7,«)cos(f7ar)cos(j?/?  —  (k  ■  x  +  <{>0(t},k))),  (75) 

where,  a  =  /?  =  k  •  t? ,  k  is  die  wavenumber  of  the  propagating  wave,  and  3  is  the  angle 

between  the  object  and  die  reference  source  as  shown  in  Figure  1  on  page  4.  Note  that  the  only 
difference  between  the  reference  phase  and  die  object  phase  is  the  >7/?  term  contained  in  the 
argument  of  die  cos  in  Eq.  (75).  The  77/?  term  represents  additional  information  regarding  the 
phase  of  die  wavefronts  not  used  by  conventional  correction  algorithms.  If  the  atmospheric 
turbulence  was  composed  of  a  single  phase  screen  with  a  single  frequency  component,  its 
easy  to  see  that  a  wavefront  correction  algorithm  would  achieve  perfect  correction  if  a  phase 
shift  r)f3  is  added  to  PTef{v, *, %)  bef  applying  die  correction  to  P0b:(v,  x,  1?).  Since 
the  atmospheric  turbulence  is  not  a  -  .id  of  a  single  phase  screen,  a  separate  shift  would 
need  to  be  added  to  each  phase  contribution  from  every  phase  screen  to  obtain  die  perfect 
correction.  This  would  require  the  ability  to  determine  die  contribution  of  every  phase  screen  to 
the  measured  cumulative  phase  distortion.  Johnston’s  work  with  multi -conjugate  adaptive- 
optics  (MCAO)  faced  a  similar  challenge  (23,  24,  25).  However,  determining  the  altitude 
related  contributions  to  a  measured  wavefront  is  currently  a  theoretical  exercise  requiring  the 
use  of  an  array  of  laser  guide  stars.  Another  alternative  is  to  choose  a  specific  77,  given  by 
Tj0,  and  use  rj0/3  as  the  shifting  parameter  for  all  phase  screen  altitudes.  Note  that  die  tj0/3 
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shift  parameter  is  not  a  function  of  altitude,  but  is  still  a  function  of  spatial  frequency  since 

[3  =  k  ■  d 

With  an  t)0(3  shift  added  to  die  reference  phase,  the  new  reference  phase,  Prej  is  written 
as 

Prefiv, *)  =  k)  cos(7/a)  cos(i]0^  -  (k  ■  x  +  <f>0(r},  k))).  (76) 

The  residual  phase  of  the  adaptive  optical  system  is  determined  by  subtracting  the  phase  of 
die  new  reference  wavefront  from  the  wavefront  phase  of  the  object  This  residual  phase  is 
now  given  by 


APp(r),  «,ir,i9)  =  -Acos(T)a)cos(Tj/S—(K-x  +  <f>0(K))) 

+A  cos (rfct)  cos(i?0/?  -  k  ■  x  +  <£„(«))•  (77) 

Following  a  similar  analysis  as  given  in  Chapter  n,  the  resulting  phase  OTF  is  found  to  be 

Hv(pJ)  =  expj  -  87t  k2  Jj  ${k,ti)  cos2(^-)[1  -  cos(k  ■  p)\ 

x  [1  —  cos ((ri  -  t]0)k  ■  ^)]dj?dK  j.  (78) 

To  implement  this  correction  algorithm  in  an  actual  adaptive-optics  system  would 
require  the  measured  phase  to  first  be  decomposed  into  Fourier  components.  Each  frequency 
component  of  the  resulting  decomposed  wavefront  would  then  be  multiplied  by  exp(jrfo0) 
to  accomplish  the  shift  and  then  inverse  Fourier  transformed.  The  resulting  wavefront  would 
then  be  used  for  correction  of  the  object  wavefront 

53  Numerical  results 

This  section  presents  an  analysis  of  the  new  correction  algorithm  through  an  evaluation 
of  the  phase  OTF  given  in  Eq.  (78).  The  evaluation  of  Eq.  (78)  uses  continuous  turbulence 
profiles  rather  than  the  layered  profiles  as  discussed  in  Chapter  H  A  discussion  of  die  reasons 
for  using  a  continuous  profile  is  reserved  for  die  end  of  this  section. 
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Note  that  die  new  phase  OTF  expression  of  Eq.  (78)  matches  die  conventional  phase  of 
Eq.  (27)  on  page  17.  if  we  let  r)0  —  0.  Note  also,  that  if  the  atmosphere  was  composed  of  a  single 
phase  screen  at  altitude  ij0,  the  phase  OTF  of  Eq.  (78)  would  equal  1  (Le.,  ideal  correction). 
However,  when  phase  screens  are  located  away  from  tj0,  a  less  than  ideal  correction  will 
result  (i.e.  Hp  <  1).  Therefore,  for  an  initial  test,  r?0  is  chosen  equal  to  the  center  of  mass 
of  die  turbulence  profile.  For  this  test,  two  continuously  distributed  turbulence  profiles  are 
considered.  Model  SLC-N  represents  a  fit  to  die  AMOS  night  data  as  shown  in  Figure  S  on 
page  21  and  model  HV-54  represent  the  Hufnagel- Valley  turbulence  profile  calculated  with 
a  54  mph  upper  atmospheric  wind,  also  shown  in  Figure  5  on  page  21.  Model  SLC-N  has  a 
center  of  mass  equal  to  1534.25  meters  and  model  HV-54  has  a  center  of  mass  equal  to  3408.5 
meters. 

Figures  22  and  23  show  the  phase  OTFs  obtained  using  a  shift  parameter  of  tj0/3  where 
r/0  is  set  equal  to  the  center  of  mass  of  the  respective  profiles.  The  object  angles  for  these 
plots  are  the  isoplanatic  angles  defined  by  Fried  (10)  where  0O  for  model  SLC-N  is  13/xrad 
and  0o  for  model  HV-54  is  2.4/xrad.  Figures  22  and  23  show  approximately  a  1 1 .5  %  increase 
in  the  high  frequency  OTF  for  model  SLC-N  and  about  a  3  %  increase  for  model  HV-54. 
Note  also  that  for  model  HV-54,  a  significant  trade-off  consideration  must  be  made  between 
the  improvement  in  the  high  frequency  response  and  the  attenuation  in  the  low  frequency 
response  as  shown  in  Figure  23.  The  difference  in  correction  improvement  attained  by  the 
two  turbulence  profiles  is  determined  to  be  related  to  the  shape  of  the  turbulence  profile. 
By  looking  at  the  turbulence  profiles  in  Figure  5  on  page  21,  a  large  amount  of  far-field 
turbulence  in  model  HV-54  can  be  seen.  Also,  it  is  noted  that  relative  to  the  center  of  mass 
at  an  altitude  of  3408.5  meters,  there  is  a  large  amount  of  turbulence  located  a  significant 
distance  away  from  the  crater  of  mass.  The  greater  amount  of  far-field  turbulence  results  in 
a  greater  attenuation  in  the  amplitude  portion  of  the  OTF  which  is  not  corrected  for  by  the 
adaptive-optics  correction  algorithm.  The  distribution  of  turbulence  away  from  file  crater 
of  mass  results  in  a  less  effective  correction  when  the  correction  algorithm  corrects  to  the 
center  of  mass  (i.e.  r)0  =  center  of  mass).  In  contrast,  model  SLC-N  contains  less  far-field 
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turbulence  and  a  greater  proportion  of  turbulence  located  near  the  center  of  mass  at  1534.25 
meters.  These  two  properties  of  the  SLC-N  turbulence  profile  result  in  an  increase  in  the 
effectiveness  of  the  correction  algorithm  when  using  file  center  of  mass  of  the  C\  profile  for 
the  r/0/3  shift  parameter. 

The  next  step  is  to  determine  if  the  center  of  mass  is  the  best  choice  for  use  as  a  shifting 
parameter.  In  Figures  24  and  25  the  p  =  1  meter  value  of  the  phase  OTF  is  plotted  as  a 
function  of  tjq  .  These  figures  show  that  a  slightly  better  correction  is  achieved  by  shifting 
downward  in  altitude  from  the  center  of  mass. 

Figures  22  through  25  give  a  good  qualitative  look  at  the  FOV  improvement  obtained 
through  the  new  correction  algorithm.  A  more  quantitative  analysis  of  the  FOV  is  presented  in 
Figures  26  and  27.  The  FOV  of  the  adaptive-optics  system  is  defined  by  the  isoplanadc  angle. 
The  isoplanadc  angle  angle  is  defined  by  Fried  (10)  as  the  angle  where  the  high  frequency 
limit  of  the  phase  OTF  has  been  attenuated  to  a  factor  of  e~l  of  its  value  at  i?  =  0.  Using 
this  definition  of  the  FOV,  the  SLC-N  turbulence  profile  has  an  FOV  of  19.9/irad  using  the 
conventional  correction  algorithm.  The  FOV  using  file  new  correction  algorithm  with  a  shift 
corresponding  to  r/0  =  1000  meters  is  24.33/irad  or  an  increase  of  22.2%  in  the  FOV.  For  the 
HV-54  turbulence  profile,  the  increase  in  FOV  is  more  modest  at  3.6/xrad  for  file  conventional 
correction  algorithm  and  3.9/xrad  for  the  new  correction  algorithm  using  tj0  =  2000  meters. 
Use  of  the  new  correction  algorithm  results  in  an  8.3%  increase  in  the  FOV  for  model  HV-54. 
The  reasons  for  the  difference  in  the  improvement  between  the  two  turbulence  profiles  are  file 
same  as  discussed  above  for  the  high  frequency  OTF  improvement 

During  the  correction  algorithm  analysis,  tests  were  made  using  several  combinations 
of  finite  numbers  of  phase  screens  as  well  as  a  continuous  turbulence  profile.  Figure  28  shows 
the  phase  OTFs  using  a  shift  parameter  of  7?0/3  where  t]0  is  set  equal  to  the  center  of  mass 
for  the  atmospheric  model  SLC-N.  The  two  curves  plotted  in  Figure  28  are  a  comparison 
of  the  OTF  calculated  using  a  4-layered  atmospheric  model  and  file  OTF  calculated  using 
a  continuous  atmospheric  profile.  Similarly,  Figure  29  plots  the  p  =  1  meter  value  of  the 
phase  OTF  as  a  function  of  rj0.  For  r/0  =  0  (Le.  the  conventional  correction  algorithm). 
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Figure  22.  Comparison  between  the  phase  OTFs  for  die  new  correction  algorithm  and  the 
conventional  correction  algorithm.  The  OTFs  are  calculated  with  a  A  =  0.5/rni, 
an  AMOS  night  turbulence  profile,  t?  =  13.0/zrad,  and  a  Von  Karman  refractive 
index  power  spectrum  with  Lm  =  1mm  and  L0  =  5m.  The  amplitude  OTF  is  not 
affected  with  the  new  correction  algorithm  but  is  shown  here  for  reference. 

the  calculation  methods  produce  a  2%  difference  in  OTF  value.  However,  at  the  peak  of 
rj0  =  lOOOmeters,  the  calculation  methods  produce  a  7%  difference  in  OTF  value.  Results  for 
the  conventional  algorithm  calculations  were  also  verified  in  Chapter  n.  The  results  of  this 
analysis  are  a  reminder  that  a  layered  atmospheric  model  is  only  an  approximation.  All  real 
turbulence  profiles  are  continuously  distributed  with  respect  to  altitude.  Many  measurements 
of  atmospheric  turbulence  have  been  made  (3,  29,  37,  41,  50)  and  are  commonly  used  to 
justify  the  use  of  a  layered  atmosphere  for  all  applications.  The  measured  turbulence  data  only 
indicates  that  some  turbulence  profiles  may  be  modeled  as  discrete  layers  better  than  others. 
In  Chapter  n  and  Chapter  HI,  it  was  determined  that  a  layered  atmospheric  profile  was  a  good 
approximation  for  use  in  OTF  and  SNR  calculations.  However,  before  these  conclusions  were 
reached,  the  results  using  a  layered  profile  were  compared  to  a  continuous  atmospheric  profile. 
It  is  advised  that  for  any  new  correction  algorithm,  theory  should  be  developed  to  account  for 
the  continuous  atmospheric  profile  before  simplifications  to  a  layered  atmosphere  are  made. 
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Figure  23.  Comparison  between  the  phase  OTFs  for  the  new  correction  algorithm  and  the 
conventional  correction  algorithm.  The  OTFs  are  calculated  with  a  A  =  0.5/rm, 
a  Hufhagel- Valley  turbulence  profile  with  a  54  mph  upper  altitude  wind  (21), 
t?  =  2.4/rrad,  and  a  Von  Karman  refractive  index  power  spectrum  with  Lm  = 
1mm  and  L0  =  5m.  The  amplitude  OTF  is  not  affected  with  die  new  correction 
algorithm  but  is  shown  here  for  reference. 


68 


Figure  24.  The  phase  OTF  plotted  against  shift  parameter  rj0  for  model  SLC-N.  The  com¬ 
puted  value  at  each  altitude  represents  the  phase  OTF  at  p  —  lm  with  a 
A  =  0.5/im,  an  AMOS  night  turbulence  profile,  d  =  13.0/xrad,  and  a  Von 
Karman  refractive  index  power  spectrum  with  Lm  =  1mm  and  L0  =  5m. 


ho  (m) 

Figure  25.  The  phase  OTF  plotted  against  shift  parameter  tj0  for  model  HV-54.  The 
computed  value  at  each  altitude  represents  the  phase  OTF  at  p  =  lm  with  a 
A  =  0.5/im,  a  Hufnagel- Valley  turbulence  profile  with  a  54  mph  upper  altitude 
wind  (21),  i?  =  2.4/trad,  and  a  Von  Karman  refractive  index  power  spectrum 
with  Lm  =  1mm  and  L0  =  5m. 
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Figure  26.  Isoplanatic  angle  for  atmospheric  model  SLC-N.  The  computed  value  at  each 
angular  separation  represents  the  OTF  at  p  =  oom  with  a  A  =  0.5/zm,  an  AMOS 
night  turbulence  profile,  r)0  =  1000.0m,  and  a  Von  Karman  refractive  index 
power  spectrum  with  Lm  =  1mm  and  L0  =  5m. 


Figure  27.  Isoplanatic  angle  for  atmospheric  model  HV-54.  The  computed  value  at  each  an¬ 
gular  separation  represents  file  OTF  at  p  =  oom  with  a  A  =  0.5/xm,  a  Hufnagel- 
Valley  turbulence  profile  with  a  54  mph  upper  altitude  wind  (21),  tj0  =  2000.0m, 
and  a  Von  Karman  refractive  index  power  spectrum  with  Lm  =  1mm  and  L0  = 
5m. 
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Figure  28.  Comparison  with  layered  atmospheric  model  for  new  correction  algorithm  OTF 
calculation  using  atmospheric  model  SLC-N.  The  OTFs  are  calculated  with  a 
A  =  0.5/zm,  an  AMOS  night  turbulence  profile,  =  13.0/zrad,  and  a  Von 
Karman  refractive  index  power  spectrum  with  Lm  =  1mm  and  L0  =  5m.  The 
layered  atmospheric  calculation  uses  file  4-layered  atmospheric  model  discussed 
in  Chapter  II 
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Figure  29.  Comparison  with  layered  atmospheric  model  for  new  correction  algorithm  OTF 
calculation  plotted  against  shift  parameter  t]0  for  model  SLC-N.  The  computed 
value  at  each  altitude  represents  the  phase  OTF  at  p  =  lm  with  a  A  =  0.5/xm, 
an  AMOS  night  turbulence  profile,  d  =  13.0/xrad,  and  a  Von  Karman  refractive 
index  power  spectrum  with  Lm  =  1mm  and  L0  =  5m.  The  layered  atmospheric 
calculation  uses  the  4-layered  atmospheric  model  discussed  in  Chapter  n 
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5.4  Summary  and  Conclusions 

A  new  wavefront  correction  algorithm  using  the  diffraction-based  model  developed 
in  Chapter  n  is  developed  and  presented  in  this  chapter.  This  new  correction  algorithm 
makes  use  of  additional  information  regarding  the  phase  of  the  propagating  wavefronts.  The 
improvement  in  die  correction  algorithm  is  made  through  a  spatial  frequency  dependent  shift 
parameter.  r/0/9,  where  rjQ  is  thought  of  as  an  altitude.  The  best  choice  of  rj0  is  found  to  be 
an  altitude  slightly  lower  than  the  calculated  center  of  mass  of  die  atmospheric  turbulence 
profile.  By  choosing  the  optimum  value  of  r?0,  the  correctable  FOV  increased  by  22.2%  for  the 
SCL-N  turbulence  profile  and  by  8.3%  for  die  HV-54  profile.  The  difference  in  die  correctable 
FOV  increase  is  determined  to  be  due  to  die  distribution  of  die  Cl  profiles.  A  turbulence 
profile  with  turbulence  distributed  closer  to  a  single  altitude  achieves  greater  benefit  with  the 
new  correction  algorithm  presented  in  this  chapter. 

It  is  also  determined  through  this  analysis  that  care  must  be  taken  regarding  modeling 
a  continuous  atmospheric  turbulence  profile  with  discrete  layers.  A  significant  difference  in 
the  correction  results  is  noted  between  using  a  continuous  atmospheric  turbulence  profile  and 
using  a  layered  model  ‘approximation’  of  the  turbulence  profile. 
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VI.  Conclusions  and  recommendations 


6.1  Introduction 

This  research  effort  takes  a  new  approach  to  the  analysis  of  performance  measures  of  an 
adaptive-optics  system.  In  a  typical  adaptive-optics  system,  the  object  and  reference  beacon 
are  separated  by  an  object  angle.  The  correction  in  an  adaptive-optics  system  is  accomplished 
by  subtracting  die  measured  phase  of  the  reference  beacon  from  die  phase  of  the  object  The 
non-homogeneous  nature  of  atmospheric  turbulence  results  in  performance  measures  that  are 
a  function  of  die  object  angle.  In  this  research  a  diffraction-based  propagation  model  is  used 
to  track  both  the  amplitude  and  phase  variations  in  a  propagating  wavefront  The  diffraction- 
based  model  allows  for  proper  accounting  of  the  amplitude  and  phase  contributions  to  the 
performance  of  the  adaptive-optics  system.  Several  analysis  tools  are  also  developed  which 
allowed  a  numerical  evaluation  of  die  derived  performance  expressions. 

This  chapter  presents  a  summary  of  the  contributions  of  this  research  along  with  the 
significance  of  the  contributions.  These  contributions  are  unique  and  are  seen  to  have  a  signif¬ 
icant  impact  on  the  theoretical  thinking  of  die  adaptive-optics  community.  Recommendations 
for  further  research  are  also  given. 

6.2  Contributions 

Propagation  Model:  A  diffraction-based  model  for  plane-wave  propagation  at  an 
arbitrary  angle  through  atmospheric  turbulence  is  developed  in  Chapter  H  This  model  is 
based  on  the  work  of  Lee  and  Harp  (30)  for  propagation  along  the  optic  axis  and  extends 
their  results  to  account  for  propagation  at  an  arbitrary  angle.  Results  for  propagation  at  an 
arbitrary  angle  are  required  in  adaptive-optics  analysis  since  wavefront  correction  is  based  on 
the  phase  difference  between  waves  propagating  at  different  angles.  The  diffraction-based 
propagation  model  is  essential  for  calculation  of  the  amplitude  and  phase  variations  of  a 
propagating  wave.  The  conventional  geometric  optics  propagation  model  neglects  diffraction 
effects  and  therefore  can  not  account  for  the  amplitude  variations  in  a  propagating  wave. 
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Neglecting  amplitude  variations  results  in  an  inaccurate  prediction  of  adaptive-optics  system 
performance. 

Amplitude  and  Phase  Optical  Transfer  Function  (OTF)  Derivation:  An  object  angle 
dependent,  average  amplitude  and  phase  OTF  of  an  adaptive-optics  system  is  developed  and 
presented  in  Chapter  n.  The  overall  OTF  is  a  product  of  the  amplitude  and  phase  OTFs. 
Computing  the  amplitude  and  phase  OTFs,  requires  the  previously  developed  diffraction 
based  propagation  model  to  characterize  die  amplitude  and  phase  transmittance  functions  for 
an  adaptive-optics  system.  The  transmittance  functions  are  converted  to  transfer  functions 
through  definitions  taken  from  Goodman  (14).  The  OTF  expressions  presented  in  Chapter  n 
represent  the  first  time  that  angle  dependent  amplitude  and  phase  transfer  functions  have  been 
developed  for  an  adaptive-optical  system.  The  transfer  functions  are  in  the  form  of  a  three 
dimensional  integration  over  the  two  dimensions  of  the  turbulence  spectrum  and  the  single 
dimension  of  the  propagation  path. 

Layered  Turbulence:  A  moment  matching  method  of  modeling  atmospheric  turbu¬ 
lence  through  an  arbitrary  number  of  turbulent  layers  located  at  pre-determined  altitudes  is 
developed  and  presented  in  Chapter  H  The  three  dimensional  integration  for  die  OTF  is 
computationally  intensive  and  requires  recomputation  for  every  new  atmospheric  turbulence 
profile.  The  layered  turbulence  model  presented  in  this  dissertation,  uses  a  set  of  phase 
screens  located  at  predetermined  altitudes  to  represent  file  atmospheric  turbulence  profile. 
This  atmospheric  model  using  predetermined  altitudes  is  required  to  develop  a  method  of 
rapidly  calculating  the  OTF  under  changing  atmospheric  conditions.  The  layered  turbulence 
model  using  moment  matching  phase  screens  at  predetermined  altitudes  as  presented  in  this 
dissertation  is  the  first  time  such  a  model  has  been  developed. 

OTF  Evaluation:  A  method  of  rapidly  evaluating  the  OTF  of  an  adaptive-optics 
system  which  allows  scaling  «o  different  atmospheric  conditions  is  developed  and  presented 
in  Chapter  EL  This  method  requires  file  development  of  new  theory  to  represent  the  OTF  by  a 
series  of  normalized  phase  transfer  functions  and  amplitude  correlation  functions  for  each  of  the 
layers  of  the  turbulence  model  The  normalization  is  with  respect  to  the  strength  of  turbulence. 
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thus  allowing  scaling  to  various  turbulence  conditions.  This  method  of  computation  is  applied 
to  a  typical  atmospheric  profile  to  demonstrate  the  significant  contribution  of  the  amplitude 
OTF  to  the  performance  of  an  adaptive-optics  system.  It  is  also  shown  that  the  isoplanatic 
angle  is  larger  than  predicted  by  previous  theory  using  a  geometric  optics  calculation.  This 
indicates  that  the  adaptive-optics  system  degradation  with  respect  to  object  angle  is  less  severe 
than  stated  by  previous  geometric  optics  calculations. 

Signal-to-Noise  Ratio  SNR  Derivation:  The  object  angle  dependent  SNR  of  an 
adaptive-optics  system  OTF  is  developed  and  presented  in  Chapter  DOL  The  average  OTF 
developed  in  Chapter  n  gives  no  indication  as  to  how  a  specific  realization  of  the  OTF  varies 
from  the  average  OTF.  This  information  is  gained  through  the  OTF  SNR.  Knowledge  of  the 
SNR  is  also  required  whenever  image  reconstruction  is  used  in  conjunction  with  adaptive- 
optics  to  determine  the  spatial  frequencies  available  for  the  reconstruction  process.  The  SNR 
derivation  in  this  dissertation  represents  the  first  presentation  of  an  angle  dependent  OTF  SNR 
for  an  adaptive-optics  system.  The  resulting  SNR  is  in  the  form  of  a  seven  dimensional 
integration. 

SNR  Evaluation:  A  method  of  rapidly  evaluating  the  SNR  of  an  adaptive-optics  system 
OTF  which  allows  scaling  to  different  atmospheric  as  well  as  geometric  conditions  is  developed 
and  presented  in  Chapter  HI.  The  seven  dimensional  integral  result  of  the  SNR  derivation,  is 
actually  a  four  dimensional  integral  over  a  function  of  three  dimensional  phase  and  amplitude 
correlations.  For  the  evaluation  of  the  SNR,  the  layered  atmospheric  model  is  used  to  calculate 
an  array  of  normalized  correlation  function  values  over  a  wide  range  of  object  angles  A  least 
squares  fitting  routine  is  then  developed  to  determine  approximating  functions  to  approximate 
the  three  dimensional  integration  of  the  correlation  functions.  The  approximating  functions 
are  used  to  evaluate  the  SNR  for  a  typical  atmospheric  turbulence  profile.  This  evaluation 
demonstrates  the  utility  of  the  SNR  to  adaptive-optics  image  reconstruction  by  identifying  an 
adaptive-optics  cutoff  frequency.  The  SNR  evaluation  in  this  dissertation  represents  the  first 
presentation  of  an  evaluation  of  the  angle  dependent  OTF  SNR  for  an  adaptive-optics  system. 
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Point  Spread  Function  (PSF)  and  Strehl  Ratio  Evaluation:  The  angle  dependent 
average  PSF  and  Strehl  ratio  of  an  adaptive-optics  system  is  developed  and  presented  in 
Chapter  IV.  The  PSF  is  calculated  by  first  computing  an  array  of  OTF  values  for  a  particular 
atmospheric  profile.  The  OTF  values  are  computed  using  the  average  OTF  expressions 
and  the  OTF  evaluation  method  developed  in  Chapter  H  The  OTF  array  is  then  inverse 
Fourier  transformed  using  standard  FFT  techniques.  The  non- symmetric  nature  of  the  PSF 
is  demonstrated  with  an  observed  elongation  of  the  PSF  in  the  direction  of  the  object  angle. 
The  Strehl  ratio  is  computed  as  the  ratio  of  the  peak  of  file  PSF  to  the  peak  of  the  diffraction 
limited  PSF.  This  is  the  first  presentation  of  an  angle  dependent  PSF  that  accounts  for  both 
amplitude  and  phase  effects  in  an  adaptive-optics  system. 

Correction  Algorithm:  A  new  adaptive-optics  wavefront  correction  algorithm  which 
extends  the  correctable  FOV  of  the  adaptive-optics  system  is  developed  and  presented  in 
Chapter  V.  Information  regarding  the  wavefront  phase  differences  between  the  beacon  and 
the  object  is  gained  through  the  propagt!  on  model  presented  in  Chapter  EL  This  additional 
information  is  used  to  develop  the  new  correction  algorithm.  An  analysis  is  performed  using 
two  different  atmospheric  profiles  to  demonstrate  the  utility  of  this  new  correction  algorithm. 
It  is  shown  that  for  an  atmospheric  profile  containing  only  a  small  amount  of  high  altitude 
turbulence  (the  SLC-N  model  of  Figure  5,  page  21 ),  an  increase  of  over  20%  in  the  FOV 
is  achieved.  When  the  turbulence  profile  is  more  distributed  (the  HV-57  model  of  Figure  5, 
page  21 ),  an  increase  of  approximately  8%  in  the  FOV  is  achieved. 

6.3  Recommendations  for  future  research 

The  first  recommendation  is  for  a  more  extensive  analysis  using  the  average  OTF 
developed  in  Chapter  n.  The  analysis  should  determine  how  atmospheric  profiles  affect 
the  level  of  the  amplitude  and  phase  OTF.  An  attempt  should  be  made  to  quantify  when  a 
geometric  optics  calculation  method  is  adequate  and  when  the  diffraction-based  calculation 
method  becomes  necess&i  y.  This  same  analysis  should  be  performed  on  the  SNR  calculations. 
More  research  is  necessary  to  determine  some  simplifications  that  could  be  used  for  the  SNR, 
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and  under  what  atmospheric  conditions  the  simplifications  would  apply.  For  example,  it  should 
be  determined  what  atmospheric  conditions  allow  the  amplitude  correlations  to  be  neglected. 
It  is  suspected  that  simplifications  to  die  SNR  such  as  neglecting  the  amplitude  correlations 
would  be  valid  over  a  wider  range  of  atmospheric  conditions  than  for  the  average  OTF. 

The  atmospheric  profiles  such  as  the  ones  shown  in  Figure  S,  page  21,  are  averaged 
profiles  may  not  be  a  good  representation  of  actual  profiles.  These  profiles,  are  standard 
throughout  the  community,  but  are  created  using  measured  data  processed  through  a  smoothing 
function.  Measured  profiles  may  vary  significantly  from  die  averaged  profile.  From  the 
analysis  of  Chapter  n,  it  is  known  that  the  transformation  from  turbulence  strength  to  OTF 
value  is  a  non-linear  operation.  Therefor,  standard  smoothing  techniques  for  the  turbulence 
profile  may  be  a  poor  method  of  obtaining  a  good  representation  of  the  turbulence  for  die 
purpose  of  adaptive-optics  performance  analysis.  Two  suggested  solutions  to  this  problem 
are  proposed.  First,  determine  a  better  smoothing  operation  on  the  turbulence  profile  data. 
This  would  require  determining  die  nonlinear  transformation  relation  between  the  turbulence 
profile  and  the  performance  measures.  This  relationship  could  then  be  used  to  develop  a 
properly  weighted  smoothing  function  for  the  turbulence  data.  A  second  suggestion  is  to 
simply  treat  the  turbulence  profile  as  a  random  process  with  a  given  mean  and  standard 
deviation.  The  adaptive-optics  performance  measures  could  then  be  evaluated  at  the  mean 
value  of  the  turbulence  profile  as  well  as  at  plus  and  minus  one  or  two  standard  deviations. 
This  may  yield  a  clearer  picture  of  the  expected  adaptive-optics  performance. 

The  new  correction  algorithm  presented  in  Chapter  V  appears  to  have  significant  promise 
in  extending  the  correctable  FO  V.  More  work  could  be  done  with  the  correction  algorithm  in  its 
present  form  to  determine  what  atmospheric  conditions  yield  the  best  increase  in  the  correctable 
FOV.  It  is  suspected  that  results  of  this  analysis  would  indicate  that  as  a  profile  becomes  more 
concentrated  in  a  single  layer,  performance  is  increased.  Also,  as  turbulence  shifts  toward 
lower  altitudes,  an  increase  in  performance  should  be  seen.  Finally,  die  correction  algorithm 
indicates  that  perfect  correction  could  be  obtained  if  a  received  signal  could  be  decomposed 
into  specific  altitude  contributions.  More  research  needs  to  be  done  in  this  area.  Some  help  in 
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this  solution  may  be  obtained  through  Johnston’s  (22)  work  with  multi-conjugate  adaptive- 
optics.  It  is  also  possible  that  different  frequencies  have  a  different  optimum  correction 
altitude.  Currently,  the  correction  algorithm  of  Chapter  V  adds  a  shift  to  each  frequency 
component  that  is  based  on  a  single  altitude.  Further  research  in  this  area  could  result  in 
additional  improvement  in  the  correctable  FOV. 
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Appendix  A.  Atmospheric  model 

Id  this  appendix,  a  model  of  atmospheric  turbulence  is  presented.  The  important 
concepts  include  the  idea  of  pockets  or  ‘eddies’  of  turbulence,  the  dimensions  of  these  eddies, 
and  how  these  dimensions  are  incorporated  into  a  refractive  index  power  spectrum. 

The  atmospheric  model  begins  by  considering  the  power  spectrum  of  the  refractive 
index  variations.  The  refractive  index  variations  is  the  cause  of  atmospheric  turbulence.  It  is 
common  to  think  of  the  refractive  index  variations  as  pockets  of  air  ex  ‘eddies’,  each  with  a 
characteristic  refractive  index  as  shown  in  Figure  30.  The  power  spectrum  of  the  refractive 
index  turbulence,  given  by  $(/cx,  k„,  kz),  may  be  regarded  as  a  measure  of  the  relative 
abundance  of  eddies  with  dimensions  Lx  =  2it/kx,  Ly  =  2ir/Ky,  and  Lz  =  2tt/kz  where  L 
is  the  length  of  an  eddy  and  k  is  the  wavenumber.  The  atmospheric  model  presented  in  this 
appendix  and  used  throughout  this  dissertation  is  one  with  locally  isotropic  turbulence.  This 
means  that  at  any  point,  the  distribution  of  eddy  dimensions  is  the  same  in  all  three  directions. 
Isotropic  turbulence  results  in  a  power  spectrum  given  by  $(/c)  which  is  a  function  of  a  single 
wavenumber  k  with  a  corresponding  eddy  size  of  L  =  2-k/k. 

The  refractive  index  power  spectrum  can  be  thought  of  as  divided  into  three  significant 
regions.  The  boundaries  of  these  regions  are  defined  by  the  inner  and  outer  scale  size  of  the 
turbulence.  The  scale  size  L0  =  2x/k0  is  called  the  outer  scale  of  the  turbulence  and  the  scale 
size  of  Lm  =  2ir/Km  is  called  the  inner  scale  of  the  turbulence.  The  region  between  the  outer 
scale  size  and  the  inner  scale  size  is  called  the  inertial  subrange. 

The  power  spectrum  of  the  inertial  subrange  is  determined  simply  by  the  physical  laws 
that  govern  the  breakup  of  large  eddies  into  smaller  ones  (Le„  the  science  of  turbulent  flow). 
The  power  spectrum  of  eddies  larger  than  the  outer  scale  size,  La,  will  be  affected  by  things 
like  large  scale  geographic  and  meteorological  conditions.  Turbulent  eddies  smaller  than  the 
inner  scale  size,  Lm ,  dissipate  energy  as  a  result  of  viscous  forces  with  die  result  being  that  die 
power  spectrum  beyond  this  point  falls  off  very  rapidly.  The  sizes  of  eddies  of  most  concern 
in  adaptive  optics  are  those  in  the  inertial  subrange  which  are  modeled,  through  the  work  of 
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Figure  30.  A  representation  of  pockets  or  eddies  of  turbulence  in  the  atmosphere.  Each 
pocket  can  be  thought  of  being  composed  of  a  single  refractive  index.  The  vari¬ 
ations  in  die  refractive  index  result  in  perturbations  in  a  propagating  wavefront. 


Kolmogorov  (27),  with  a  power  spectrum  defined  by 

$(K,7?)  =  0.033Cn2(^)«-11/3,  (79) 

where  rj  is  the  altitude  and  C2(r/)  is  called  the  structure  constant  of  the  refractive  index 
fluctuations  and  is  a  measure  of  the  strength  of  the  fluctuations.  The  value  of  the  structure 
constant  depends  on  local  atmospheric  conditions  and  height  above  the  ground.  Typical 
values  near  die  ground  vary  from  10-13m-2/3  for  strong  turbulence  to  10_17m~2/3  for  weak 
turbulence  with  10-15m~2/3  as  an  average  value.  This  power  spectrum  has  appropriately  been 
called  the  Kolmogorov  power  spectrum.  The  Kolmogorov  spectrum  is  very  useful,  but  care 
must  be  taken  outside  of  the  inertial  subrange.  In  particular,  note  the  pole  in  die  spectrum  at 

K  =  0. 

In  the  analysis  presented  in  this  dissertation,  it  is  necessary  to  integrate  the  power 
spectrum  over  all  values  of  k.  To  avoid  the  singularity  at  zero  a  modification  of  the  Kolmogorov 
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spectrum  called  die  Von  Karman  spectrum  is  used.  This  spectrum  is  given  by 


#(>;,*) 


0.033C„!(>,) 

- —  e  \  "m  / 

(/C2  +  K2)11/6 


(80) 


where  k0  and  nm  are  the  outer  and  inner  scale  wavenumbers,  respectively.  Not  only  is  the 
Von  Karman  spectrum  of  a  nicer  form  for  numerical  integration,  it  is  also  more  representative 
of  actual  turbulence.  Any  refractive  index  power  spectrum  should  level  off  as  k  =»  0  since 
there  can  only  be  a  finite  amount  of  air  in  the  earths  atmosphere.  A  normalized  spectrum  of 
pfjff  with  *0  =  3%  and  Km  =  “  Pitted  in  Figure  31.  These  scale  parameters  are 

used  throughout  the  body  of  this  work  with  more  details  of  how  the  scale  parameters  effect 
system  performance  given  in  Section  2.S.  Note  that  the  Von  Karman  spectrum  with  an  outer 
scale  wavenumber  of  0,  and  an  inner  scale  wavenumber  set  at  infinity,  is  equivalent  to  the 
Kolmogorov  power  spectrum. 


Figure  31.  The  normalized  Von  Karman  refractive  index  power  spectrum,  ,  as  plotted 
against  wavenumber  k.  This  spectrum  is  plotted  for  an  outer  scale  size  of  3 
meters  and  an  inner  scale  size  equal  to  1  millimeter. 
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Appendix  B.  Phase  screen  model 

This  appendix  presents  the  concept  of  an  atmospheric  phase  screen  necessary  to  under¬ 
stand  die  model  of  atmospheric  propagation.  The  goal  is  to  model  a  section  of  the  atmosphere 
as  a  thin  phase  screen  such  that  a  wavefront  propagating  through  die  phase  screen  is  identical 
to  a  wavefront  passing  through  an  actual  section  of  atmosphere.  Lee  and  Harp  (30)  used  the 
method  of  phase  screens  to  calculate  the  phase  and  amplitude  perturbations  of  a  plane  wave 
and  their  results  were  validated  with  the  wave  equation  method  by  Clifford  and  Strohbehn  (4). 

As  a  wave  propagates  through  a  section  of  the  atmosphere,  it  experiences  phase  delays 
related  to  the  index  of  refraction  of  that  particular  section.  Since  the  atmosphere  is  non- 
homogeneous  with  respect  to  the  index  of  refraction,  different  portions  of  the  wave  will 
experience  different  phase  delays.  A  section  of  atmosphere  can  therefore  be  modeled  as  a 
phase-only  transmittance  function  or  phase  screen.  Let  the  transmittance  of  the  phase  screen 
be  represented  by 

t(Ti,x)  =  e-jWr>'*),  (81) 

where,  r/  is  the  altitude  of  the  phase  screen  and  A<^>(t?,  x)  is  the  phase  of  the  screen  at  a 
location  given  by  x.  As  a  wave  passes  through  this  screen,  a  phase  equal  to  A x)  will  be 
added  to  the  wavefront  In  order  to  properly  model  the  phase  screen,  A^(r/,x)  must  represent 
the  actual  phase  changes  experienced  in  propagation  through  a  section  of  turbulence.  Let  an 
individual  spatial  frequency  component  k,  of  A<f>(rf,x)  be  given  by 

A<j>(r],  x,  ic)  =  k  d h  a(r],  k)  cos(k  •  x  +  <t>0{ih  «)),  (82) 

where  k  defines  the  wavenumber  of  the  propagating  wave,  dh  is  a  differential  thickness  of 
the  phase  screen,  <2(77,  k)  and  <j>0(rh  «)  are  the  amplitude  and  initial  phase  offset  of  die  spatial 
Fourier  component  of  the  index  of  refraction  variations  having  a  spatial  frequency  k.  Within 
a  plane,  the  index-of-reff  action  fluctuations  are  isotropic  but  with  a  power  spectrum  that  is, 
in  general,  dependent  on  the  altitude  of  the  screen. 
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Let  $(77,  k)  be  the  complex  amplitude  of  the  Fourier  transform  component  of  the 
refractive  index  variations,  0(77, it)  cos(k-x+4>0(ti, «)) .  where  ^(77,  it )  =  0(77, 

A  relationship  between  the  complex  Fourier  components,  ^(77,  it),  and  the  refractive  index 
power  spectrum  is  given  by  Tatar  ski  (46)  as 


<^(771,«1)^*(772,/C2))  =  S(i£i  -  K2)F(r,i  -  772,  K1)dKldK2,  (83) 

where  F(t)i  -  772,  «i )  represents  the  two  dimensional  refractive  index  power  spectrum  and  6 
is  the  Dirac  delta  function.  Notice  that  Eq.  (83)  indicates  that  different  frequencies  within  a 
phase  screen  are  independent  of  each  other.  This  fact  is  essential  in  die  calculations  presented 
in  the  body  of  this  dissertation.  It  is  now  convenient  to  represent  the  phase  screen  of  Eq.  (82) 
as 

A<t>{rj,  x,  /c)  =  A(t),  it)  cos(«  •  x  +  ^(77,  /c)),  (84) 

where,  ^(77,  it)  =  k  d*>  a(r/,  it).  The  complex  phasor  representation  of  the  phase  screen  is 
then  given  by 


A^(t7,  x,  it)  =  A(7},it)e^Me^ 

=  A(  77,  /c)e**  (85) 

The  relation  between  the  phase  screen  variations  and  the  refractive  index  power  spectrum  is 
then  given  by 

(A(tji,ki)A*(t}2,k2))  =  (2k  dfti  ^(771,K1)2fc  dh2  $!’ (rj2it2)) 

4A:2(5(/c1  -  it2)F(rft  -  rj2,iti)ditidK2.  (86) 

This  relation  will  be  used  for  the  derivation  of  the  residual  amplitude  and  phase  correlation 
functions  derived  in  Appendix  D. 
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Appendix  C.  Derivation  of  amplitude  and  phase  perturbation  equations 

In  this  appendix,  die  amplitude  and  phase  perturbation  equations  are  derived  for  a 
wavefront  that  has  propagated  through  the  atmosphere  at  some  arbitrary  angle  as  shown  in 
Figure  32.  In  the  body  of  this  dissertation,  die  object  angle,  d,  is  defined  as  the  angle  between 
the  beacon  and  the  object  To  maintain  consistency,  the  optic  axis  is  defined  here  as  the  axis 
from  the  optics  to  the  beacon  and  die  turbulent  layer  and  observation  plane  are  perpendicular 
to  the  propagation  direction  from  the  beacon.  Any  propagating  wave  can  be  considered  to 
be  propagating  in  a  direction  determined  by  the  propagating  wavenumber  k,  defined  by  the 
individual  kx  and  ky  components  of  k.  In  our  system,  the  turbulent  layers  are  in  the  x,y 
plane,  and  the  optic  axis  defines  die  2  axis.  For  initial  ease  of  calculations,  die  'off  axis’ 
wave  from  the  object  is  designated  as  a  wave  having  an  arbitrary  kx  wavenumber  component 
This  arbitrary  kx  wavenumber  component  will  result  in  a  propagation  angle  with  respect  to 
the  optic  axis.  The  analysis  is  not  restricted  by  only  allowing  a  kx  wavenumber  component 
since  the  direction  of  propagation  in  the  x,  y  plane  will  define  the  x  axis  of  die  system.  The 
first  step  in  the  analysis  is  to  determine  the  measured  wavefront  obtained  from  a  point  source 
as  propagated  through  a  single  frequency  component  of  a  phase  screen  placed  at  an  arbitrary 
distance  away  from  the  aperture  plane  of  die  system.  A  discussion  of  a  phase  screen  is  given  in 
Appendix  B.  The  phase  screen  will  be  a  single  frequency  sinusoidal  grating  with  an  arbitrary 
orientation  in  die  x,  y  plane,  with  frequency  and  orientation  defined  by  ic,  located  at  an  altitude 
77,  and  written  as 

x,  k)  =  A(tj,  k)  cos(ic  •  x  +  (f>0( 77,  «)),  (87) 

where  A(r],  k)  and  <f>0( 77, «)  are  defined  as  in  Appendix  B.  In  the  following  analysis  we  assume 
that  ^(77,  ic)  <  1  which  is  consistent  with  the  weak  turbulence  assumption.  It  is  convenient 
to  represent  this  phase  screen  in  exponential  form  as 

A  <£(77,  x,  k)  =  +  e-H*z+Mr>,*))y  (88) 
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Lines  of  Constant  Phase 


Figure  32.  Scattering  of  a  wave  through  a  phase  screen 


The  incoming  wave  is  written  with  an  x  frequency  component  as 


Uref(r,x)  =  e-*k'r+k*x\  (89) 

where  kz  and  kx  are  related  to  the  propagation  wavenumber  k  by  k  =  yjk\  +  k*.  and  r  is  the 
distance  from  the  object  to  the  phase  screen  as  shown  in  Figure  32.  The  dependence  on  r  will 
eventually  be  discarded.  However,  plane-wave  propagation  is  assumed.  Therefore,  r  must 
be  sufficiently  far  away  to  yield  a  good  plane-wave  approximation. 

After  the  wave  passes  through  file  phase  screen,  file  transmitted  field  is  given  by 


Ut{r,x,K )  = 


(90) 
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Expanding  the  term  using  a  Taylor  series,  yields 


=  i  _  jA<f>(ri,x,K)  + 


(jA<^(  17,1,/c))2  (jA(f>(r),  x,  k))3 


2! 


3! 


+ 


(91) 


Assuming  that  A <^(77,  x,  k)  C  2tt,  as  implied  by  A(  17, k)  C  1  the  higher  order  terms  in 
Eq.  (91)  can  be  dropped.  Rewriting  Eq.  (90)  gives 


Ut{x,  r,  k)  = 


e  —«-)(!  -  jA<f>(r),  x, «)) 


_  c-J(fcir+fcxI)[l  -  jA(rj,  k)  cos (k  •  x  +  «))] 

=  j  ^(*1’ K)  ^cl(X  i+4>0(ri,X))  |  c-j(>t-g+»o(»?,it))^j  (92) 


This  equation  states  that  die  effect  of  passing  die  field  UTeJ  through  die  phase  screen  is 
that  there  are  now  three  plane  waves  -  one  traveling  in  the  original  direction,  and  two  scattered 
waves  traveling  at  directions  determined  by  die  k  frequency  component  of  the  phase  screen. 
In  order  to  propagate  Ut  to  die  observation  plane,  at  a  distance  r  +  tj  from  die  object,  it  is 
necessary  to  determine  the  2  components  of  the  wave  numbers  of  the  scattered  waves.  These 
two  components  are  given  by 


k'z  =  yjk2  —  (kx  +  k  cos  9)2  —  (k  sin  0)2,  (93) 

and, 

k"  =  \Jk2  —  (kx  —  k  cos  O')2  —  (/esin0)2,  (94) 

where  0  is  die  angle  that  k  makes  with  the  x  axis.  Ut  is  now  propagated  to  die  observation 
plane,  resulting  in 
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Ut(r  _|_  f  ?  £)  =  e-i(M»-+*))+***)  _  j  *)  (c-j((k*r+kxx)+k',ri-(il-g+40(v,it)))'} 

2 

j  A{ri,  k)  ( c-j((k,T+kxx)+k','T)+m-2+<j>o(Ti,il )))^ 

2 

2 

j  Mjh  £)  ^+^o(»i.<i)))^j  (95) 

2 

Since  the  variations  in  die  phase  and  amplitude  are  of  interest,  the  constant  phase  term  e~ikxT 
is  dropped  resulting  in 


Ut(r],  x,  k)  = 


_  e-j(k,i)+kxx)  _  j  AjVi  *)  ^-j(kxx+k'zT)-(it-g+<t>o(TI,i?)))^j 

_j  Ajlli  *)  ^e~j(kxx+k"n+(X-£+^o(v,>0))J 
=  e~iks  c~ikx1)  ^  (c-j(ktv-{Z-?+<i>0(Z)))J 


Making  substitutions  to  ease  die  algebraic  manipulation,  let 


K  •  x  4-  <f>0(ri,  ic), 
k"rj. 


Equation  97  is  substituted  into  Eq.  (96)  to  yield 
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Ut(rj,x,K)  =  e 


_  I  *>- 


_  e~ik*x  I 


=  e~jk*x 


_  e-j(c+< 6)^ 

^e~ik,v  -  jd(^Sci(^)c-i(*t*)(e-i(-*)  +  e-'<c+6>)^ 

_  jdi|LSc-i(2^)(e-i(«-»-(2f£))  +  ei((^)-c-fr))j 

(e--’*'”  -  ;A(7,  #c)e~-,(2^£)  cos  ((— 2“)  “  6))  •  (98) 


Substituting  die  expressions  for  a,  6,  and  c  from  Eq.  (97)  back  into  Eq.  (98)  gives 

Ut(v,x,K)  = 

x  cos  ^  2  - j  —  (k  •  x  +  <j>0(y,  «))^  ^ 

=  e-jV'*x+k'1>)  (l  -  jA(ri,K)e-MilT1-k') 

X  COS  ^7  2  ^  j  -  ( %•  2+  Mv,  «))jj  (99) 


Again,  drop  the  constant  phase  term  e~jk,v.  The  term  e~jkxX  represents  a  linear  phase  term 
or  die  unaberrated  wavefront  phase.  This  linear  phase  term  will  result  in  an  overall  tilt  to  die 
wavefront  The  final  expression  fra:  die  wavefront  is  written  as 


UM h  x,k)  =  1  -  jA(ri,  «)e  ^  '*  *  k,)  cos  ^7  ^**— -*^  -  (k  •  x  +  </>0(t],  k))^  . 


(100) 


Before  die  equations  become  too  lengthy,  the  variables  a  and  0  are  introduced,  where 


•-((W 


(101) 
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and. 


P  =  I -l-2~l  1  •  (102) 

To  find  die  amplitude  and  phase  perturbations.  Ut  is  written  in  terms  of  real  and 
imaginary  parts. 

Ut(T),x,it)  =  1  -  jA(Tj,K)e~lva  cos(r)P  -  (k  •  x  +  <&,(»?,£))) 

=  1  -  jA(rj,  k)[cos(t/q)  -  j  sin(»?a)]  cos {rjp  -  (k  •  x  +  Mv, «))) 

=  (1  -  A(t),  k)  sin^a)  cos (rfP  -  (k  ■  x  +  <f>0(rj , «)))]  - 

j[A( V,  *)  cos(j?q )  cos(i//?  -(«•!  +  </>0(tj,  «)))].  (103) 

The  amplitude  perturbations,  Pa,  can  now  easily  be  calculated  as  the  magnitude  of  Ut  minus 
one  ( jf/f  |  -  1).  The  magnitude  minus  one  is  considered,  because  propagation  began  with 
a  uniform  amplitude  wave  and  any  perturbations  of  this  unity  amplitude  is  file  quantity  of 
interest  The  amplitude  perturbation  equation  is  written  as: 

Pa  =  \Ut\~l 

=  [&(Ut)  +  &(Ut)}$-  1 

=  {[!  -  M7!, «)  sin(?7a)  cos(t?/3  -  (£  •  x  +  <f>0(Tj,  it)))]2 

+[-i4(»7,ic)cos(»7a)cos(77^  -  («  •  x  +  <^0(»7,  k)))]2}*  -  1.  (104) 

Since  the  assumption  is  made  that  the  maximum  amplitude  of  the  phase  screen  is  much  less 
than  unity  (te.,  A(t],  k )  <  1),  the  higher  order  term  of  [-A{r),  k)  cos{rja  cos(»7/?  -  (it  •  x  + 
it)))]*  is  neglected.  Therefore, 

Pa  =  {[1  -  A(tj,  it)  sin(r)a)  cos(i] p  ~(k  -x  +  «)))]2}*  -  1.  (105) 

By  expanding  the  square  term,  dropping  the  higher  order  term,  and  using  the  binomial  expan¬ 
sion  mi  the  remaining  square  root,  the  following  final  expression  for  the  amplitude  perturbation 


90 


is  obtained: 


Pa  =  k)  sin(»7a)  cos (r)fl  -  («  •  x  +  <f>0(r],  «))). 


(106) 


Now,  in  calculating  the  phase  perturbations,  recall  that  the  initial  plane  wave  is  defined 
by  a  uniform  phase  across  the  wavefront  All  constant  phase  terms  have  been  dropped  in  the 
final  wavefront  expression  of  Eq.  (100).  Therefore,  any  remaining  phase  is  considered  to  be 
a  phase  perturbation.  The  phase  perturbations,  Pp,  can  be  represented  as 


=  tan  1 

[9W)1 

*((/,)] 

=  tan-1 

‘  -A(  T) 

, «)  cos(r?a)  cos -  (/c  •  x  +  <f>0(rj, «)))  ] 

1  -  A(t},  k)  sin(f/a)  cos (i]/3  -  (k  ■  x  +  <f>0(ri ,  /c)))j 

(107) 


A  Taylor  series  expansion  of  tan-1(x)  for  die  case  when  x1  <  1  is  used,  where 


tan  *(x)  =  x  — 


(108) 


Again,  neglecting  die  higher  order  terms,  die  final  form  of  the  phase  perturbations  is  written 
as 


p  -  ~A C77,  *)  cos(7Q)  c°s(7/3  -{k-x  +  <t>0{r),  k))) 

P  1  -  A(tj,  k)  sin(?7a)  cos (r//3  -  (k  •  x  +  <f>0(v , «))) 

~  —A(t),  k)  cos(?/a)  cos(tj(3  —  (k  •  x  +  4>0{Vj  *)))•  (109) 

This  final  approximation  is  made  using  A(tj,  k)  <  1. 

Equations  (106)  and  (109)  form  the  final  results  of  this  appendix.  The  remainder  of 
this  appendix  determines  simplified  expressions  for  a  and  0  which  are  used  in  the  numerical 
analysis.  The  simplification  of  a  begins  by  substituting  Eqs.  (93)  and  (94)  into  Eq.  (101)  and 
using  a  binomial  expansion  on  the  square  root  functions. 
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a 


■  (BfH 

=  ^^/Ar2  —  (kx  +  kcos0)2  —  (/csinfl)2 

+  ^ifc2  -  (kx  -Kcosey  -  {Ksiney  -  y/&  -  k 2 
=  ^^/A:2  —  A:2  —  k2  —  2kxKcos  9 

+  ^^A:2  —  A:2  —  k2  +  2fcr>c  cos  6  —  \Jk2  —  W* 


1  rZi — Z?/'  i  *2  +  2kxKcos6 

=  2V*-*V~  "2-A:2 


+ 


/c2  —  2A:xk  cos  0 

A;2  —  A:2  V' 


(110) 


Now  use  the  binomial  expansion  where, 


X  X2  X3 


-  >-§• 


dll) 


This  approximation  is  very  good  for  the  case  when  x  <  1 .  Applying  the  binomial  expansion 
to  hie  square  roots  of  Eq.  (1 10),  hie  following  approximation  is  obtained: 


(112) 
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Finally,  for  small  values  of  9  we  use  the  approximation  that  kz  «  k  which  yields  the  final 
form  for  a  as 

a  =  (ir)-  <113) 

Notice  that  in  this  final  approximation,  the  kx  dependency  has  been  lost  Under  these  con¬ 
ditions,  a  for  the  reference  beacon  wave  is  equal  to  the  a  for  the  object  wave.  This  adds  a 
significant  simplification  to  the  analysis  in  the  body  of  this  dissertation.  Using  the  same  type 
of  analysis  as  used  for  she  simplification  of  a,  a  simplified  expression  for  fi  is  determined  as: 


K  -  K 


\\jk2  —  ( kx  +  /ccos  0)2  —  («sin0)2 
-  \\l&  ~  (*.  —  k  cos  0)2  —  ( k  sin  0)2 

- 1  Fiftp5) 

1  /«2  —  2fcx/ccos0\\ 

2  V  k^k2  )) 


-  1  + 


kxK 


cos  9. 


(114) 


To  this  point  the  location  of  the  off  axis  source  has  been  specified  through  the  kx  component 
of  the  propagating  wave.  A  more  convenient  notation  is  tme  that  directly  specifies  the  angle 
relative  to  the  beacon  direction.  The  relation  between  this  angle  designated  9  and  the  kx 
component  of  the  propagating  wave  is  given  by 

kx 

=  tan  t),  (115) 


or. 


(116) 
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For  small  arguments,  we  use  die  approximation,  tan-1  (ar)  «  x  to  yield 


(117) 


This  allows  /?  to  be  written  as 

/?  =  cos  0.  (118) 

— • 

The  vector  d  is  defined  to  be  a  vector  in  the  x,  y  plane  with  a  magnitude  equal  to  the  object 
angle,  t),  and  an  orientation  along  the  x  axis.  With  die  vector  definition  of  d,  ft  is  written  as 

/?  =  /?•£  (119) 
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Appendix  D.  Derivation  of  amplitude  and  phase  correlation  functions 

The  purpose  of  this  appendix  is  to  derive  the  correlation  function  for  the  residual 
phase  and  amplitude  of  a  wavefront  in  die  pupil  of  an  adaptive-optics  system.  Much  of  this 
derivation  uses  theoretical  development  given  by  Tatarski  for  the  case  of  a  wave  propagating 
through  turbulence  (46).  hi  this  derivation,  Tatarski’s  theory  is  applied  to  the  angle  dependent 
amplitude  variations  and  the  ‘corrected’  angle  dependent  phase  variations  as  discussed  in 
Appendix  C. 

The  amplitude  perturbation  equation  was  derived  in  Appendix  C  and  defined  in  complex 
phasor  notation  in  Chapter  II  as 

Pa(x,  t?)  =  //d,d«A(r,,ic)e-^sin(,a)e^.  (120) 

The  correlation  function  is  therefore  defined  by 
I a{p,*)  =  (Pa(xJ)Pam(x~P^)) 

=  (////  d^ldJ72dK1dK2A(J71,K1)e■•'(i^,'^)sin(7/1al)e■7,’,/}, 

x  i‘(j/2,«2)e-i(;?j(i?-^)  sinOfea^e^  V  (121) 

Note  that  the  only  random  quantities  are  A  and  A*.  The  expectation  is  passed  through  file 
integration  and  Eq.  (86)  in  Appendix  B  is  used  to  relate  file  complex  representation  of  the 
phase  screen  variations,  A(r/,  k),  to  file  2-D  refractive  index  power  spectrum,  F(t)X  -  7/2,  k). 
The  result  is  die  correlation  function  being  written  as 

£a(p,  $)  =  4fc2  J I y<l^id?72d«F(^i  -  ^2,«)e--,'j7’sin(^iQ)sin(772a)e-'/3^T?1-’72),  (122) 

where  the  integration  over  dr?i  and  dr?2  are  from  0  to  the  upper  level  of  turbulence  given  by  L 
and  k  =  |k|  =  ^k2x  +  k*.  Next,  a  change  of  variables  is  performed  in  order  to  represent  the 
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2-D  power  spectrum  as  a  function  of  only  one  variable.  Let 


Af,  =  T}1  -  T)2, 


(123) 


2Et/  =  Vi  +  f?2‘ 


(124) 


Note  that  the  Jacobian  for  this  transformation  is  1,  and  that  the  following  equalities  hold: 


A  Tj  +  2Et/ 


(125) 


2£t/  -  A r) 


(12  6) 


The  dAr;  d  Et/  integration  is  now  performed  over  a  region  D  in  the  (At/,  Et/)  plane  where  D 
is  a  rhombus  with  vertices  of  (0, 0),  ( L ,  §),  (0,  L),  and  (-L,  ~).  The  amplitude  correlation 
is  now  written  as 


TJ^p^d)  =  4At2  Jdiie  J j  dAT/ dET/F(AT/,  k) 

.  (  An  -f  2Et/\  .  (  2Et/ 
x  sin  - - - J  sin  - 


sin  (127) 


Use  the  identity, 


_  .  f  a  +  b\  .  fa  —  b\ 

—  2 sin  I  — - —  I  sin  I  — - —  I  =  cos(a)  —  cos(o), 


(128) 


to  show  that 


At,  +  227/  ^  gin  ^q2Et/^  At /)  =  1  (cos(aA7?)  _  cos(2aET/)).  (129) 
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Substituting  Gq.  (129)  into  Eq.  (127)  to  yield 


TJipJ)  =  2k2  jdZe-1*  *  JJd dAr?  dEi^Ai*,  /c)e-'/,(A,,)(cos(QA»?)-cos(2aE»?))  (130) 

The  integrand  is  even  with  respect  to  Ar?,  so  die  integration  over  die  region  D  is  twice  the 
value  of  die  integration  over  die  right  half  plane  of  D.  Therefore, 

r„ (/?,j?)  =  4k2  JdKe~i*p^J  2J  dAt]  dEjjF( At;, /c)e^^A^(cos(aA»7)  —  cos(2aErj)) 
4-  Jl  J  dAr)  dLi]F(AT},  /c)eJ^A,J)(cos(aA»?)  —  cos(2q£»/))|  .  (131) 


It  is  known  that  the  power  spectrum  will  only  have  significant  values  for  kAtj  <  1  (i.e., 
F(  At),  k)  — ►  0  for  kAt)  >  1)  (46).  Therefore,  in  die  region  where  0  <  /cAr/  «  1,  where 
F(At),  k)  has  a  significant  value,  <  %  <  1  and  using  a  =  £  and  fi  =  ic  -3,  the 
following  approximations  can  be  used: 


cos(aA^)  =  cos(  — ^ )  «  1, 


(132) 


ej0A>7  _  iJcosfi  pj,  J 


(133) 


The  correlation  function  is  now  written  as 


ra(p)  =  4k2  jdUe  p\^j *(1  —  cos(2aEr/))  J  dAr/ dE»/F(A»7, «) 

4-  —  cos(2aS?7))  J  dAi]  <$Et)F(At),  «)j.  (134) 


Note  that  the  dependence  on  t?  has  been  lost  through  the  loss  of  /3  in  the  approximation  of 
Eq.  (133).  The  upper  limits  on  dA??  can  be  taken  to  oo  without  significantly  changing  die 
result  Now  use  the  relation  given  by  Tatarski  that  relates  the  2-D  power  spectrum  to  the  3-D 
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power  spectrum,  4>(/c)  (46): 


<1At)F(Ati,k)  =  7T$(«). 


The  amplitude  correlation  function  is  now  written  as 


Tg(p)  =  4 k2ir  Jdice  2  dE»?(l  —  cos(2aE»7)) 

+  /  dE,(l  —  COS  (2aE,))) 

=  4k2ir  JdKe~’*'ty(it)  J  d£d(l  —  cos(2a£*j)). 


Use  the  identity. 


1  —  cos(2a)  =  2sin2(a), 


or. 


1  —  cos(2o:Ej/)  =  2sin2(aSrj), 


to  write  Eq.  (136)  as 


Ta(/i)  =  8fc27r  jdUe  J*'^4>(#c)  J  dErj sin2 (at'Erf). 


(135) 


(136) 

(137) 

(138) 

(139) 


Now  use  the  property  that  the  real  valued  correlation  function  is  equal  to  one  half  the  real 
portion  of  die  complex  valued  correlation  function,  (Le.,  Ta(p)  =  ffi{rg(/>)}  ).  This  yields 
a  final  correlation  function  of 


r a(p)  =  4fc27r  JdK  cos(k  ■  p)$(k)  J  d£r?  sin2(a£»7) 
=  4^//d,d«cos(«.^)sin2(a,). 


(140) 


Equation  (140)  represents  the  final  form  of  the  amplitude  correlation  function  used  in  the  body 
of  this  dissertation.  Note  that  at  this  point  Eq.  (140)  could  be  changed  to  polar  coordinates 
where  die  =  k  d/edd  and  die  integration  over  d 9  goes  from  0  to  r.  The  integration  over  d 8 
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can  then  be  performed  analytically  using  the  relation 


cos(a  cos(0))d0  =  7rJo(a), 


(141) 


where  J0  is  die  zeroth  order  Bessel  function  of  the  first  kind.  The  correlation  function  is  now 
written  as 


Ta(p)  =  4k2  n  J J J  dQdnd  r)  k  cos(«/9  cos(0))$(*c)  sin 2(arJ?) 
=  4k2ir2  J J  drjdn  kJ0(kp)$(k)  sin2(a»/). 


(142) 


This  is  die  exact  form  of  die  amplitude  correlation  function  given  by  Tatarski  and  also  used 
by  Lee  and  Harp  (30,  46).  Note  that  the  amplitude  correlation  function  is  not  a  function  of 
object  angle.  Recall  that  an  assumption  was  made  regarding  the  size  of  the  object  angle  in 
this  derivation.  However,  the  condition  for  the  assumption  will  always  be  meet  in  practical 
applications. 

Next,  consider  the  correlation  of  the  residual  phase  in  the  pupil  of  an  adaptive-optics 
system.  The  residual  phase  perturbation  equation  was  derived  in  Appendix  C  and  defined  in 
complex  phasor  notation  in  Chapter  n  as: 


APp(x,t?)  =  J J drjdK A(t),  n)e  cos(»/a)(l  —  ejv0).  (143) 


The  correlation  of  the  residual  phase  in  then  given  by 


W,4)  =  (Ppix^P^ix-pJ)) 

=  (  //// dj/idd2dKidK2v4(i/,, K 1)e",(*1'£)  cos(ftai)(l  -  e3Vl 01) 

x  K2)e-J^j  ^-^)  cos(»72a2)(l  -  e-***)}.  (144) 
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Again,  note  that  the  only  random  quantities  are  A  and  A*.  The  expectation  is  passed  through 
die  integration  and  use  Eq.  ( 86)  in  Appendix  B  is  used  to  yield 

£p(p,  t9)  =  4 k2  J  l&r)x&q2F{rix  -  f?2, «) 

x  cos(jjiq)cos(t?jq)(1  —  e^Vl)(\  —  (145) 

where  the  integration  over  diji  and  d r?2  are  from  0  to  the  upper  level  of  turbulence  given  by 
L  and  k  =  |  ac|  =  sJk2x  4-  The  same  change  of  variables  given  by  Eqs.  (123)  and  (124)  is 
now  used  to  yield 

r p(p,0)  =  4k2  JdKe-j*  >  JJD<i^V  dEi?F(AT?, k) 

X  cos  cos  (l  -  e****)  (l  -  e^2^)  . 

(146) 

Use  the  identity, 

.a-f-6  a  —  b. 

—  2cos( — - — )cos( — - — )  =  cos(a)  +  cos(6),  (147) 

to  show  that 

(  At?  +  2Et?\  (  2E?7-At?\  1,  „  .  ,  . 

cos  [a - - - J  cos  - - - J  =  -(cos(2aE?7)  +  cos(aA??)).  (148) 

The  residual  phase  correlation  correlation  is  now  written  as 

Trip,#)  =  2k2  JdKe-jX?  JJDdArid^riF{AV,K) 

x  (1  —  ejp  ',)(1  —  e-J/?  ^  ")(cos(2qE»7)  +  cos(aAj/)),  (149) 

where  again  the  dA 7?  dEi?  integration  is  over  the  region  D  in  the  (At?,  Et?)  plane.  Next, 
consider  the  term 

(1  _  e>f>^)(l  _  e-iP^).  (150) 
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First,  let  a  =  /3T,t)  and  b  =  Z?^2.  The  terms  of  Eq.  (150)  can  then  be  simplified  as 

(1  _  eJ(a+I>))(l  _  g-j(a-b))  _  j  _  e-j{a-b)  _  ej(a+k)  +  gj26 

=  1  +  e-'26  —  ej6(e-’°  +  e—  ja) 

=  1  +  ej2b  -  ^(2  cos  (a)).  (151) 

Therefore, 

(1  —  e-"3  n+i  ”)(1  —  e~30  " )  =  1  -|-  eJ<3A,)  —  e-^-^  (2  cos(Z?E»/)).  (157 


Substituting  Eq.  (152)  into  Eq.  (149)  yields 

F P(p,ti)  =  2k2  JdKe~jit'p  JJ  dArj  dEr/F(  Ar?,  /c) 

x  (1  +  e-,/3Ar!  -  e;^^a(2cos(i3E»?)))(cos(2aE»/)  +  cos(aA»/)).  (153) 


The  integration  over  die  region  D  is  again  taken  to  be  twice  the  value  of  the  integration  over 
the  right  half  plane  of  D.  Therefore, 

rp(p,  i?)  =  \k2  2J  dArj  dT,rjF(  At?,  k) 

x  (1  +  e20Av  —  e20~*L(2co$(/3TlT})))(cos(2a'ET})  +  cos(oAi/)) 

+ a  dA t)  dHTiF{Arj,  k ) 

x  (1  +  ej/3£kv  —  e;/3^(2cos(/?E?/)))(cos(2a£77)  +  cos(oA»7))J.(154) 


Again,  F( Ai?, «)  — >  0  for  kAtj  »  1.  Therefore,  in  the  region  where  0  <  kAtj  w  1,  the 
following  approximations  are  made: 


,  A  ,  ,k2A»7. 

cos(aATj)  =  c°s(  — — )  ~  1, 


e>y3Arj  _  e>A* IK0CO8  0  ^  j 
gi5#2  _  ej^co<ff?»l. 


(155) 
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With  hie  above  approximations,  the  correlation  function  is  now  written  as 

r p(p,d)  =  4k*  JdKe-J*  p~\JJ  -  2cos(/3£»?))(cos(2a£T?)  +  l)jfE  dAiyF(Ai?,  k) 
-f-  j dSi/(2  —  2cos(^Ei/))(cos(2aSr/)  +  1)^  6At]F(Ati, «)j.  (156) 

The  upper  limits  on  dA t?  can  be  taken  to  oc  without  significantly  changing  the  result  Now  use 
die  relation  given  by  Tatarski  that  relates  the  2-D  power  spectrum  to  the  3-D  power  spectrum, 
4>(k)  (46),  as  given  in  Eq.  (135),  to  yield 

rp(p,t))  =  8irk2  JdKe~i* ^$(/c)  J  dE»?(l  —  cos(/?Er/))(cos(2aET/)  +  1) 

=  8vk2  JJ  dT)<lKe~*x'?$(K)(l  —  cos(/3r)))(cos(2ar])  +  1).  (157) 

Use  the  relation 

cos(2  ar/)  +  1  =  2cos2(a»7),  (158) 

to  write  the  final  form  of  the  complex  representation  of  the  residual  phase  correlation  function 
as 

rp(p,  i?)  =  16 xk2  J J  d?7dKe~J*  ^$(/c)(l  —  cos(/?j/))cos2(a»/).  (159) 

Again,  use  the  property  that  the  real  valued  correlation  function  is  equal  to  one  half  die  real 
portion  of  die  complex  valued  correlation  function  to  yield  a  final  correlation  function  of 

rp(p,r?)  =  8nk2  J  J  dtjdit  cos(«  •  /7)$(k)(1  —  cos(^rj))  cos2(arj).  (160) 

This  is  die  form  of  die  residual  phase  correlation  function  used  in  die  body  of  this  dissertation. 

Note  that  the  residual  phase  correlation  can  not  be  analytically  integrated  over  d 0  due  to  the 

— * 

cos (fty)  =  cos(k  •  dr})  =  cos^ndcos (0))  term. 
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Appendix  E.  The  geometric  optics  approximation 

The  purpose  of  this  appendix  is  to  show  a  direct  comparison  between  the  average  phase 
OTF  expression,  ( Hp(p ,  i?)),  derived  using  the  diffraction  method  given  in  Chapter  n  and  the 
OTF  expression  given  by  Fried  using  the  geometric  optics  approximation  (10).  The  usefulness 
of  this  comparison  is  in  giving  confidence  that  the  analysis  of  the  OTF  as  presented  in  the 
body  of  this  dissertation  has  been  properly  performed,  and  in  identifying  die  portion  of  the 
OTF  equation  that  defines  the  diffraction  effects.  This  comparison  is  made  by  beginning  with 
the  diffraction  based  OTF  and  determining  what  assumptions  need  be  be  made  to  yield  the 
equation  given  by  Fried. 

In  developing  properties  of  atmospheric  turbulence,  Goodman  (14)  addresses  the  con¬ 
cept  of  geometric  optics  by  defining  the  ‘near  field’  region  as  die  region  where  77  < 
Note  that  this  is  a  relation  between  die  altitude  of  turbulence  and  the  wavenumber  of  the 
turbulence  for  a  given  propagating  wavenumber.  If  die  turbulence  is  close  to  the  optics,  ‘near 
field’  conditions  exist  When  17  <  we  can  also  say  that  ^  <  1  and  the  assumption 
cos2(^-)  w  1  is  used.  The  following  derivation  shows  that  the  ‘near  field’  assumption  along 
with  a  specific  refractive  index  power  spectrum  is  the  only  difference  between  the  diffraction 
calculation  for  {Hp(p,  9))  and  Fried’s  geometric  optics  OTF  calculations. 

With  the  assumption  of  cos2  ( ^ )  «  1 ,  the  phase  correlation  of  Eq.  (26)  is  written  as 


rp(p  ,t?)  =  Sxk2  JJ  $(k,r))cos(k  •  p  )[1  —  cos(tj£  •  9)}6k6i] 
=  8^2[//$(«,,)cos(«.p)d«d, 

—  JJ  $(«,  rj)  cos(«  •  p  )  cos (rjk  ■  t?)d/cdr?j . 


(161) 


The  integration  over  k  is  broken  into  polar  coordinates  where  6k  =  k  6k  69.  The  integration 
limits  are:  k  -  from  0  to  00, 9  -  from  0  to  n,  and  r?  -  from  0  to  L,  the  altitude  of  the  source. 

Next,  expand  the  dot  products  of  Eq.  (161)  noting  that  9  defines  the  orientation  of  k, 
9P  defines  the  orientation  of  p  and  the  orientation  of  d  defines  the  x  axis  (Le.,  6#  =  0). 
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Equation  (161)  is  now  written  as 


r„(p,  0„,  t?)  =  8nk2  J  J  v)  cos(k  p  cos  (9  -  0„))d0d/cd»? 

-  /;/;/>,  ,  Tj)  cos(np  cos  (8— Op))  cos(f?«t?  cos(0))d0d*cd??j(162) 


The  integration  over  9  in  the  first  integration  term  of  Eq.  (162)  is  solved  using  the  following 
property: 

/  cos(acos(0  +  b))dO  =  irJo(a)  (163) 

Jo 

where,  J0(a)  represents  the  zero  order  Bessel  function  of  the  first  kind.  Therefore, 


rL  roo  ri r  rL  roo 

II  K$(k,Ti)cOs(kPCOs(0  —  9p))d9dK&q  =  X  /  /  K$(K,T})J0{Kp)dK&T). 

Jo  Jo  Jo  Jo  Jo 

(164) 

To  evaluate  die  integration  over  0  in  the  second  integration  term  of  Eq.  (162),  first  write 
the  cosine  terms  cos(/cpcos(0  -  9P))  cos(r)itti  cos(0))  as 

cos(a  cos  c)  cos  (6  cos  d)  =  i  cos(a  cos  c  —  b  cos  d)  4-  ^  cos  (a  cos  c  +  b  cos  d),  (165) 

4 U  Li 

where  a  —  up,  b  =  TfK0,  c  =  0  —  0P,  and  d  =  0.  Note  that  d  can  be  written  as  d  =  c  +  d 
where  d  =  0P.  Now  consider  the  term  a  cos  c  +  6  cos  d. 


rL  roo 


a  cos  c  +  b  cos  d  = 


|(e* +  «-*)  + +  e“*) 
^(ejc  +  e~3C)  +  ^(ej(c+c,)  + 


oJC 


o-JO 


y(a  +  be]C  )  +  — (o  +  be  ]C) 


pjo 


z*e~JC 


ze 

~2  2 
|z|e;argi  \z\e-jdsgf 


eJC  + 


-JO 


2  2 
=  \z\  cos(c  +  argz), 


(166) 
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with  d  =  c  +  d  and  z  =  a  +  beJC>.  So, 


cos(acosc)cos(6cosd)  =  ^[cos(|z|  cos(c  +  argz))  +  cos(|2/|  cos(c  +  argz'))],  (167) 

where  z'  =  a-  bejc'.  Therefore,  the  second  integration  term  of  Eq.  (162)  can  be  written  as 

rL  roo  fir  1 

/  /  /  »7)-[cos(|2|cos(c  +  arg(2)))  +  cos(|z'|cos(c  +  arg(2')))]d0dKd77  (168) 

Jo  Jo  Jo  2 

where. 


z  =  Kp  +  (j7Kt?)ej9‘>, 

z'  =  Kp  —  (T)Kti)ei9p, 

1*1  =  Kyjp2  +  2pf)d  cos  Op  +  (r/i))2, 

\z'\  =  Kyjp 2  -  2prj0  cos  Op  +  (tj t?)2.  (169) 


The  d0  integration  of  Eq.  (168)  is  accomplished  using  Eq.  (163).  The  result  of  the  integration 
over  9  yields: 

Jo  fo  **(*’  7)  (Mk\/p2  +  2/w?t?  cos  0P  +  (i]0)2) 

+  Jq(k\Jp2  —  2pr]A  cos  0P  +  (»/t?)2)^d«d»/.  (170) 


Therefore,  die  phase  correlation  is  written  as 


Tp(p,0p,ti)  =  $ir2k2  (  f  k$(k,ti) 

Jo  Jo 


MkP ) 

L 

-  ~«/o  (^KyJ p2  +  2pTp9  cos  Op  4"  (rjd)2^ 

-  IjJo  {^\Jp2  -  2prf0  cosOp  -f  (r/t?)2^ j  d«dr; 


(171) 
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and. 


rp(O,0„t?)  = 


Sn2k2  Jq  k$(k,ti)  J0{ 0)  -  -  ^Jo(*»?tf)]d/cdr7 

noo 

k$(*c,»?)[./o(0)  -  J0(«»/^)]d«dij 

rL  roo 

=  8n2k2  /  «:$(«,  j;)[l  -  J0(«7t?)]dKd»/  (172) 

Jo  Jo 


Now  let  S(p,0„,t?)  =  rp(O,0p,tf)  -  r„(p,^,t?).and  (Hp(p,0p,  0))  =  exp[-S(p,0„t?)]. 
Therefore, 


5(p,0„,i?)  =  8 ir2k2  [  [  k$(k,ti)[1  -  J0(Krid)  -  Jq(kp) 

Jo  Jo 

+^MK\Jp2  +  2pr?i9  cos  0P  +  (rp))2) 

~^Jo(KyJp2  -  2pyd cos +  (»?r?)2)]d/cdi/  (173) 

If  die  atmospheric  refractive  index  power  spectrum  is  represented  by  the  Kolmogorov  spec¬ 
trum,  as  discussed  in  Appendix  A,  the  integration  over  k  can  be  solved  analytically.  Let  the 
power  spectrum  be  given  by 

$(k,»7)  =  0.033C*(?7)k-^.  (174) 


Now,  Eq.  (173)  is  written  as 

S{p,0p, 0)  =  8*2k20.033jQLCZ(ri)y™K-${l-  J0(«p))d« 

4-  /  k~I(1  —  J0(«^^))d« 

Jo 

+  \l  ^1  -  J0(li\jp2  +  2pT)0COS0p  -I-  (^)2)^dK 
+^jf  _  J0(KyJp 2  -  2prj0cos0p  +  (»?tf)2)^d*jdr?,  (175) 
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where  an  integration  of  /0°°  k  *  d«  has  been  added  and  subtracted  to  obtain  a  desirable  form. 
Using  die  following  integration  relation, 


r oo  g  [a*\ai  ic 

JQ  IC’*(1-  J0(aK))dK  =  2l5/3^-2 

=  1.11833  a5/3,  (176) 

the  final  expression  for  S(p,  6P,  t?)  is  given  by 

S(p,6p,ti)  =  2.91/c2  jf  C2(j/)[pS  +  (v#)*  ~  ^(p2  +  2p»?t? cos 0p  +  (r?i?)2)* 

-i(p2-2p7?t?cos0p  +  (r/i?)2)«jdr/.  (177) 

By  referring  to  Fried’s  paper  and  assuming  plane  wave  propagation,  it  is  seen  that  the 
same  results  have  been  obtained  through  completely  different  methods  of  analysis.  Note  that 
the  ‘near  field’  condition  of  r\  <  ^  that  allows  the  approximation  of  cos2(£)  *  1  is 
the  same  that  will  allow  the  approximation  of  sin2(^£)  «  0.  By  referring  to  Eq.  (18)  in 
Section  2.3,  it  is  seen  that  this  approximation  results  in  the  amplitude  OTF  being  equal  to  one 
(i.e„  (Ha(p  ))  =  1).  The  conclusion  of  the  analysis  in  this  appendix  is  that  imposing  die  ‘near 
field’  condition  on  the  diffraction  based  OTF  calculation,  results  in  the  geometric  optics  based 
OTF.  Also,  by  assuming  a  Kolmogorov  power  spectrum,  the  same  OTF  solution  is  achieved 
as  the  one  proposed  by  Fried. 
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